Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,874,924 T→A 100% V31V (GTA→GTT leuE ← leucine efflux protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,874,9240TA98.2% 724.5 / 2.0 223V31V (GTA→GTTleuEleucine efflux protein
Reads supporting (aligned to +/- strand):  ref base T (4/0);  new base A (113/106);  total (117/106)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 224 total aligned reads displayed.

ACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGTACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGC  >  CP009273/1874781‑1875067
                                                                                                                                               |                                                                                                                                               
aCAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACaa                                                                                                                                             >  2:1431367/1‑148 (MQ=255)
aCAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACaa                                                                                                                                             >  2:1645479/1‑148 (MQ=255)
aCAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACaa                                                                                                                                             <  1:2327571/148‑1 (MQ=255)
  aaTGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAAc                                                                                                                                           >  2:1784536/1‑148 (MQ=255)
  aaTGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAAc                                                                                                                                           <  1:2224775/148‑1 (MQ=255)
   aTGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGaacaaaca                                                                                                                                          >  2:652451/1‑148 (MQ=255)
     gTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACTCTATTTTTGAGAACAAACAgg                                                                                                                                        >  2:2268085/1‑148 (MQ=255)
     gTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAgg                                                                                                                                        >  2:16240/1‑148 (MQ=255)
     gTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAgg                                                                                                                                        >  1:1448051/1‑148 (MQ=255)
     gTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAgg                                                                                                                                        >  1:2754039/1‑148 (MQ=255)
     gTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAgg                                                                                                                                        <  1:1560195/148‑1 (MQ=255)
     gTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAgg                                                                                                                                        <  1:3327537/148‑1 (MQ=255)
     gTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCAGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAgg                                                                                                                                        <  1:1349820/148‑1 (MQ=255)
         aataatATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTAt                                                                                                                                    <  1:568604/148‑1 (MQ=255)
          ataataTCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTAtt                                                                                                                                   >  1:3164526/1‑148 (MQ=255)
           taataTCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTAttt                                                                                                                                  <  2:1908920/148‑1 (MQ=255)
            aataTCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTg                                                                                                                                 <  2:672484/148‑1 (MQ=255)
            aataTCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTg                                                                                                                                 >  2:3101747/1‑148 (MQ=255)
               atCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGcc                                                                                                                              >  2:3066061/1‑148 (MQ=255)
                tCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGccc                                                                                                                             >  2:1972579/1‑148 (MQ=255)
                 cGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCt                                                                                                                            <  1:1515832/148‑1 (MQ=255)
                   ggggggTCTTACGTAATGTCGCCACTCCCGCCCCTGCCCGAAACCTCCATACCGCATCGCCAATAACGACACCGCAGGCCGCAAGAGAACCGCCTGTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTgg                                                                                                                          <  2:2332965/148‑1 (MQ=255)
                   gggTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTgg                                                                                                                          >  2:2712312/1‑148 (MQ=255)
                       ggtCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCAcc                                                                                                                      >  2:3156064/1‑148 (MQ=255)
                           ttaattaaTGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAaca                                                                                                                  >  2:1373568/1‑148 (MQ=255)
                                   tGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATaaaa                                                                                                          <  2:3090042/148‑1 (MQ=255)
                                       gCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATgg                                                                                                      >  2:2551255/1‑148 (MQ=255)
                                        ccACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGc                                                                                                     >  1:1760600/1‑148 (MQ=255)
                                          aCTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGccc                                                                                                   <  1:1411030/148‑1 (MQ=255)
                                          aCTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGccc                                                                                                   >  1:1887841/1‑148 (MQ=255)
                                               aGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAAcc                                                                                              <  1:1186828/148‑1 (MQ=255)
                                                gCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGTACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCa                                                                                             >  1:3221753/1‑148 (MQ=255)
                                                 cccATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAg                                                                                            >  2:1899208/1‑148 (MQ=255)
                                                 cccATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAg                                                                                            >  2:2852382/1‑148 (MQ=255)
                                                   cATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGAt                                                                                          <  1:2516434/148‑1 (MQ=255)
                                                    aTGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATa                                                                                         >  2:883519/1‑148 (MQ=255)
                                                       ccAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGt                                                                                      >  1:801605/1‑148 (MQ=255)
                                                       ccAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGt                                                                                      >  1:290680/1‑148 (MQ=255)
                                                        cAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTc                                                                                     >  1:1528551/1‑148 (MQ=255)
                                                          gAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCa                                                                                   >  1:3479235/1‑148 (MQ=255)
                                                            aaCATCAATACCGCATCGCCAATAAATACCCCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGt                                                                                 <  1:484939/148‑1 (MQ=255)
                                                            aaCATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGt                                                                                 <  1:484993/148‑1 (MQ=255)
                                                             aCATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTa                                                                                >  1:513684/1‑148 (MQ=255)
                                                             aCATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTa                                                                                >  2:299801/1‑148 (MQ=255)
                                                             aCATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTa                                                                                >  1:2028168/1‑148 (MQ=255)
                                                              cATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTaa                                                                               <  1:2491390/148‑1 (MQ=255)
                                                              cATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTaa                                                                               <  2:850638/148‑1 (MQ=255)
                                                                 cAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTc                                                                            >  2:2154798/1‑148 (MQ=255)
                                                                  aaTACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCa                                                                           <  2:3428355/148‑1 (MQ=255)
                                                                    tACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGa                                                                         <  1:1275173/148‑1 (MQ=255)
                                                                     aCCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGaa                                                                        >  2:3473855/1‑148 (MQ=255)
                                                                           tCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGt                                                                  >  2:219536/1‑148 (MQ=255)
                                                                           tCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGt                                                                  >  2:2623072/1‑148 (MQ=255)
                                                                           tCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGt                                                                  <  2:1593614/148‑1 (MQ=255)
                                                                            cGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTa                                                                 <  1:1118253/148‑1 (MQ=255)
                                                                             gCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTAt                                                                <  1:2983554/148‑1 (MQ=255)
                                                                             gCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTAGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTAt                                                                >  1:1197249/1‑148 (MQ=255)
                                                                              ccAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTAtt                                                               <  2:1298679/148‑1 (MQ=255)
                                                                               cAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTc                                                              >  1:1972061/1‑148 (MQ=255)
                                                                                   aaataCACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCg                                                          >  1:2629627/1‑148 (MQ=255)
                                                                                    aataCACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGa                                                         <  2:3278367/148‑1 (MQ=255)
                                                                                      taCACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAac                                                       >  2:178806/1‑148 (MQ=255)
                                                                                        cacCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAacac                                                     <  2:1464722/148‑1 (MQ=255)
                                                                                        cacCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAacac                                                     <  1:1645479/148‑1 (MQ=255)
                                                                                        cacCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAacac                                                     >  1:1039064/1‑148 (MQ=255)
                                                                                         acCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAAcaca                                                    <  1:2790537/148‑1 (MQ=255)
                                                                                         acCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAAcaca                                                    <  2:3257983/148‑1 (MQ=255)
                                                                                         acCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAAcaca                                                    >  2:2810481/1‑148 (MQ=255)
                                                                                          ccGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACAt                                                   >  2:763163/1‑148 (MQ=255)
                                                                                              aGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAAc                                               <  2:2897698/148‑1 (MQ=255)
                                                                                                 ccGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACctc                                            >  1:1851380/1‑148 (MQ=255)
                                                                                                 ccGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACctc                                            <  2:16759/148‑1 (MQ=255)
                                                                                                     aaGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTa                                        >  1:1783367/1‑148 (MQ=255)
                                                                                                     aaGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTa                                        >  1:1377179/1‑148 (MQ=255)
                                                                                                     aaGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTa                                        <  1:1147770/148‑1 (MQ=255)
                                                                                                     aaGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTa                                        >  2:323277/1‑148 (MQ=255)
                                                                                                       gATAACCGCCTTTCATACCGCTAATGACGCTATTTTTGAGTACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAAt                                      >  1:3407040/1‑148 (MQ=255)
                                                                                                         tAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTa                                    <  1:3066061/148‑1 (MQ=255)
                                                                                                         tAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTa                                    <  1:2146165/148‑1 (MQ=255)
                                                                                                          aaCCGCCTTTCATTCCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTAt                                   <  2:1834587/148‑1 (MQ=255)
                                                                                                           aCCGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATc                                  <  1:2033105/148‑1 (MQ=255)
                                                                                                            ccGCCTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCt                                 <  2:3522304/148‑1 (MQ=255)
                                                                                                                cTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCg                             <  2:2644080/148‑1 (MQ=255)
                                                                                                                cTTTCATACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCg                             >  1:2719410/1‑148 (MQ=255)
                                                                                                                     aTACCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCAt                        <  2:3201066/148‑1 (MQ=255)
                                                                                                                       aCCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATgc                      <  2:3292230/148‑1 (MQ=255)
                                                                                                                       aCCGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATgc                      >  1:328753/1‑148 (MQ=255)
                                                                                                                        ccGCTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATgcg                     <  2:1057717/148‑1 (MQ=255)
                                                                                                                           cTACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGa                  >  2:2058366/1‑148 (MQ=255)
                                                                                                                            tACTGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGAt                 >  2:670000/1‑148 (MQ=255)
                                                                                                                              cTGACGCTATTTTTGAGAACAACCAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTa               <  2:3127131/148‑1 (MQ=255)
                                                                                                                               tGACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTaa              <  1:2858494/148‑1 (MQ=255)
                                                                                                                                gACGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTaaa             <  1:2338166/148‑1 (MQ=255)
                                                                                                                                 aCGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAg            >  2:1214792/1‑148 (MQ=255)
                                                                                                                                 aCGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAg            >  2:1214859/1‑148 (MQ=255)
                                                                                                                                  cGCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAgg           >  1:2209443/1‑148 (MQ=255)
                                                                                                                                   gCTATTTTTGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGt          >  2:289452/1‑148 (MQ=255)
                                                                                                                                        ttttGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAata     <  1:908033/148‑1 (MQ=255)
                                                                                                                                        ttttGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAata     <  1:700283/148‑1 (MQ=255)
                                                                                                                                           tGAGAACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGc  >  2:910678/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                               
ACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACCGCTACTGACGCTATTTTTGAGTACAAACAGGGTATTTGGCCCTGGCACCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGC  >  CP009273/1874781‑1875067

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