Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,056,374 G→A G45D (GGC→GAC)  fau → putative 5‑formyltetrahydrofolate cyclo‑ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,056,3740GA96.5% 276.8 / NA 87G45D (GGC→GAC) fauputative 5‑formyltetrahydrofolate cyclo‑ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (36/48);  minor base C (1/2);  total (37/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAATTGGATGTGCGTGACGTGCTACCCCTTTCCCGATTAGACGTGCTGATCACACCGCTGGTCGCCTTTGATGAGTACGGTCAGCGCC  >  NC_000913/3056127‑3056619
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
gttcaagggttACAGCCTGCGGCGGCATCTCGGAGATGCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGa                                                                                                                                                                                                                                                    <  1:650654/251‑1 (MQ=255)
    aagggttACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGCCGCTGTATTCCTCTCTTTTGATGCCGAACTc                                                                                                                                                                                                                                                >  2:650178/1‑251 (MQ=255)
     agggttACAGCCTGCGGCGGCATCTCGGAGATTCCCTTTTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCg                                                                                                                                                                                                                                               <  1:257391/251‑1 (MQ=255)
      gggttACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGCAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGGGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCg                                                                                                                                                                                                                                               >  1:255538/1‑250 (MQ=255)
          tACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACAcc                                                                                                                                                                                                                                          >  1:216424/1‑251 (MQ=255)
            cAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCa                                                                                                                                                                                                                                        <  1:549903/251‑1 (MQ=255)
                   cggcggCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTc                                                                                                                                                                                                                                 >  1:372989/1‑251 (MQ=255)
                   cggcggCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACCACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATAATTCGGCAACGTCGTCGTGCGGTAACGCCGGAACAACAGCAGGAAATGGGTCACCAAGCCGCTACCCGGATGATGACTTATCCCCAGGTGGTGATCGCACATACGGTCGCTGTATTACTATCTTGTGATAACGAACTCGACACCCAGCCaatc                                                                                                                                                                                                                                 >  1:428972/1‑248 (MQ=255)
                     gcggcatacgcgaaaTTCCCTTCTTACCGGGNCCAGACCGTGACGCAACCTCCCGAACCCCCACCGCCATGATCCCGAAAGGAAATCCGCAAAATGTTTCGGCACCGGTGTCTTGCGTTAACGCCGGAACAACCGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGACACTctt                                                                                                                                                                                                                               <  2:366879/237‑3 (MQ=255)
                      ggtcaccccGGAGATCCCCTTCTTATCGGACACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCTACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATTACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCGCGTTTGATGACGAACTCGACACCCAGCCACTCATa                                                                                                                                                                                                                              <  2:424948/243‑1 (MQ=255)
                              cGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGTAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGaac                                                                                                                                                                                                                          <  1:440804/247‑1 (MQ=255)
                               ggAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAActc                                                                                                                                                                                                                     >  1:653104/1‑251 (MQ=255)
                                      cccTTCTTATCTGGCATCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGgcgc                                                                                                                                                                                                               <  1:519463/250‑1 (MQ=255)
                                        cttcttATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGTTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACANCCAGCCACTCATAGACCAACTCTGGCGCGcc                                                                                                                                                                                                            >  2:372386/1‑251 (MQ=255)
                                         ttcttATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCATCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCg                                                                                                                                                                                                           >  1:455545/1‑251 (MQ=255)
                                          ttttaacTGGCACCAGCCAGGACTGAACTACCAGACCTCCCACTGACATTAGCCCGACAAGAAAGCCGAAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAAAAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGcc                                                                                                                                                                                                            <  1:570941/243‑1 (MQ=255)
                                               cccgccACCAGCAATGACGCAACTACCAGAACCCCCACTGACATTACCCCGACAAGAAATCCGCAAAAGTATTCGGCAACTTCGTCGGTCGTTAACCCCGGAACAACAGCAGTAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGGGGTGATTGCACATACGGTCGCGGTATTCCTCGCTTGTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTTGCGCGCCGGTAAgc                                                                                                                                                                                                     <  2:672658/246‑1 (MQ=255)
                                                cTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTACCCCGACAAGAAATCCGCAAAATGCTTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACACGGCTGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCANTCATAGAACAACTCTGGCGCGCCGGAAAgcg                                                                                                                                                                                                    >  2:357187/1‑251 (MQ=255)
                                                 tGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAgcgc                                                                                                                                                                                                   >  1:404834/1‑251 (MQ=255)
                                                  ggCACCAGCCATGACGCAATTACCAGACCCCCCACTGACATTATCCCGACAAAAAACCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTACCCCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACGCATAGAACAACTCTGGCGCGCCGGTAAgcg                                                                                                                                                                                                    <  2:455545/249‑1 (MQ=255)
                                                  ggCAACAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTAATAGAACAACTCTGGCGCGCCGGTAcgcgcg                                                                                                                                                                                                  >  2:143794/1‑251 (MQ=255)
                                                    cACCAGCCATGACGCAACCACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAAGGATTCGGCAACGTCGTCGTGCGTTACCGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGta                                                                                                                                                                                                <  1:262087/251‑1 (MQ=255)
                                                    cACCAGCCATGACGAAACTCCCAGAACTCCCACTGACATGATCCCGAAAAGAAAACCGAAAAATGATTCGGCAACGTCGTCGGGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCACCAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGta                                                                                                                                                                                                <  2:252250/251‑1 (MQ=255)
                                                        cagcctggACGCAACTACCAGACCTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCTTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATATAACAACTCTGGCGCGCCGGTAAGCGCGTatat                                                                                                                                                                                             <  1:16212/244‑1 (MQ=255)
                                                         gCCATGACGCAACTACCAGAACTCCCACTGACATTACCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCt                                                                                                                                                                                           <  1:175009/251‑1 (MQ=255)
                                                            cATGACGCCACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGACGGTAAGCGCGTCTATCTTc                                                                                                                                                                                         >  2:438005/1‑251 (MQ=255)
                                                            cATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTc                                                                                                                                                                                         >  2:540054/1‑251 (MQ=255)
                                                             aTGACGCAACTACCAGAACTCCCACTGACATTATCCCGAAAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTcc                                                                                                                                                                                        <  1:438005/251‑1 (MQ=255)
                                                               gACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCa                                                                                                                                                                                       <  1:181497/250‑1 (MQ=255)
                                                                aCGCAACTCCCAGAACCCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGGCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGt                                                                                                                                                                                     <  1:96207/251‑1 (MQ=255)
                                                                 cGCAACTACCAGACCTCCCGCTGACATTATCCCGATAAGAAATCCGCAAGAGGATTCGGCAACGTCTGCGTGCGTTAACGCCGGAACACCACCAGGAAAGGGGTCAACAAGCCGCTACCCGGATTATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAGCTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAATCGCGTATATCTTCCAGtt                                                                                                                                                                                    <  2:615957/251‑1 (MQ=255)
                                                                   cAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGtttt                                                                                                                                                                                  >  1:212318/1‑251 (MQ=255)
                                                                    accTACCAGAACTCCCACTGACATTACCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACCCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTg                                                                                                                                                                                 <  1:403934/249‑1 (MQ=255)
                                                                          cAGAACCCCCACTGAAATTACCCCGACAAGAAATCCGCAAAATGATTCGGCAACGGCGTCGTGCGTTAACGCCGGAACAACAGAAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATccc                                                                                                                                                                           <  2:638720/251‑1 (MQ=255)
                                                                           aGAACTCCCACTGAAATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTAACCGGATGCTGACTTATCCCCCGGTGGTGATGGCACCTACGGTCGCTGTATTCCTCTATTTTGATGAAGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTACCGTTTTGCATCCCt                                                                                                                                                                          >  2:260654/1‑251 (MQ=255)
                                                                             aacccccACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCttt                                                                                                                                                                        <  1:40169/247‑1 (MQ=255)
                                                                                   aaCTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCAGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGATTTGCATCCCTTTAGTGcc                                                                                                                                                                  >  1:45060/2‑251 (MQ=255)
                                                                                     cTGACATTAGCCCTACAAGAAAGCCGAAAAAGGAGTCGGCAACTTCGTCGGGCGTTAACGCCGGACAAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTGCCAGTTTTGCATCCCTTTAGTGcc                                                                                                                                                                  <  2:636016/249‑1 (MQ=255)
                                                                                       gACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGt                                                                                                                                                               >  2:93263/1‑250 (MQ=255)
                                                                                       gACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGt                                                                                                                                                               >  1:322535/1‑250 (MQ=255)
                                                                                       aaCATTAACCCTACAAGAAACCAGCAAAATGATGCGGCAACGTCTTCGTGCTTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTa                                                                                                                                                              <  2:383525/250‑1 (MQ=255)
                                                                                             accccGACAAGAAACCCGAAAAATGATTCGGCAACGTCGTCGGGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGc                                                                                                                                                        <  2:488191/249‑1 (MQ=255)
                                                                                                ccGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGt                                                                                                                                                     <  1:290603/251‑1 (MQ=255)
                                                                                                   aCAAGAAATCCGCAAAAGGATGCGCAAACGTNGTCTGGCGTTAACTCCGGAACAACAGCAGGAATTGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATTCGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTcc                                                                                                                                                  <  2:372989/251‑1 (MQ=255)
                                                                                                   aCAAGAAACCCGCAAAATGATTCGGCAACGTCGTCGGGCGTTAACGCCGGAACAACAGCGGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTcc                                                                                                                                                  <  2:402372/251‑1 (MQ=255)
                                                                                                            ccGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTGGAGTGCCGGTAATTTGCTGTTCCTGAATTAc                                                                                                                                          <  1:639175/250‑1 (MQ=255)
                                                                                                             cGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTAcc                                                                                                                                         >  1:715664/1‑250 (MQ=255)
                                                                                                                  aaTGATTCGGCAACGTCGTCGAGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGc                                                                                                                                   >  1:712125/1‑251 (MQ=255)
                                                                                                                      aTTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCaa                                                                                                                                 >  1:726167/1‑249 (MQ=255)
                                                                                                                          ggCAACGTCGTCGTGCGTTAACGCCGGAACCACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTCTCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACACAGCCACTCATAGCACAACTCTGGCGCGCCGGGAAGCGAGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTCCCATCCGCAAAgaga                                                                                                                            >  2:650765/1‑247 (MQ=255)
                                                                                                                           gCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAAc                                                                                                                          <  1:267804/251‑1 (MQ=255)
                                                                                                                                  cgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTGTATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGtgat                                                                                                                   >  2:312555/1‑251 (MQ=255)
                                                                                                                                  cgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGAAACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGtgat                                                                                                                   >  2:74796/1‑251 (MQ=255)
                                                                                                                                   gtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGtgatg                                                                                                                  >  2:43695/1‑251 (MQ=255)
                                                                                                                                      gtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGAAACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAAc                                                                                                               <  1:74796/251‑1 (MQ=255)
                                                                                                                                       tGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAAc                                                                                                               >  1:723731/1‑250 (MQ=255)
                                                                                                                                        gCGTTACCGCGGGAAAAACAGCAGGAAAGGGGCTAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCGCTTTGGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACa                                                                                                              <  2:43380/250‑1 (MQ=255)
                                                                                                                                         cGTTAACGCCGGAACAACAGCCGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCAATCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACAATCCGCAAAGCGCACTGGTGATGAACAgg                                                                                                            >  1:570798/1‑251 (MQ=255)
                                                                                                                                             aaCGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTg                                                                                                         >  1:659422/1‑250 (MQ=255)
                                                                                                                                               cGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAg                                                                                                      <  1:568155/251‑1 (MQ=255)
                                                                                                                                                 ccgaaccaccaacagGAAATGGGTCAAAAAGCCGCTACCCGGATGATGACTTATCCCCCGTTGGTGATGGCCCATACGTTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCTCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGAt                                                                                                    <  2:393835/239‑1 (MQ=255)
                                                                                                                                                  cGGAAAAACAGCAGGAAATGGGGAAACAAGCCGCTACCCGGATGATGCCTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCGGGCGCGACGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATc                                                                                                   <  1:373568/251‑1 (MQ=255)
                                                                                                                                                   ggAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATcc                                                                                                  >  1:529028/1‑251 (MQ=255)
                                                                                                                                                      acaacaGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACGCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATg                                                                                               <  1:372386/251‑1 (MQ=255)
                                                                                                                                                               ggAAATGGGTCAACAAGCGGCTACCCGGATGCTGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAAt                                                                                      >  2:604297/1‑251 (MQ=255)
                                                                                                                                                                gAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAAtt                                                                                     <  1:93263/251‑1 (MQ=255)
                                                                                                                                                                   aTGGGTCAACAAGCCGCTACTCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATACCTCTCTTTTGCTGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCATTTAGCGCCGGTAATTTGCTGTTCCTGCATTACCATCCGCACAGCGTACTGGAGATGACCAGGTTGAAAATCCATGAGCCAAAATTgg                                                                                   >  2:485718/1‑250 (MQ=255)
                                                                                                                                                                             aaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAATTGGATGTGCGTGAc                                                                        <  1:685315/251‑1 (MQ=255)
                                                                                                                                                                                  gCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAATTGGATGTGCGTGACGTGCt                                                                   <  1:650178/251‑1 (MQ=255)
                                                                                                                                                                                      cccGGCTGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAATTGGATGTGCGTGACGTGCTAccc                                                               >  1:336240/1‑251 (MQ=255)
                                                                                                                                                                                             gatgaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTGTATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAATTGGATGTGCGTGACGTGCTACCCCTTTcc                                                         <  1:312555/250‑1 (MQ=255)
                                                                                                                                                                                               tgaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGCAGATCCATGAGCCAAAATTGGATGTGCGTGACGTGCTACCCCTTTCCCgt                                                      >  2:218952/1‑250 (MQ=255)
                                                                                                                                                                                                    tATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATTCATGAGCCAAAATTGGATGTGCGTGACGTGCTACCCCTTTCCCGATTAg                                                  <  1:398267/250‑1 (MQ=255)
                                                                                                                                                                                                     aTCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAATTGGATGTGCGTGACGTGCTACCCCTTTCCCGATTAGAc                                                >  1:203079/1‑251 (MQ=255)
                                                                                                                                                                                                     aTCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTGTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAATTGGATGTGCGTGACGTGCTACCCCTTTCCCGATTAGAc                                                <  2:216424/251‑1 (MQ=255)
                                                                                                                                                                                                     aTCCCCCGGTGGTGATGGCACATACGGTCGCTGTATGCCCGCTTTTGATGCCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAATTGGATGTGCGTGACGTGCTACCCCTTTCCCGATTAGAc                                                <  2:129563/250‑1 (MQ=255)
                                                                                                                                                                                                       cccccGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAATTGGATGTGCGTGACGTGCTACCCCTTTCCCGATTAGACg                                               <  1:473741/250‑1 (MQ=255)
                                                                                                                                                                                                         cacGGTGGTGAGGGCAAATACGGTCGCTGTATTCCTCTCTTTTGATGAAGAACTCGACACCCAGCCTCTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCGGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAATTGGATGTGCGTGACGTGCTACCCCTTTCCCGATTAGACGTGc                                            <  2:397886/249‑1 (MQ=255)
                                                                                                                                                                                                            gatggttATGGCACCTACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCTTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTCCCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGTGCCAAAATTGGATGTGCGTGTCGTGCTACCCCTTTCCCGATTAGACGTGCTGa                                         <  2:712125/249‑1 (MQ=255)
                                                                                                                                                                                                               ggtgATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAATTGGATGTGCGTGACGTGCTACCCCTTTCCCGATTAGACGTGCTGATca                                      <  1:451837/251‑1 (MQ=255)
                                                                                                                                                                                                                  gATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAATTGGATGTGCGTGACGTGCTACCCCTTTCCCGATTAGACGTGCTGATCacac                                   >  1:469504/1‑251 (MQ=255)
                                                                                                                                                                                                                      gaccaTACGGCCGCTGTATTCCTCTCTTTTGACGCTGAACTCGACACCCAGCCCCTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGGAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAATTGGATGTGCGTGACGTGCTACCCCTTTCCCGATTAGACGTGCTGATCACACCGCt                               <  2:300450/248‑1 (MQ=255)
                                                                                                                                                                                                                        acaTACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATACCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCCTCCGCAAAGCGAACTGGTGATGAACAGGCTGAAGATCCATGAGCCAAAATTGGATGTGCGTGACGTGCTACCCCTTTCCCGATTAGCCGTGCTGCTAACACCGCTgg                             >  2:580867/1‑251 (MQ=255)
                                                                                                                                                                                                                         caTACGGTCGCTGTATTCCCCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAATTGGATGTGCGTGACGTGCTACCCCTTTCCCGATTAGACGTGCTGATCACACCGCTGGt                            <  2:212318/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cGCTGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAATTGGATGTGCGTGACGTGCTACCCCTTTCCCGATTAGACGTGCTGATCACACCGCTGGTCGCCttt                     <  2:723731/250‑1 (MQ=255)
                                                                                                                                                                                                                                    tGTATTCCTCTCTTTTGATGACGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAATTGGATGTGCGTGACGTGCTACCCCTTTCCCGATTAGACGTGCTGATCACACCGCTGGTCGCCTTtgatg                 >  2:223511/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aTTCCTCTCTTTTGATGACGAAATCTACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCGTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAATTGGATGTTCGTGACGTGCTACCCCTTTCCCGATTAGACGTGCTGATCACGCCGCTGGTCGCCTTTGATGAGt              <  2:481856/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   tgatgaCGCACTCGACACCCAGCCACTCATAGAAAAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGGTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAGGATCCATGAGCCAAAATTGGATGGGCGTGACGTGCTACCCCTTTCCCGATTAGACGTGCTGATCACACCGCTGGTCGCCTTTGATGAGTACGGTCCGCGcc  >  2:23559/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
GTTCAAGGGTTACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAATTGGATGTGCGTGACGTGCTACCCCTTTCCCGATTAGACGTGCTGATCACACCGCTGGTCGCCTTTGATGAGTACGGTCAGCGCC  >  NC_000913/3056127‑3056619

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: