Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,376,178 C→A 46.2% E1004* (GAA→TAA)  mscM ← mechanosensitive channel protein, miniconductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,376,1780CA46.2% 180.5 / inf 292E1004* (GAA→TAA) mscMmechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base C (73/84);  new base A (64/71);  total (137/155)
Fisher's exact test for biased strand distribution p-value = 9.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 295 total aligned reads displayed.

ATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGATGGTGGTGGCGCGGG  >  CP009273/4376030‑4376327
                                                                                                                                                    |                                                                                                                                                     
aTACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCg                                                                                                                                                      >  1:3006565/1‑150 (MQ=255)
      tGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTg                                                                                                                                                >  1:676783/1‑150 (MQ=255)
       gACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTgg                                                                                                                                               <  1:3286008/150‑1 (MQ=255)
         cccATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCa                                                                                                                                             >  1:1946523/1‑150 (MQ=255)
            aTCTCAGCGGCGTAAATACTCAGCTCGAAAATCTGAATACCCTGTTGCAGATCCACCATGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGTATTACCGTCACTTCTTAGCTATTGGCATCg                                                                                                                                          <  2:943665/150‑1 (MQ=255)
            aTCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGACTCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCg                                                                                                                                          <  1:2643054/150‑1 (MQ=255)
             tctcAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCgg                                                                                                                                         <  1:712277/150‑1 (MQ=255)
              ctcAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGc                                                                                                                                        >  1:576599/1‑150 (MQ=255)
                 aGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAgg                                                                                                                                     >  2:906179/1‑150 (MQ=255)
                     gcgTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGcc                                                                                                                                 <  2:3195741/150‑1 (MQ=255)
                       gTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCgg                                                                                                                               >  1:3194689/1‑150 (MQ=255)
                        tAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGt                                                                                                                              >  1:3833387/1‑150 (MQ=255)
                         aaaTACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTa                                                                                                                             >  1:1093756/1‑150 (MQ=255)
                          aaTACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTAt                                                                                                                            >  2:3128408/1‑150 (MQ=255)
                          aaTACGCAGCTCGAAAATCTGAATCCCCTGGTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTAt                                                                                                                            <  1:3304667/150‑1 (MQ=255)
                              cGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTc                                                                                                                        <  2:14866/150‑1 (MQ=255)
                               gCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCa                                                                                                                       <  1:4266490/150‑1 (MQ=255)
                                cAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCaa                                                                                                                      >  1:836927/1‑150 (MQ=255)
                                 aGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAac                                                                                                                     >  1:342515/1‑150 (MQ=255)
                                   cTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAacac                                                                                                                   >  2:4135711/1‑150 (MQ=255)
                                   cTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAacac                                                                                                                   >  2:2197460/1‑150 (MQ=255)
                                    tCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAacagc                                                                                                                  <  1:4320756/150‑3 (MQ=255)
                                        aaaTCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCAcgc                                                                                                              >  1:3695905/1‑150 (MQ=255)
                                          aTCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATCTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCAcgcgc                                                                                                            >  2:2531011/1‑150 (MQ=255)
                                            cTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGt                                                                                                          <  2:2626261/150‑1 (MQ=255)
                                            cTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGt                                                                                                          >  1:3021253/1‑150 (MQ=255)
                                             tGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTg                                                                                                         >  2:1904313/1‑150 (MQ=255)
                                                agccccTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGAcc                                                                                                      <  1:2974256/148‑1 (MQ=255)
                                                aTCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGAcc                                                                                                      <  1:1906395/150‑1 (MQ=255)
                                                aTCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGAcc                                                                                                      <  1:2979162/150‑1 (MQ=255)
                                                aTCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGAcc                                                                                                      >  1:2854854/1‑150 (MQ=255)
                                                 tCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCg                                                                                                     >  1:1866840/1‑150 (MQ=255)
                                                   cccTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGATGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAg                                                                                                   >  2:160347/1‑150 (MQ=255)
                                                      tGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCa                                                                                                <  2:2982924/150‑1 (MQ=255)
                                                      tGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCa                                                                                                >  1:2545481/1‑150 (MQ=255)
                                                      tGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCa                                                                                                >  2:2304524/1‑150 (MQ=255)
                                                       gTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCag                                                                                               >  1:4100556/1‑150 (MQ=255)
                                                       gTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCag                                                                                               <  1:2558905/150‑1 (MQ=255)
                                                        ttGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCaga                                                                                              <  1:2795989/150‑1 (MQ=255)
                                                           cAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCagagag                                                                                           <  2:408866/150‑1 (MQ=255)
                                                            aGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGc                                                                                          <  2:565818/150‑1 (MQ=255)
                                                            aGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGc                                                                                          >  1:1201111/1‑150 (MQ=255)
                                                            aGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGc                                                                                          >  2:1387206/1‑150 (MQ=255)
                                                              aTCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGa                                                                                        <  2:3444598/150‑1 (MQ=255)
                                                               tCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGAc                                                                                       <  2:2805556/150‑1 (MQ=255)
                                                               tCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGAc                                                                                       >  2:184561/1‑150 (MQ=255)
                                                               tCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGAc                                                                                       >  2:1256039/1‑150 (MQ=255)
                                                               tCCACCAGGAAGACTTCCGGTGCAGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGAc                                                                                       <  2:1287136/150‑1 (MQ=255)
                                                                 caccaGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCa                                                                                     >  2:437984/1‑150 (MQ=255)
                                                                 caccaGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCa                                                                                     <  2:3866031/150‑1 (MQ=255)
                                                                    caGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGtt                                                                                  >  2:3994299/1‑150 (MQ=255)
                                                                    caGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGtt                                                                                  <  2:2347163/150‑1 (MQ=255)
                                                                         aGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATaa                                                                             <  2:2676370/150‑1 (MQ=255)
                                                                          gACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATaaa                                                                            <  2:3084738/150‑1 (MQ=255)
                                                                          gACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATaaa                                                                            >  2:4130734/1‑150 (MQ=255)
                                                                            cTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACt                                                                          <  1:1514746/150‑1 (MQ=255)
                                                                               ccGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCCGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCt                                                                       >  2:2735176/1‑150 (MQ=255)
                                                                               ccGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCt                                                                       <  2:2097123/150‑1 (MQ=255)
                                                                                cGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTCTCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTc                                                                      >  2:2291006/1‑150 (MQ=255)
                                                                                cGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTc                                                                      >  2:116169/1‑150 (MQ=255)
                                                                                  gTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGGGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCgg                                                                    <  1:1019922/150‑1 (MQ=255)
                                                                                   tGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGt                                                                   >  1:778490/1‑150 (MQ=255)
                                                                                    gCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGta                                                                  <  1:1621741/150‑1 (MQ=255)
                                                                                      cGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGtaat                                                                >  1:3515623/1‑150 (MQ=255)
                                                                                       gggTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGtaata                                                               >  2:1539908/1‑150 (MQ=255)
                                                                                        ggTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGtaataa                                                              <  2:1011259/150‑1 (MQ=255)
                                                                                        ggTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGtaataa                                                              <  2:3638322/150‑1 (MQ=255)
                                                                                           tGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACg                                                           <  2:428442/150‑1 (MQ=255)
                                                                                           tGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTATTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACg                                                           <  1:4185200/150‑1 (MQ=255)
                                                                                            gTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGc                                                          <  2:3239481/150‑1 (MQ=255)
                                                                                             tCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGcc                                                         <  2:3540553/150‑1 (MQ=255)
                                                                                              cGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCt                                                        <  1:3716306/150‑1 (MQ=255)
                                                                                               gATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCtt                                                       <  1:1706761/150‑1 (MQ=255)
                                                                                               gATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCtt                                                       >  1:586523/1‑150 (MQ=255)
                                                                                               gATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCtt                                                       <  2:824492/150‑1 (MQ=255)
                                                                                               gATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCtt                                                       <  2:1787085/150‑1 (MQ=255)
                                                                                                  cACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGtt                                                    >  1:2244024/1‑150 (MQ=255)
                                                                                                   aCCAGCGAGCCGCGACGCGCTGCGGTGAGCTGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGGGACCGAGCCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTc                                                   <  1:1602114/150‑1 (MQ=255)
                                                                                                   aCCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTc                                                   <  1:2661389/150‑1 (MQ=255)
                                                                                                      aGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGc                                                >  1:2360984/1‑150 (MQ=255)
                                                                                                         gAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACg                                             <  1:2770584/150‑1 (MQ=255)
                                                                                                            cagcGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGAtt                                          >  2:2225594/1‑150 (MQ=255)
                                                                                                            cagcGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGAtt                                          <  1:905903/150‑1 (MQ=255)
                                                                                                             agcGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTa                                         >  1:970608/1‑150 (MQ=255)
                                                                                                             agcGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTa                                         >  2:3777113/1‑150 (MQ=255)
                                                                                                                gACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTAtct                                      >  1:2459031/1‑150 (MQ=255)
                                                                                                                    cgcTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCttt                                  <  1:1382843/150‑1 (MQ=255)
                                                                                                                      cTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTAc                                >  2:257793/1‑150 (MQ=255)
                                                                                                                        gCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACgg                              <  1:496056/150‑1 (MQ=255)
                                                                                                                         cGGTGAGCAGGATGTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGt                             <  2:2178596/150‑1 (MQ=255)
                                                                                                                          ggTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTc                            >  2:290484/1‑150 (MQ=255)
                                                                                                                          ggTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTc                            >  2:49015/1‑150 (MQ=255)
                                                                                                                               gCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGGGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGt                       <  2:1991782/150‑1 (MQ=255)
                                                                                                                                       tCCGTCAATTCTTCGCTTTTGGCATCGGCAGGGGCCGTTATCGTCAACACCACGCGCGTGACCGCGTAAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtg               <  1:1818084/150‑1 (MQ=255)
                                                                                                                                        ccGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtgg              <  1:2071750/150‑1 (MQ=255)
                                                                                                                                        ccGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtgg              <  2:2413394/150‑1 (MQ=255)
                                                                                                                                        ccGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtgg              >  2:691959/1‑150 (MQ=255)
                                                                                                                                         cGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggt             <  1:2304524/150‑1 (MQ=255)
                                                                                                                                             aCTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggtggtg         <  1:2636551/150‑1 (MQ=255)
                                                                                                                                                  ttcGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGATGGTGGTGGcgcg    <  2:3393954/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGATGGTGGTGGCGCGGG  >  CP009273/4376030‑4376327

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: