Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 86,729 G→T 15.3% R50L (CGT→CTT)  rsmH → 16S rRNA m(4)C1402 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927386,7290GT15.3% 887.3 / 92.6 275R50L (CGT→CTT) rsmH16S rRNA m(4)C1402 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (127/105);  new base T (17/26);  total (144/131)
Fisher's exact test for biased strand distribution p-value = 7.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 281 total aligned reads displayed.

AATGATGGAAAACTATAAACATACTACGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCC  >  CP009273/86580‑86878
                                                                                                                                                     |                                                                                                                                                     
  tgatgGAAAACTATAAACACACTACGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGtt                                                                                                                                                     <  2:3890820/150‑1 (MQ=255)
    atgGAAAACTATAAACATACTACGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCTTTTg                                                                                                                                                   <  1:3435872/150‑1 (MQ=255)
     tgGAAAACTATAAACATACTACGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCTTTTGc                                                                                                                                                  >  2:1693076/1‑150 (MQ=255)
      ggAAAACTATAAACATACTACGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCTTTTGCt                                                                                                                                                 >  1:360086/1‑150 (MQ=255)
       gAAAACTATAAACATACTACGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTg                                                                                                                                                >  1:90146/1‑150 (MQ=255)
       gAAAACTATAAACATACTACGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTg                                                                                                                                                >  2:3425678/1‑150 (MQ=255)
        aaaaCTATAAACATACTACGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTgg                                                                                                                                               <  2:222654/150‑1 (MQ=255)
        aaaaCTATAAACATACTACGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTgg                                                                                                                                               <  2:150879/150‑1 (MQ=255)
        aaaaCTATAAACATACTACGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTgg                                                                                                                                               <  1:1514143/150‑1 (MQ=255)
          aaCTATAAACATACTACGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGcg                                                                                                                                             <  2:764714/150‑1 (MQ=255)
           aCTATAAACATACTACGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGcga                                                                                                                                            >  2:927082/1‑150 (MQ=255)
             tataAACATACTACGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGcgatc                                                                                                                                          >  1:3117767/1‑150 (MQ=255)
             tataAACATACTACGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGcgatc                                                                                                                                          <  1:4214625/150‑1 (MQ=255)
             tataAACATACTACGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGACCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGcgatc                                                                                                                                          >  1:1843659/1‑150 (MQ=255)
               taAACATACTACGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCTTTTGCTGGcgatcga                                                                                                                                        <  1:3689939/150‑1 (MQ=255)
                aaaCATACTACGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGcgatcgat                                                                                                                                       >  2:4121820/1‑150 (MQ=255)
                aaaCATACTACGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGcgatcgat                                                                                                                                       >  1:4124384/1‑150 (MQ=255)
                aaaCATACTACGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGcgatcgat                                                                                                                                       <  1:691693/150‑1 (MQ=255)
                    aTACTACGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATcgcg                                                                                                                                   >  2:3270503/1‑150 (MQ=255)
                     tactacGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGa                                                                                                                                  >  2:1593451/1‑150 (MQ=255)
                     tactacGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGa                                                                                                                                  >  2:3275699/1‑150 (MQ=255)
                     tactacGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGa                                                                                                                                  <  2:463684/150‑1 (MQ=255)
                      actacGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGAc                                                                                                                                 <  1:3803327/150‑1 (MQ=255)
                       ctacGGTGCTGCTGGATGAAGCCGTTAATGGCCTCACTATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGCAGAGGGGCGTTTGCTGGCGATCGATCGCGAcc                                                                                                                                >  2:1975250/1‑150 (MQ=255)
                       ctacGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATAGATCGCGAcc                                                                                                                                >  2:4582770/1‑150 (MQ=255)
                             tgctgcTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAgg                                                                                                                          <  2:2792922/150‑1 (MQ=255)
                             tgctgcTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAgg                                                                                                                          <  2:2792834/150‑1 (MQ=255)
                               ctgctgGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCTTTTGCTGGCGATCGATCGCGACCCGCAGGCt                                                                                                                        >  2:3090469/1‑150 (MQ=255)
                               ctgctgGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCt                                                                                                                        <  2:2195872/150‑1 (MQ=255)
                                tgctgGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCTTTTGCTGGCGATCGATCGCGACCCGCAGGCTa                                                                                                                       >  1:1803542/1‑150 (MQ=255)
                                tgctgGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTa                                                                                                                       <  1:2669307/150‑1 (MQ=255)
                                tgctgGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTa                                                                                                                       >  2:508068/1‑150 (MQ=255)
                                  ctgGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATc                                                                                                                     <  1:2916858/150‑1 (MQ=255)
                                  ctgGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATc                                                                                                                     >  1:2706993/1‑150 (MQ=255)
                                   tgGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCg                                                                                                                    <  2:1999175/150‑1 (MQ=255)
                                   tgGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCg                                                                                                                    >  1:2557479/1‑150 (MQ=255)
                                    ggATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGc                                                                                                                   >  2:969316/1‑150 (MQ=255)
                                       tGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGt                                                                                                                >  1:1028468/1‑150 (MQ=255)
                                        gAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGtt                                                                                                               >  1:2824379/1‑150 (MQ=255)
                                         aaGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTg                                                                                                              >  1:3027470/1‑150 (MQ=255)
                                                  aTGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCGTGGCGTAGAGGGGCTTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTa                                                                                                     <  1:109098/150‑1 (MQ=255)
                                                   tGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTAt                                                                                                    >  2:3970480/1‑150 (MQ=255)
                                                      ccTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATtga                                                                                                 >  2:3141399/1‑150 (MQ=255)
                                                       cTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCTTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATtgat                                                                                                <  1:2967319/150‑1 (MQ=255)
                                                       cTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCTTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATtgat                                                                                                <  1:2967327/150‑1 (MQ=255)
                                                       cTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCTTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATtgat                                                                                                <  2:451898/150‑1 (MQ=255)
                                                           atatCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATc                                                                                            <  2:1012152/150‑1 (MQ=255)
                                                             atCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCcg                                                                                          <  1:581363/150‑1 (MQ=255)
                                                                  tCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCtt                                                                                     <  1:1556877/150‑1 (MQ=255)
                                                                  tCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCtt                                                                                     <  2:2455545/150‑1 (MQ=255)
                                                                   ccTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCTTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTtc                                                                                    >  2:4065014/1‑150 (MQ=255)
                                                                    cTGATGGCATCTACAGTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTtct                                                                                   >  1:2306683/1‑150 (MQ=255)
                                                                      gATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTcc                                                                                 <  1:170066/150‑1 (MQ=255)
                                                                        tGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCTTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCcat                                                                               >  2:2681534/1‑150 (MQ=255)
                                                                        tGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCcat                                                                               >  1:1935706/1‑150 (MQ=255)
                                                                          gCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCTTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCcatca                                                                             <  1:328936/150‑1 (MQ=255)
                                                                          gCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCcatca                                                                             <  1:4007883/150‑1 (MQ=255)
                                                                          gCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCcatca                                                                             <  1:1198185/150‑1 (MQ=255)
                                                                           cATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCcatcat                                                                            <  2:131338/150‑1 (MQ=255)
                                                                                  aTTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCTTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGa                                                                     <  1:2663255/150‑1 (MQ=255)
                                                                                  aTTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGa                                                                     >  2:1654981/1‑150 (MQ=255)
                                                                                  aTTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGa                                                                     >  1:769889/1‑150 (MQ=255)
                                                                                    tGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGAcc                                                                   >  2:1764886/1‑150 (MQ=255)
                                                                                     gATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCt                                                                  <  1:3939926/150‑1 (MQ=255)
                                                                                     gATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCt                                                                  <  1:3541082/150‑1 (MQ=255)
                                                                                       tGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCttt                                                                <  1:4136862/150‑1 (MQ=255)
                                                                                       tGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCttt                                                                <  2:3783806/150‑1 (MQ=255)
                                                                                       tGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCttt                                                                <  2:1463359/150‑1 (MQ=255)
                                                                                        gggACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTtc                                                               >  1:1407789/1‑150 (MQ=255)
                                                                                                tGGTCGCGGTTGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGTGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGTCTATTGATGATCCGCGCTTCTCCATAATCCACGGACCTTTCTCCGCGCt                                                       <  1:4328986/150‑1 (MQ=255)
                                                                                                tGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCt                                                       >  1:4522363/1‑150 (MQ=255)
                                                                                                  gTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTgg                                                     >  2:2961435/1‑150 (MQ=255)
                                                                                                    cgcgGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCTTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGc                                                   >  2:3481462/1‑150 (MQ=255)
                                                                                                    cgcgGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGc                                                   >  2:4620198/1‑150 (MQ=255)
                                                                                                     gcgGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCTTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCg                                                  >  2:3340168/1‑150 (MQ=255)
                                                                                                      cgGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGa                                                 >  1:1763250/1‑150 (MQ=255)
                                                                                                            tcactcacGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCTCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGt                                           >  2:4238198/1‑150 (MQ=255)
                                                                                                             cactcacGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGtt                                          >  2:2906510/1‑150 (MQ=255)
                                                                                                                tcacGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGcc                                       <  1:3172918/150‑1 (MQ=255)
                                                                                                                tcacGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGcc                                       <  2:4459639/150‑1 (MQ=255)
                                                                                                                  acGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCTTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACGATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGa                                     >  1:1610615/1‑150 (MQ=255)
                                                                                                                  acGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGa                                     >  1:2282327/1‑150 (MQ=255)
                                                                                                                    gTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAgc                                   >  2:3738745/1‑150 (MQ=255)
                                                                                                                       tGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGcgcg                                >  2:3392775/1‑150 (MQ=255)
                                                                                                                        gATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGa                               >  1:4261699/1‑150 (MQ=255)
                                                                                                                          tCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATc                             >  1:921997/1‑150 (MQ=255)
                                                                                                                           ccTCTCGCAGCTTGGCGAAGAGGGGCTTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCGTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGTGCGCGATCt                            <  2:1190785/150‑1 (MQ=255)
                                                                                                                              ctcGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTAt                         <  2:3028685/150‑1 (MQ=255)
                                                                                                                                cGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCg                       >  1:3699072/1‑150 (MQ=255)
                                                                                                                                      ttGGCGAAGAGGGGCTTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGa                 <  1:2076707/150‑1 (MQ=255)
                                                                                                                                      ttGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGa                 <  2:318684/150‑1 (MQ=255)
                                                                                                                                       tGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCGTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGAt                <  1:867769/150‑1 (MQ=255)
                                                                                                                                        ggCGAAGAGGGGCTTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATc               <  1:3279143/150‑1 (MQ=255)
                                                                                                                                          cGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGa             <  2:2115714/150‑1 (MQ=255)
                                                                                                                                           gAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGAc            <  1:503786/150‑1 (MQ=255)
                                                                                                                                            aaGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACg           <  2:192292/150‑1 (MQ=255)
                                                                                                                                            aaGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACg           >  2:1799746/1‑150 (MQ=255)
                                                                                                                                             agagGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACgg          >  2:3662105/1‑150 (MQ=255)
                                                                                                                                              gagGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCGTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGc         <  2:261021/150‑1 (MQ=255)
                                                                                                                                                     gTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTctcc  >  2:848889/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATGATGGAAAACTATAAACATACTACGGTGCTGCTGGATGAAGCCGTTAATGGCCTCAATATCCGTCCTGATGGCATCTACATTGATGGGACTTTTGGTCGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCC  >  CP009273/86580‑86878

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: