Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,499,568 C→A 88.1% D182Y (GAC→TAC)  ypdF ← Xaa‑Pro aminopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,499,5680CA88.1% 403.0 / 55.2 168D182Y (GAC→TAC) ypdFXaa‑Pro aminopeptidase
Reads supporting (aligned to +/- strand):  ref base C (7/13);  new base A (69/79);  total (76/92)
Fisher's exact test for biased strand distribution p-value = 3.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 170 total aligned reads displayed.

ACCAGCAAGGTGCGCGTCATATCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTCAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGATCAGCAATCCCACAGGCGAGGCGGATTTTC  >  CP009273/2499423‑2499710
                                                                                                                                                 |                                                                                                                                              
aCCAGCAAGGTGCGCGTCATATCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTaaa                                                                                                                                              >  2:2453072/1‑148 (MQ=255)
 ccAGCAAGGTGCGCGTCATATCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTaaaa                                                                                                                                             <  2:2615026/148‑1 (MQ=255)
 ccAGCAAGGTGCGCGTCATATCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTCaaa                                                                                                                                             <  1:1175026/148‑1 (MQ=255)
 ccAGCAAGGTGCGCGTCATATCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTCaaa                                                                                                                                             <  1:3106246/148‑1 (MQ=255)
  cAGCAAGGTGCGCGTCATATCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAaaaa                                                                                                                                            <  1:2263375/148‑1 (MQ=255)
       aGGTGCGCGTCATATCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTCAAAAGAGGc                                                                                                                                       >  1:1495704/1‑148 (MQ=255)
         gTGCGCGTCATATCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCtt                                                                                                                                     <  2:700181/148‑1 (MQ=255)
          tGCGCGTCATATCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCttt                                                                                                                                    <  1:703947/148‑1 (MQ=255)
          tGCGCGTCATATCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCttt                                                                                                                                    <  2:1263194/148‑1 (MQ=255)
          tGCGCGTCATATCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCttt                                                                                                                                    <  2:2865399/148‑1 (MQ=255)
          tGCGCGTCATATCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCttt                                                                                                                                    <  2:2776569/148‑1 (MQ=255)
           gcgcGTCATATCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCtttt                                                                                                                                   <  1:699545/148‑1 (MQ=255)
            cgcgTCATATCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCttttt                                                                                                                                  >  2:1690061/1‑148 (MQ=255)
               gTCATATCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTCAAAAGAGGCTTTTTCTg                                                                                                                               <  2:1344076/148‑1 (MQ=255)
               gTCATATCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTg                                                                                                                               <  2:3131874/148‑1 (MQ=255)
               gTCATATCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTg                                                                                                                               <  2:1035594/148‑1 (MQ=255)
                tCATATCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCCCTGGCGACAATGGTGGAAAAAGAGGCTTTTTCTgc                                                                                                                              >  1:1713971/1‑148 (MQ=255)
                tCATATCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTc                                                                                                                                 >  1:2693496/1‑145 (MQ=255)
                tCATATCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTc                                                                                                                                 <  2:2693496/145‑1 (MQ=255)
                   tatCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGcc                                                                                                                           >  1:1639865/1‑148 (MQ=255)
                    atCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGccc                                                                                                                          <  1:1323651/148‑1 (MQ=255)
                      cAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCctg                                                                                                                        >  1:2127217/1‑148 (MQ=255)
                         agCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCctgctg                                                                                                                     >  2:1628328/1‑148 (MQ=255)
                             gTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGgcgc                                                                                                                 <  2:2077596/148‑1 (MQ=255)
                              tAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCa                                                                                                                >  1:377011/1‑148 (MQ=255)
                              tAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCa                                                                                                                >  1:3062900/1‑148 (MQ=255)
                              tAGCCCTGATACAGCGCACCGAAATCGAGAGAGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCa                                                                                                                >  2:519990/1‑148 (MQ=255)
                                 cccTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATaa                                                                                                             >  1:209899/1‑148 (MQ=255)
                                 cccTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATaa                                                                                                             >  1:209905/1‑148 (MQ=255)
                                 cccTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATaa                                                                                                             >  2:735013/1‑148 (MQ=255)
                                 cccTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTGGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATaa                                                                                                             <  2:1384963/148‑1 (MQ=255)
                                  ccTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACTATCTTGTCGCTGGCTTTGCCGTGCGGCATCGACCCACGCAAGCCACTGGCGACAATGGTGTATAAAGATTCTTTTTCTTCGCCCTGCTGGCGCATaaa                                                                                                            <  2:295698/148‑1 (MQ=255)
                                  ccTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATaaa                                                                                                            <  2:1056080/148‑1 (MQ=255)
                                  ccTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATaaa                                                                                                            >  1:2502036/1‑148 (MQ=255)
                                  ccTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATaaa                                                                                                            <  2:3062838/148‑1 (MQ=255)
                                  ccTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATaaa                                                                                                            <  2:3062901/148‑1 (MQ=255)
                                    tGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAAcc                                                                                                          >  1:2930517/1‑148 (MQ=255)
                                     gATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCa                                                                                                         <  2:1498253/148‑1 (MQ=255)
                                           gcgcACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAg                                                                                                   >  1:2870467/1‑148 (MQ=255)
                                           gcgcACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAg                                                                                                   <  2:3010919/148‑1 (MQ=255)
                                            cgcACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTCAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGt                                                                                                  <  2:2337635/148‑1 (MQ=255)
                                            cgcACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGt                                                                                                  >  2:3024545/1‑148 (MQ=255)
                                                ccGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGATGGGGCATAAACCACTCCAGTTCAg                                                                                              <  1:1379091/148‑1 (MQ=255)
                                                ccGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAAGGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAg                                                                                              <  2:2112186/148‑1 (MQ=255)
                                                 cGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGc                                                                                             >  2:351961/1‑148 (MQ=255)
                                                 cGAAATCGAGAGTCACAAACTCGCCCGCTGCAACACTCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATACACCACTCCAGTTCAGc                                                                                             <  1:1014000/148‑1 (MQ=255)
                                                   aaaTCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTCAAAAGAGGTTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCgg                                                                                           >  2:1924132/1‑148 (MQ=255)
                                                   aaaTCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCgg                                                                                           <  1:1655250/148‑1 (MQ=255)
                                                    aaTCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGc                                                                                          >  1:75886/1‑148 (MQ=255)
                                                     aTCGAGAGTGACAAACTCGCCCGCTCCAACATTGTTGTCGCGGGCTTTGCCGTGCGGCAGCGCCCCCCGCCAGCNAGTGGCGACAATGGTTTAAAACGAGGCTTTTTCTGCGCCCTGCGGGCGAATAAACCACTCCAGTTCAGCGGCt                                                                                         <  1:1683186/148‑1 (MQ=255)
                                                      tCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTa                                                                                        >  2:2742708/1‑148 (MQ=255)
                                                               aCAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGcgct                                                                               <  1:2948230/148‑1 (MQ=255)
                                                                 aaaCTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGcgctcg                                                                             <  2:2453057/148‑1 (MQ=255)
                                                                  aaCTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTCAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGcgctcgc                                                                            <  1:1349755/148‑1 (MQ=255)
                                                                  aaCTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTCAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGcgctcgc                                                                            <  1:2996856/148‑1 (MQ=255)
                                                                     tCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGAAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCa                                                                         >  1:922300/1‑148 (MQ=255)
                                                                      cGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCAt                                                                        <  2:2710884/148‑1 (MQ=255)
                                                                        cccGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATcc                                                                      <  1:463433/148‑1 (MQ=255)
                                                                         ccGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATccc                                                                     <  2:1338270/148‑1 (MQ=255)
                                                                             tGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGcc                                                                 <  2:1670686/148‑1 (MQ=255)
                                                                             tGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGcc                                                                 <  1:2453072/148‑1 (MQ=255)
                                                                               caacaaTCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTg                                                               >  1:1483886/1‑148 (MQ=255)
                                                                                 acaaTCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTCAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGaa                                                             <  2:2859623/148‑1 (MQ=255)
                                                                                   aaTCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCCTAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATa                                                           >  1:391331/1‑148 (MQ=255)
                                                                                   aaTCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATa                                                           >  1:2555525/1‑148 (MQ=255)
                                                                                     tCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATaaa                                                         >  1:831399/1‑148 (MQ=255)
                                                                                      cTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAg                                                        <  1:2743441/148‑1 (MQ=255)
                                                                                      cTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACCCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCAGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAg                                                        >  1:2997530/1‑148 (MQ=255)
                                                                                       ttGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAgc                                                       <  2:677469/148‑1 (MQ=255)
                                                                                        tGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAgcg                                                      >  1:2491424/1‑148 (MQ=255)
                                                                                           cGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAgcggcg                                                   >  1:998953/1‑148 (MQ=255)
                                                                                             cTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGaa                                                 <  2:1082603/148‑1 (MQ=255)
                                                                                             cTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGaa                                                 <  1:2073428/148‑1 (MQ=255)
                                                                                             cTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGaa                                                 >  2:459530/1‑148 (MQ=255)
                                                                                                  tttGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGc                                            >  2:1644303/1‑148 (MQ=255)
                                                                                                  tttGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGc                                            >  1:1213886/1‑148 (MQ=255)
                                                                                                      ccGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTg                                        <  1:160287/148‑1 (MQ=255)
                                                                                                           cGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACcg                                   <  1:798151/148‑1 (MQ=255)
                                                                                                              cAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGa                                >  1:1572947/1‑148 (MQ=255)
                                                                                                              cAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAATTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGa                                >  1:1311111/1‑148 (MQ=255)
                                                                                                                 cgcCCCACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGATCa                             >  2:1823867/1‑148 (MQ=255)
                                                                                                                      cACGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGATCAGCAAt                        >  1:2826833/1‑148 (MQ=255)
                                                                                                                        cGCCAGCCACTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGATCAGCAATcc                      <  2:2229488/148‑1 (MQ=255)
                                                                                                                         gccagccaCTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGATCAGCAATccc                     <  1:2581548/148‑1 (MQ=255)
                                                                                                                            agccaCTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGATCAGCAATCCcaca                  <  1:1890759/148‑1 (MQ=255)
                                                                                                                            agccaCTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGATCAGCAATCCcaca                  <  1:1055511/148‑1 (MQ=255)
                                                                                                                             gccaGTGGCGACATTGGTGTAAAAAGAGGGTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGATCAGCAATCCCACag                 <  1:1399722/148‑1 (MQ=255)
                                                                                                                              ccaCTGGCGACAATGGTGTCAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGATCAGCAATCCCACagg                <  2:1410335/148‑1 (MQ=255)
                                                                                                                               caCTGGCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGATCAGCAATCCCACaggc               <  2:914446/148‑1 (MQ=255)
                                                                                                                                 cgggCGACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGATCAGCAATCCCACaggcga             <  2:597265/146‑1 (MQ=255)
                                                                                                                                 cTGGCGACAATGGTGTCAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTCTCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGATCAGCAATCCCACaggcga             >  2:1992371/1‑148 (MQ=255)
                                                                                                                                  tGGCGACAATGGTGTCAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGATCAGCAATCCCACaggcgag            <  2:3137444/148‑1 (MQ=255)
                                                                                                                                  tGGCGACAATGGTGTCAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGATCAGCAATCCCACaggcgag            <  2:1495704/148‑1 (MQ=255)
                                                                                                                                   ggCGACAATGGTGTCAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGATCAGCAATCCCACaggcgagg           <  1:2356771/148‑1 (MQ=255)
                                                                                                                                      gACAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGATCAGCAATCCCACAGGCGAGGCgg        <  2:328536/148‑1 (MQ=255)
                                                                                                                                       aCAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGATCAGCAATCCCACAGGCGAGGCGGa       >  2:520130/1‑148 (MQ=255)
                                                                                                                                        cAATGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGATCAGCAATCCCACAGGCGAGGCGGAt      >  1:850362/1‑148 (MQ=255)
                                                                                                                                          aTGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGATCAGCAATCCCACAGGCGAGGCGGAttt    >  1:997532/1‑148 (MQ=255)
                                                                                                                                          aTGGTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGATCAGCAATCCCACAGGCGAGGCGGAttt    <  1:1159065/148‑1 (MQ=255)
                                                                                                                                            ggTGTAAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGATCAGCAATCCCACAGGCGAGGCGGATTTtc  <  2:374193/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
ACCAGCAAGGTGCGCGTCATATCAGAGCAGTAGCCCTGATACAGCGCACCGAAATCGAGAGTGACAAACTCGCCCGCTGCAACAATCTTGTCGCTGGCTTTGCCGTGCGGCAGCGCCCCACGCCAGCCACTGGCGACAATGGTGTCAAAAGAGGCTTTTTCTGCGCCCTGCTGGCGCATAAACCACTCCAGTTCAGCGGCTATCTCGCGCTCGCTCATCCCCGCCTGAATAAAGCGGCGAATATGCTCTGCACCGCGATCAGCAATCCCACAGGCGAGGCGGATTTTC  >  CP009273/2499423‑2499710

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