Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,503,813 A→C 5.0% A235A (GCA→GCC yfeO → putative ion channel protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,503,8130AC 5.0% 895.5 / 34.7 279A235A (GCA→GCCyfeOputative ion channel protein
Reads supporting (aligned to +/- strand):  ref base A (140/125);  new base C (2/12);  total (142/137)
Fisher's exact test for biased strand distribution p-value = 5.36e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 293 total aligned reads displayed.

CTCTTTGCGCCGTTAATGGCGGCAGCAGCTGGTGCACTTACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGG  >  CP009273/2503664‑2503961
                                                                                                                                                     |                                                                                                                                                    
ctctTTGCGCCGTTAATGGCGGCAGCAGCTGGTGCACTTACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCgca                                                                                                                                                      >  2:3196298/1‑150 (MQ=255)
 tctTTGCGCCGTTAATGGCGGCAGCAGCTGGTGCACTTACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCgcag                                                                                                                                                     >  2:543818/1‑150 (MQ=255)
  ctTTGCGCCGTTAATGGCGGCAGCAGCTGGTGCACTTACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCgcagc                                                                                                                                                    >  1:2689245/1‑150 (MQ=255)
     tGCGCCGTTAATGGCGGCAGCAGCTGGTGCACTTACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAgg                                                                                                                                                 >  2:3300047/1‑150 (MQ=255)
        gcCGTTAATGGCGGCAGCAGCTGGTGCACTTACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGAt                                                                                                                                              <  2:3313809/150‑1 (MQ=255)
         ccGTTAATGGCGGCAGCAGCTGGTGCACTTACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGTGATg                                                                                                                                             <  2:3630418/150‑1 (MQ=255)
           gTTAATGGCGGCAGCAGCTGGTGCACTTACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGt                                                                                                                                           >  2:1283574/1‑150 (MQ=255)
           gTTAATGGCGGCAGCAGCTGGTGCACTTACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGt                                                                                                                                           <  1:1561116/150‑1 (MQ=255)
            ttAATGGCGGCAGCAGCTGGTGCACTTACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTc                                                                                                                                          >  2:359456/1‑150 (MQ=255)
             tAATGGCGGCAGCAGCTGGTGCACTTACTACCGGATTATTTGTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCg                                                                                                                                         <  1:3870177/150‑1 (MQ=255)
               aTGGCGGCAGCAGCTGGTGCACTTACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGcc                                                                                                                                       <  1:2249739/150‑1 (MQ=255)
                tGGCGGCAGCAGCTGGTGCACTTACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCg                                                                                                                                      <  2:1031653/150‑1 (MQ=255)
                    ggcAGCAGCTGGTGCACTTACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATg                                                                                                                                  >  1:836832/1‑150 (MQ=255)
                    ggcAGCAGCTGGTGCACTTACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCACATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCATGGATGGTCTCCGTATg                                                                                                                                  >  2:172591/1‑150 (MQ=255)
                         cagcTGGTGCACTTACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCt                                                                                                                             >  2:1664137/1‑150 (MQ=255)
                         cagcTGGTGCACTTACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCt                                                                                                                             <  2:3270492/150‑1 (MQ=255)
                         cagcTGGTGCACTTACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCt                                                                                                                             >  1:2210127/1‑150 (MQ=255)
                             tGGTGCACTTACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTAcc                                                                                                                         >  2:2139056/1‑150 (MQ=255)
                             tGGTGCACTTACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCCCAATGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATAGTCGCGGCGATTGCCAACGCAGCAGGGATGGTCGCCGTATGGTGCTTAcc                                                                                                                         >  2:482191/1‑150 (MQ=255)
                                tGCACTTACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACg                                                                                                                      >  1:1662728/1‑150 (MQ=255)
                                   aCTTACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGtt                                                                                                                   <  1:3043917/150‑1 (MQ=255)
                                   aCTTACTACCGGATTATATTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGtt                                                                                                                   <  2:697659/150‑1 (MQ=255)
                                    cTTACGACCGGATTATTTTTCCATCCTCCTTGTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTg                                                                                                                  <  2:2224584/150‑1 (MQ=255)
                                     ttACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGc                                                                                                                 <  1:3332044/150‑1 (MQ=255)
                                          acCGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCAcgcg                                                                                                            >  2:3875903/1‑150 (MQ=255)
                                           ccGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCga                                                                                                           >  1:457004/1‑150 (MQ=255)
                                             ggATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCgatg                                                                                                         <  2:1464113/150‑1 (MQ=255)
                                             ggATTATTTTTCCATCCTCATGTTTCACTGCCCAGTGCTCATTACGGACAGATGGAAATGACCGATAGTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCgatg                                                                                                         <  2:3388032/150‑1 (MQ=255)
                                              gATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCgatga                                                                                                        >  1:1198848/1‑150 (MQ=255)
                                              gATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCgatga                                                                                                        >  1:3292084/1‑150 (MQ=255)
                                               attattTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGTCGATTGCCATCGCAGCAGGGATGGTCGCAGTATGGTGCTTACCACGGTTGCACGCgatgat                                                                                                       >  1:676494/1‑150 (MQ=255)
                                               attattTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCgatgat                                                                                                       <  1:3059832/150‑1 (MQ=255)
                                                ttattTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGatgatg                                                                                                      <  1:1015145/150‑1 (MQ=255)
                                                ttattTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGatgatg                                                                                                      <  2:1836173/150‑1 (MQ=255)
                                                ttattTTTCCATCCTCATTTTTCACTGCCCATTGATCATTACGGACAGATTGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGatgatg                                                                                                      <  1:3106756/150‑1 (MQ=255)
                                                  attTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCa                                                                                                    <  1:2849626/150‑1 (MQ=255)
                                                       tCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCaaa                                                                                               >  1:3735671/1‑150 (MQ=255)
                                                        ccATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAAt                                                                                              >  1:3271606/1‑150 (MQ=255)
                                                         cATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATg                                                                                             >  1:3148908/1‑150 (MQ=255)
                                                         cATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATg                                                                                             <  1:1140115/150‑1 (MQ=255)
                                                             cTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGaaaa                                                                                         >  2:1001130/1‑150 (MQ=255)
                                                             cTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGaaaa                                                                                         >  2:3665227/1‑150 (MQ=255)
                                                              tCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGaaaaa                                                                                        >  1:1809391/1‑150 (MQ=255)
                                                              tCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGaaaaa                                                                                        >  1:490308/1‑150 (MQ=255)
                                                                 tttttCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATcc                                                                                     <  2:2603433/150‑1 (MQ=255)
                                                                   tttCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCgg                                                                                   >  1:135914/1‑150 (MQ=255)
                                                                   tttCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCgg                                                                                   <  1:701251/150‑1 (MQ=255)
                                                                        cTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTc                                                                              >  2:2720920/1‑150 (MQ=255)
                                                                          gCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGt                                                                            >  1:2550196/1‑150 (MQ=255)
                                                                           cccATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTg                                                                           <  1:3638652/150‑1 (MQ=255)
                                                                            ccATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGc                                                                          >  2:1374858/1‑150 (MQ=255)
                                                                            ccATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAAGTGTCGCGGCGATTGCCATCGCAGCCGGGAGGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGc                                                                          >  2:3981097/1‑150 (MQ=255)
                                                                             cATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCt                                                                         <  1:4089341/150‑1 (MQ=255)
                                                                             cATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCt                                                                         >  2:2356927/1‑150 (MQ=255)
                                                                             cATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCt                                                                         <  2:921531/150‑1 (MQ=255)
                                                                              aTTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTg                                                                        >  1:3389654/1‑150 (MQ=255)
                                                                                tGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTggg                                                                      >  1:2120701/1‑150 (MQ=255)
                                                                                tGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATACGGTGCTCGTGCTggg                                                                      >  1:2666333/1‑150 (MQ=255)
                                                                                  cTCAGTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTa                                                                    <  1:151597/150‑1 (MQ=255)
                                                                                   tCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTAt                                                                   <  1:979121/150‑1 (MQ=255)
                                                                                   tCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTAt                                                                   <  1:2456250/150‑1 (MQ=255)
                                                                                   tCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTAt                                                                   >  2:1092350/1‑150 (MQ=255)
                                                                                    cATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCCGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTAtt                                                                  <  1:3314107/150‑1 (MQ=255)
                                                                                      ttACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTgg                                                                <  1:2903131/150‑1 (MQ=255)
                                                                                           gACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGAt                                                           <  1:559813/150‑1 (MQ=255)
                                                                                            aCAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCCTCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGAtt                                                          >  1:1061290/1‑150 (MQ=255)
                                                                                             cAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCCGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGAttt                                                         <  2:351309/150‑1 (MQ=255)
                                                                                             cAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCAGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGAttt                                                         >  1:2343713/1‑150 (MQ=255)
                                                                                               gATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCCGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTAt                                                       <  2:3002832/150‑1 (MQ=255)
                                                                                                  ggAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATtct                                                    <  2:2985526/150‑1 (MQ=255)
                                                                                                    aaaTGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCg                                                  >  2:1635337/1‑150 (MQ=255)
                                                                                                     aaTGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCgg                                                 <  1:1250382/150‑1 (MQ=255)
                                                                                                        gACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTAt                                              >  1:4015999/1‑150 (MQ=255)
                                                                                                        gACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTAt                                              <  2:4017819/150‑1 (MQ=255)
                                                                                                        gACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTAt                                              >  1:3842122/1‑150 (MQ=255)
                                                                                                         aCCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTAtt                                             <  1:1409707/150‑1 (MQ=255)
                                                                                                           cGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCCGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCt                                           <  2:810838/150‑1 (MQ=255)
                                                                                                            gATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGGTTTATTCTCGGGATTCTg                                          >  1:2140671/1‑150 (MQ=255)
                                                                                                            gATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGAGTTATTCTCGGTATTCTg                                          >  2:112745/1‑150 (MQ=255)
                                                                                                               atTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTgggg                                       >  2:3178284/1‑150 (MQ=255)
                                                                                                               atTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTgggg                                       <  2:145630/150‑1 (MQ=255)
                                                                                                                    cAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTAt                                  <  2:3585888/150‑1 (MQ=255)
                                                                                                                     aGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTAtt                                 >  2:2804672/1‑150 (MQ=255)
                                                                                                                       cGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATtgg                               >  2:3694298/1‑150 (MQ=255)
                                                                                                                       cGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATtgg                               <  1:3549906/150‑1 (MQ=255)
                                                                                                                         gTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATtggtg                             <  1:1278906/150‑1 (MQ=255)
                                                                                                                          tGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATtggtgg                            <  1:2512333/150‑1 (MQ=255)
                                                                                                                             aaTTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTCTTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGAcc                         >  1:2201745/1‑150 (MQ=255)
                                                                                                                                tGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGt                      >  2:2408716/1‑150 (MQ=255)
                                                                                                                                   cgcgGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGTTTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTc                   <  2:594742/150‑1 (MQ=255)
                                                                                                                                   cgcgGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTc                   <  1:3966503/150‑1 (MQ=255)
                                                                                                                                     cggcgATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGc                 >  2:19630/1‑150 (MQ=255)
                                                                                                                                      ggcgATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCt                >  2:3370616/1‑150 (MQ=255)
                                                                                                                                          aTTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGttt            >  2:2930705/1‑150 (MQ=255)
                                                                                                                                           ttGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTa           <  2:2046576/150‑1 (MQ=255)
                                                                                                                                               cATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAgg       <  2:2816579/150‑1 (MQ=255)
                                                                                                                                                 tCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGc     >  2:2842123/1‑150 (MQ=255)
                                                                                                                                                   gcagcaGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTg   >  2:3829075/1‑150 (MQ=255)
                                                                                                                                                   gcagcaGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTg   <  1:2531348/150‑1 (MQ=255)
                                                                                                                                                    cagcagGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTgg  <  1:1950957/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTCTTTGCGCCGTTAATGGCGGCAGCAGCTGGTGCACTTACTACCGGATTATTTTTCCATCCTCATTTTTCACTGCCCATTGCTCATTACGGACAGATGGAAATGACCGATATTCTCAGCGGTGCAATTGTCGCGGCGATTGCCATCGCAGCAGGGATGGTCGCCGTATGGTGCTTACCACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGG  >  CP009273/2503664‑2503961

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: