Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,571,658 C→A 25.8% A73A (GCG→GCT yddB ← putative TonB‑dependent outer membrane receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,571,6580CA25.8% 727.9 / 192.7 325A73A (GCG→GCTyddBputative TonB‑dependent outer membrane receptor
Reads supporting (aligned to +/- strand):  ref base C (135/106);  new base A (46/38);  total (181/144)
Fisher's exact test for biased strand distribution p-value = 8.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 328 total aligned reads displayed.

GACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTT  >  CP009273/1571509‑1571807
                                                                                                                                                     |                                                                                                                                                     
gACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTAcc                                                                                                                                                       >  2:4291834/1‑150 (MQ=255)
     aTTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAgg                                                                                                                                                  >  2:3031734/1‑150 (MQ=255)
     aTTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAgg                                                                                                                                                  >  2:4064005/1‑150 (MQ=255)
          cTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCAATGCGTACCGCAGggttgg                                                                                                                                             >  1:62938/1‑150 (MQ=255)
           tGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGggttggt                                                                                                                                            >  1:1577463/1‑150 (MQ=255)
           tGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGggttggt                                                                                                                                            >  1:794767/1‑150 (MQ=255)
             ggTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTtc                                                                                                                                          >  1:1531133/1‑150 (MQ=255)
             ggTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTtc                                                                                                                                          >  1:3129593/1‑150 (MQ=255)
              gTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTtct                                                                                                                                         <  2:20481/150‑1 (MQ=255)
               ttCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTtctc                                                                                                                                        <  2:3364263/150‑1 (MQ=255)
                    gTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTTCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTcagcag                                                                                                                                   >  2:4508989/1‑150 (MQ=255)
                     ttattaGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGa                                                                                                                                  <  2:1297834/150‑1 (MQ=255)
                      tattaGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGAt                                                                                                                                 >  1:4046417/1‑150 (MQ=255)
                      tattaGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGAt                                                                                                                                 >  2:2858291/1‑150 (MQ=255)
                        ttaGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCg                                                                                                                               <  2:1516098/150‑1 (MQ=255)
                         taGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGc                                                                                                                              <  2:4568916/150‑1 (MQ=255)
                         taGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGc                                                                                                                              <  2:3321545/150‑1 (MQ=255)
                              gCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGCTGata                                                                                                                         <  2:1171835/150‑1 (MQ=255)
                               cACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGatat                                                                                                                        <  2:3724218/150‑1 (MQ=255)
                                aCTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATAtt                                                                                                                       >  2:670437/1‑150 (MQ=255)
                                aCTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATAtt                                                                                                                       >  1:2521886/1‑150 (MQ=255)
                                 cTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGCTGATATTg                                                                                                                      <  2:1984666/150‑1 (MQ=255)
                                   aaTACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGcc                                                                                                                    >  2:3002135/1‑150 (MQ=255)
                                   aaTACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGcc                                                                                                                    >  1:271503/1‑150 (MQ=255)
                                     tACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCAt                                                                                                                  <  1:549681/150‑1 (MQ=255)
                                             atcaaATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATGTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCcggtc                                                                                                          >  2:2721107/1‑150 (MQ=255)
                                              tcacagcggCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCcggtcg                                                                                                         <  2:44219/143‑1 (MQ=255)
                                              tcaaATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCcggtcg                                                                                                         <  1:3629151/150‑1 (MQ=255)
                                               caaATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCcggtcgg                                                                                                        >  2:1382294/1‑150 (MQ=255)
                                                aaaTAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCcggtcggt                                                                                                       >  2:3884929/1‑150 (MQ=255)
                                                aaaTAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCcggtcggt                                                                                                       >  1:2109100/1‑150 (MQ=255)
                                                  aTAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTaa                                                                                                     >  1:4323271/1‑150 (MQ=255)
                                                    aGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAAc                                                                                                   <  2:2154845/150‑1 (MQ=255)
                                                      gCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGc                                                                                                 <  1:881348/150‑1 (MQ=255)
                                                             ggTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGc                                                                                          >  2:3284730/1‑150 (MQ=255)
                                                             ggTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGc                                                                                          >  2:7874/1‑150 (MQ=255)
                                                              gTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCt                                                                                         <  1:5131367/150‑1 (MQ=255)
                                                                aGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCttt                                                                                       <  1:5004551/150‑1 (MQ=255)
                                                                    cGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACt                                                                                   >  1:3871470/1‑150 (MQ=255)
                                                                    cGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACt                                                                                   <  1:2277240/150‑1 (MQ=255)
                                                                    cGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACt                                                                                   <  2:1306280/150‑1 (MQ=255)
                                                                    cGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACt                                                                                   <  1:3280099/150‑1 (MQ=255)
                                                                         cACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGgtgt                                                                              >  1:2935134/1‑150 (MQ=255)
                                                                         cACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGgtgt                                                                              >  2:1017118/1‑150 (MQ=255)
                                                                         cACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGgtgt                                                                              <  1:1288449/150‑1 (MQ=255)
                                                                          aCCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTa                                                                             >  1:3016512/1‑150 (MQ=255)
                                                                           ccGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTaa                                                                            >  2:2561816/1‑150 (MQ=255)
                                                                             gTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATg                                                                          >  2:1481284/1‑150 (MQ=255)
                                                                             gTGAATAGAGATTGTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATg                                                                          <  2:3124085/150‑1 (MQ=255)
                                                                                aaTAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCgt                                                                       >  2:2815174/1‑150 (MQ=255)
                                                                                 aTAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCgtg                                                                      >  1:5285039/1‑150 (MQ=255)
                                                                                   agagATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGtt                                                                    >  1:1851433/1‑150 (MQ=255)
                                                                                     agaTTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGc                                                                  >  2:1556593/1‑150 (MQ=255)
                                                                                       aTTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCg                                                                >  1:3343632/1‑150 (MQ=255)
                                                                                       aTTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCg                                                                >  1:1985216/1‑150 (MQ=255)
                                                                                          ttCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGttt                                                             <  2:4724484/150‑1 (MQ=255)
                                                                                           tctcAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTa                                                            >  1:4261894/1‑150 (MQ=255)
                                                                                           tctcAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTa                                                            >  1:1375048/1‑150 (MQ=255)
                                                                                           tctcAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTa                                                            <  1:2335798/150‑1 (MQ=255)
                                                                                            ctcAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTAc                                                           <  1:4905641/150‑1 (MQ=255)
                                                                                            ctcAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTAc                                                           >  2:4887910/1‑150 (MQ=255)
                                                                                              cAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACgg                                                         >  1:1650420/1‑150 (MQ=255)
                                                                                                   cGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACa                                                    >  1:3450378/1‑150 (MQ=255)
                                                                                                   cGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACa                                                    <  1:2842782/150‑1 (MQ=255)
                                                                                                    gAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAg                                                   >  2:2410502/1‑150 (MQ=255)
                                                                                                    gAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAg                                                   >  2:3760004/1‑150 (MQ=255)
                                                                                                      atatCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAggg                                                 <  2:4970547/150‑1 (MQ=255)
                                                                                                      atatCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAggg                                                 >  2:4995893/1‑150 (MQ=255)
                                                                                                      atatCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAggg                                                 <  1:1931018/150‑1 (MQ=255)
                                                                                                        atCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGAc                                               >  1:3749935/1‑150 (MQ=255)
                                                                                                        atCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGAc                                               >  2:1631753/1‑150 (MQ=255)
                                                                                                         tctcCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACa                                              >  1:1461882/1‑150 (MQ=255)
                                                                                                          ctcCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAg                                             >  1:2605834/1‑150 (MQ=255)
                                                                                                             ccTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAg                                          >  1:1053052/1‑150 (MQ=255)
                                                                                                               tGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCa                                        >  2:2958313/1‑150 (MQ=255)
                                                                                                                  ttCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAg                                     >  2:1791346/1‑150 (MQ=255)
                                                                                                                  ttCAACGAGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAg                                     >  1:2374269/1‑150 (MQ=255)
                                                                                                                   tCAACGAGGTACTGTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCAGTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGa                                    <  2:3432852/150‑1 (MQ=255)
                                                                                                                     aaCGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGAcc                                  >  2:825697/1‑150 (MQ=255)
                                                                                                                       cGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGt                                >  2:2933007/1‑150 (MQ=255)
                                                                                                                       cGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGt                                >  1:4708616/1‑150 (MQ=255)
                                                                                                                        gAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTc                               >  2:213552/1‑150 (MQ=255)
                                                                                                                        gAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTc                               <  1:3536480/150‑1 (MQ=255)
                                                                                                                        gAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTc                               >  2:4764360/1‑150 (MQ=255)
                                                                                                                         aGGTACTTTGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCa                              >  2:5088508/1‑150 (MQ=255)
                                                                                                                            tACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCAtct                           >  2:22454/1‑150 (MQ=255)
                                                                                                                             aCTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCAtctc                          <  1:3076036/150‑1 (MQ=255)
                                                                                                                               tttGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCtt                        <  1:4638412/150‑1 (MQ=255)
                                                                                                                               tttGCGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCtt                        >  2:587284/1‑150 (MQ=255)
                                                                                                                                   cGTTGAATCCATGCGTACAGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCaaa                    >  2:2708071/1‑150 (MQ=255)
                                                                                                                                       gAATCCATGCGTACCGCAGGGTTGGGTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCGCTTCaaaaaac                >  1:3564954/1‑149 (MQ=255)
                                                                                                                                         aTCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATg              <  2:2323485/150‑1 (MQ=255)
                                                                                                                                          tCCATGCGTACAGCAGGGGTGGTTCTCAGCAGATCGCTGATATAGCCATTCCCGGTCGGAAAACGCTCGAAGCTTTCACTGGTGTAATGCGTGTTGCCGTATACGCGAACAGGGACAGGAGCAAAGACCGTCATCTCTTCAAAAAAAtcc             >  2:5250736/1‑148 (MQ=255)
                                                                                                                                           ccATGAGTACAGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACTCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACATGAGCATATACCGTCATCTCTTCAAAAAAATGCa            <  1:1204522/150‑1 (MQ=255)
                                                                                                                                                cGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTAc       >  2:4866799/1‑150 (MQ=255)
                                                                                                                                                cGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTAc       >  2:3839083/1‑150 (MQ=255)
                                                                                                                                                 gTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACa      >  1:4354079/1‑150 (MQ=255)
                                                                                                                                                 gTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACa      >  1:2212658/1‑150 (MQ=255)
                                                                                                                                                  tACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACAt     <  1:207483/150‑1 (MQ=255)
                                                                                                                                                     cGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGtt  >  1:3669713/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATATCTCCCTGGTTCAACGAGGTACTTTGCGTTGAATCCATGCGTACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTT  >  CP009273/1571509‑1571807

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: