Predicted mutation
evidence seq id position mutation annotation gene description
RA AP012306 1,587,541 Δ2 bp coding (40‑41/210 nt) cspC ← stress protein, member of the CspA‑family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123061,587,5400A.91.4% 1257.0 / 139.8 315coding (42/210 nt)cspCstress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base A (16/11);  new base . (144/144);  total (160/155)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01
*AP0123061,587,5410G.91.4% 1298.1 / 139.2 315coding (41/210 nt)cspCstress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base G (16/11);  new base . (144/144);  total (160/155)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 338 total aligned reads displayed.

ACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTAT  >  AP012306/1587391‑1587687
                                                                                                                                                     ||                                                                                                                                                  
aCAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTa                                                                                                                                                     <  2:934742/150‑1 (MQ=255)
aCAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTa                                                                                                                                                     >  2:477858/1‑150 (MQ=255)
   gCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTactc                                                                                                                                                  <  2:798454/150‑4 (MQ=255)
   gCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTactc                                                                                                                                                  >  2:1373904/1‑147 (MQ=255)
    cTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACt                                                                                                                                                 >  2:1418839/1‑150 (MQ=255)
      gCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTactcgtt                                                                                                                                               <  1:3547963/150‑7 (MQ=255)
       ccGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGt                                                                                                                                              <  1:1469939/150‑1 (MQ=255)
        cGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGa                                                                                                                                           >  1:1583224/1‑150 (MQ=255)
         ggACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTg                                                                                                                                            >  2:3185491/1‑150 (MQ=255)
         ggACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGaa                                                                                                                                          <  2:2711862/150‑1 (MQ=255)
          gACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAAc                                                                                                                                         >  1:1704646/1‑150 (MQ=255)
           aCCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAAcc                                                                                                                                        >  1:1498929/1‑150 (MQ=255)
            ccTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCa                                                                                                                                       <  1:1159026/150‑1 (MQ=255)
             cTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCAc                                                                                                                                      <  1:2498974/150‑1 (MQ=255)
             cTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCAc                                                                                                                                      >  1:3855243/1‑150 (MQ=255)
              tttCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACt                                                                                                                                     <  1:2322069/150‑1 (MQ=255)
              tttCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACt                                                                                                                                     >  2:3433139/1‑150 (MQ=255)
                  tGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAAc                                                                                                                                 <  2:273126/150‑1 (MQ=255)
                     ccGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTg                                                                                                                              <  2:3599051/150‑1 (MQ=255)
                       gTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGAc                                                                                                                            <  2:378475/150‑1 (MQ=255)
                       gTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGAc                                                                                                                            <  1:893350/150‑1 (MQ=255)
                       gTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGAc                                                                                                                            >  1:3772689/1‑150 (MQ=255)
                           tGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCttt                                                                                                                        <  1:442225/150‑1 (MQ=255)
                             aaTTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTaa                                                                                                                      >  1:1015487/1‑150 (MQ=255)
                             aaTTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTaa                                                                                                                      <  1:3972417/150‑1 (MQ=255)
                                    aaCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGc                                                                                                               <  1:1812369/150‑1 (MQ=255)
                                    aaCTCAACGTTCTGACC‑TCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGcc                                                                                                              <  2:2667323/150‑1 (MQ=255)
                                     aCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGcc                                                                                                              <  2:1087960/150‑1 (MQ=255)
                                       tCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt                                                                                                            <  1:780196/150‑1 (MQ=255)
                                          aCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg                                                                                                         <  1:2663816/150‑1 (MQ=255)
                                           cGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa                                                                                                        <  1:3286393/150‑1 (MQ=255)
                                            gTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa                                                                                                       <  1:2854127/150‑1 (MQ=255)
                                             ttCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa                                                                                                      >  2:2856410/1‑150 (MQ=255)
                                                tGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt                                                                                                   >  2:3892281/1‑150 (MQ=255)
                                                tGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt                                                                                                   >  2:2564478/1‑150 (MQ=255)
                                                 gACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAtt                                                                                                  <  2:3012672/150‑1 (MQ=255)
                                                   ccTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAtt                                                                                                  >  2:1021747/1‑150 (MQ=255)
                                                   ccTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTcc                                                                                                >  2:2488869/1‑150 (MQ=255)
                                                   ccTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTcc                                                                                                >  1:4036415/1‑150 (MQ=255)
                                                    cTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTGATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTc                                                                                                 >  1:3680041/1‑150 (MQ=255)
                                                       cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAg                                                                                            <  1:2275015/150‑1 (MQ=255)
                                                        aGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGa                                                                                           >  2:3649148/1‑150 (MQ=255)
                                                        aGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGa                                                                                           >  2:3700828/1‑150 (MQ=255)
                                                          ccAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGAtt                                                                                         <  2:626064/150‑1 (MQ=255)
                                                            agagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGAtt                                                                                         >  2:303383/1‑150 (MQ=255)
                                                            agagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGt                                                                                       >  1:2900963/1‑150 (MQ=255)
                                                             gagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGtt                                                                                      >  1:4089655/1‑150 (MQ=255)
                                                             gagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGtt                                                                                      <  2:2439165/150‑1 (MQ=255)
                                                               gTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGtttt                                                                                    >  1:1368114/1‑150 (MQ=255)
                                                                   tGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTttcttc                                                                                <  2:725345/150‑1 (MQ=255)
                                                                     aaGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTttcttc                                                                                <  1:3022134/150‑1 (MQ=255)
                                                                      aGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGcc                                                                             >  1:595389/1‑150 (MQ=255)
                                                                         cATTACCCTGGATAGCGGAGAAGCGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCCCCCGc                                                                          >  2:398525/1‑150 (MQ=255)
                                                                            tACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAgg                                                                       <  2:113264/150‑1 (MQ=255)
                                                                               ccTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCAt                                                                    <  1:841249/150‑1 (MQ=255)
                                                                                cTGGATAGCGGAGAAGTGTACGAGCTCATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCGCTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATa                                                                   <  2:2344714/150‑1 (MQ=255)
                                                                                   gATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATaa                                                                  <  2:1243542/150‑1 (MQ=255)
                                                                                   gATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACa                                                                >  1:2468076/1‑150 (MQ=255)
                                                                                    aTAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACAt                                                               <  2:505301/150‑1 (MQ=255)
                                                                                      aGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAg                                                             <  1:2606542/150‑1 (MQ=255)
                                                                                      aGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAg                                                             >  1:3624010/1‑150 (MQ=255)
                                                                                                tgtACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGa                                                   >  2:176997/1‑150 (MQ=255)
                                                                                                   aCGATCACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGacac                                                <  2:158816/150‑1 (MQ=255)
                                                                                                    cGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGAcaca                                               <  2:413531/150‑1 (MQ=255)
                                                                                                     gAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACAt                                              >  2:4228524/1‑150 (MQ=255)
                                                                                                      aaCACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACAtt                                             <  2:890918/150‑1 (MQ=255)
                                                                                                       acacATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTa                                            >  2:2214134/1‑150 (MQ=255)
                                                                                                        cacaTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTAc                                           <  2:1155416/150‑1 (MQ=255)
                                                                                                           aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGc                                        <  1:1666947/150‑1 (MQ=255)
                                                                                                            tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCa                                       <  1:1174937/150‑1 (MQ=255)
                                                                                                              tttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATg                                     >  2:1492918/1‑150 (MQ=255)
                                                                                                              tttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATg                                     <  1:2809846/150‑1 (MQ=255)
                                                                                                                tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAg                                   >  2:1708132/1‑150 (MQ=255)
                                                                                                                  ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAg                                   <  1:3905865/150‑1 (MQ=255)
                                                                                                                  ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGc                                 <  1:974094/150‑1 (MQ=255)
                                                                                                                   tgcCATCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCa                                >  2:2130062/1‑150 (MQ=255)
                                                                                                                       aTCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCaa                            >  2:2617111/1‑150 (MQ=255)
                                                                                                                        tCAGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAat                           <  2:3805735/150‑1 (MQ=255)
                                                                                                                          aGCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAatat                         >  2:3357269/1‑150 (MQ=255)
                                                                                                                           gCCGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAAtata                        >  2:2104848/1‑150 (MQ=255)
                                                                                                                            ccGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATaa                       <  1:487527/150‑1 (MQ=255)
                                                                                                                            ccGGAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATaa                       <  1:41832/150‑1 (MQ=255)
                                                                                                                              ggAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAgg                     <  1:2070169/150‑1 (MQ=255)
                                                                                                                              ggAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAgg                     <  1:3148803/150‑1 (MQ=255)
                                                                                                                               gAGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGc                    >  1:643236/1‑150 (MQ=255)
                                                                                                                                aGTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCt                   >  2:4037734/1‑150 (MQ=255)
                                                                                                                                 gTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTc                  <  2:1903961/150‑1 (MQ=255)
                                                                                                                                 gTAATGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTc                  >  1:3735541/1‑150 (MQ=255)
                                                                                                                                  tAATGAATACAAAACCTTTAGACTCGTTCAACCACTTAACCTTACCTTTAAACTTTGCCATTTTAAAAAATCCTTAGATTGTTTTCCTCGCCCCCATGCATAACATAGATAAAACTGACTCATTACTGCATGAGGCACCAATATAAGGCt                   <  2:3930031/150‑1 (MQ=255)
                                                                                                                                   aaTGAAGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCgg                <  1:349726/150‑1 (MQ=255)
                                                                                                                                     tGAAGCCAAAACCTTT‑‑ACTCATTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCa              >  1:172404/1‑150 (MQ=255)
                                                                                                                                       aaGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCaga            >  2:3260536/1‑150 (MQ=255)
                                                                                                                                       aaGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCaga            >  2:1797026/1‑150 (MQ=255)
                                                                                                                                        aGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCagag           >  1:3446090/1‑150 (MQ=255)
                                                                                                                                        aGCCAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCagag           >  1:3268332/1‑150 (MQ=255)
                                                                                                                                          ccAAAACCTTT‑‑ACTCGTTGAACCACTTAACCTGGCCTTTAATCTTTCCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCCGGCATATCATAGATAAAACTGACACATTACTGCATGAGTCAGCAATATAAGGCTCGGCAGAGaa         >  2:2030166/1‑150 (MQ=255)
                                                                                                                                            aaaaCCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGc       >  2:3478246/1‑150 (MQ=255)
                                                                                                                                             aaaCCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCg      <  1:365101/150‑1 (MQ=255)
                                                                                                                                             aaaCCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCg      >  1:2392942/1‑150 (MQ=255)
                                                                                                                                             aaaCCTTT‑‑ACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCg      >  1:2147641/1‑150 (MQ=255)
                                                                                                                                                     ||                                                                                                                                                  
ACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTAT  >  AP012306/1587391‑1587687

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑