Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,981,203 T→C 5.1% intergenic (+99/‑201) rpsA → / → ihfB 30S ribosomal protein S1/integration host factor subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,981,2030TC 5.1% 494.3 / 9.8 157intergenic (+99/‑201)rpsA/ihfB30S ribosomal protein S1/integration host factor subunit beta
Reads supporting (aligned to +/- strand):  ref base T (91/56);  new base C (4/4);  total (97/60)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

Only 100 of 162 total aligned reads displayed.

ACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGTTCAAATCCTTCCCGGCGTG  >  NC_002947/1981057‑1981352
                                                                                                                                                  |                                                                                                                                                     
   accaTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGc                                                                                                                                                 <  2:1310435/150‑1 (MQ=255)
     caTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTg                                                                                                                                               <  1:639392/150‑1 (MQ=255)
     caTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTg                                                                                                                                               >  1:312444/1‑150 (MQ=255)
       tGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGc                                                                                                                                             <  1:2450681/150‑1 (MQ=255)
        ggCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAATGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCCCTGCTGGAACACGCCAGCTTCAGGCTGGCg                                                                                                                                            >  1:3660721/1‑150 (MQ=255)
           tgcgctgcTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAg                                                                                                                                         >  2:2938390/1‑150 (MQ=255)
           tgcgctgcTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAg                                                                                                                                         >  1:3590425/1‑150 (MQ=255)
             cgctgcTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCGCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCCTGAAGCTGGCgtgttc                                                                                                                                       <  2:3451648/150‑6 (MQ=255)
              gctgctGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGccc                                                                                                                                      >  2:98307/1‑150 (MQ=255)
              gctgctGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGccc                                                                                                                                      >  2:4121866/1‑150 (MQ=255)
                 gctgcGCTAAGCAATGGCCAAGCAGAACTGAGTTCTTTTTTATCGTTAAAATGGGCGACCCTCGGGTCGCCATTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGtt                                                                                                                                   <  1:1271747/150‑1 (MQ=255)
                    gcgcGACGCATTGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTa                                                                                                                                <  1:3559650/150‑1 (MQ=255)
                               tGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTc                                                                                                                     >  1:4973402/1‑150 (MQ=255)
                                 gCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAc                                                                                                                   >  1:4313766/1‑150 (MQ=255)
                                 gCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAc                                                                                                                   >  1:4308445/1‑150 (MQ=255)
                                 gCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAc                                                                                                                   >  1:4058359/1‑150 (MQ=255)
                                   cAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCggg                                                                                                                        >  1:2061229/1‑143 (MQ=255)
                                   cAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAccc                                                                                                                 >  2:625814/1‑150 (MQ=255)
                                     aGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCcgc                                                                                                               <  1:953978/150‑1 (MQ=255)
                                        cgAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGaa                                                                                                            <  2:3410169/149‑1 (MQ=255)
                                        aGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGaa                                                                                                            <  1:3896066/150‑1 (MQ=255)
                                            cTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAg                                                                                                        <  1:1384512/150‑1 (MQ=255)
                                               aGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGGTGTAGGAGCGGGTTCACCCGCGAAGCAGgcg                                                                                                     >  1:2788869/1‑150 (MQ=255)
                                                gTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGgcgc                                                                                                    <  2:4406249/150‑1 (MQ=255)
                                                   cTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCAc                                                                                                 >  2:2167868/1‑150 (MQ=255)
                                                    tGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACt                                                                                                <  2:3476279/150‑1 (MQ=255)
                                                      tttGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTga                                                                                              >  1:1482715/1‑149 (MQ=255)
                                                      tttGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCCTGAAGCTg                                                                                                                                               <  2:3574055/101‑1 (MQ=255)
                                                       ttGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGc                                                                                             <  2:142109/150‑1 (MQ=255)
                                                       ttGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGc                                                                                             <  1:3617568/150‑1 (MQ=255)
                                                         gATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCt                                                                                           >  2:1319597/1‑150 (MQ=255)
                                                         gATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCt                                                                                           <  1:578071/150‑1 (MQ=255)
                                                                 aaaaGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGc                                                                                   <  1:3760760/150‑1 (MQ=255)
                                                                 aaaaGGGCGACCCGAGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGATGCAGGCGCCACTGGCCTATCTAAGc                                                                                   >  2:2325709/1‑150 (MQ=255)
                                                                  aaaGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCg                                                                                  >  1:2674386/1‑150 (MQ=255)
                                                                  aaaGGGCGACCCGAGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCg                                                                                  >  1:3742410/1‑150 (MQ=255)
                                                                   aaGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCgg                                                                                 >  1:390967/1‑150 (MQ=255)
                                                                   aaGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCgg                                                                                 <  1:2317908/150‑1 (MQ=255)
                                                                    aGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGt                                                                                >  1:1021297/1‑150 (MQ=255)
                                                                    aGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGt                                                                                >  2:671179/1‑150 (MQ=255)
                                                                    aGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGt                                                                                >  1:1021429/1‑150 (MQ=255)
                                                                     gggCGACCCGAGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTa                                                                               >  1:3878747/1‑150 (MQ=255)
                                                                      ggCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAg                                                                              >  2:3276205/1‑150 (MQ=255)
                                                                       gCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGa                                                                             <  2:634578/150‑1 (MQ=255)
                                                                       gCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGa                                                                             >  2:3205887/1‑150 (MQ=255)
                                                                       gCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGa                                                                             <  2:1516168/150‑1 (MQ=255)
                                                                       gCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGa                                                                             <  1:982154/150‑1 (MQ=255)
                                                                        cGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGAc                                                                            >  2:4377622/1‑150 (MQ=255)
                                                                         gACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGAcc                                                                           >  2:215561/1‑150 (MQ=255)
                                                                         gACCCTCGGGTCGCCCTTTTTTTTGCTCCCCGACCAGGTGGTTCGGCCGATTGCAATGCATGAACACCCCAGCGTAAGGCTGGCGAAACCCGTTGTAGGAGCGGGTTAAACCCCCAAGCAAGCGCCACAGGCCAATAAAAGCGGTAGAca                                                                           >  1:2178740/1‑149 (MQ=255)
                                                                            ccTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTaa                                                                                     <  2:475002/137‑1 (MQ=255)
                                                                            ccTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTaa                                                                                     >  1:475002/1‑137 (MQ=255)
                                                                            ccTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCAtt                                                                        >  1:1794017/1‑150 (MQ=255)
                                                                            ccTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCAtt                                                                        >  1:1793891/1‑150 (MQ=255)
                                                                              tCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCt                                                                      >  2:207618/1‑150 (MQ=255)
                                                                              gagggTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTCAGCGGTAGACCATTCt                                                                      >  2:1121788/3‑150 (MQ=255)
                                                                              gagggTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCt                                                                      >  2:1985366/3‑150 (MQ=255)
                                                                              gagggTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCt                                                                      >  2:1985594/3‑150 (MQ=255)
                                                                              gagggTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCt                                                                      >  1:4571079/3‑150 (MQ=255)
                                                                               cGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCtt                                                                     <  1:1191914/150‑1 (MQ=255)
                                                                                 ggTCGCCCCTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCCaacg                                                                                                                                        >  1:3148711/1‑79 (MQ=255)
                                                                                 ggTCGCCCCTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCCaacg                                                                                                                                        <  2:3148711/81‑3 (MQ=255)
                                                                                  gTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGaa                                                                  <  2:4983669/150‑1 (MQ=255)
                                                                                  gTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGaa                                                                  <  2:278921/150‑1 (MQ=255)
                                                                                   tCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAAt                                                                 <  1:4413780/150‑1 (MQ=255)
                                                                                   tCGCCATTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAAt                                                                 <  1:4413800/150‑1 (MQ=255)
                                                                                    cGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAAtt                                                                <  1:3829225/150‑1 (MQ=255)
                                                                                     gCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTa                                                               >  1:3080510/1‑150 (MQ=255)
                                                                                      cccTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGGACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTAc                                                              >  2:4710581/1‑150 (MQ=255)
                                                                                        cTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACtt                                                            <  2:2257959/150‑1 (MQ=255)
                                                                                        cTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACtt                                                            <  1:2186743/150‑1 (MQ=255)
                                                                                        cTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACtg                                                            >  2:2561650/1‑149 (MQ=255)
                                                                                         ttttttttGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTc                                                           >  1:2301685/1‑150 (MQ=255)
                                                                                         ttttttttGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTc                                                           >  1:4683828/1‑150 (MQ=255)
                                                                                           ttttttGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTACTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCaa                                                         >  1:1355466/1‑150 (MQ=255)
                                                                                               ttGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTTGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCtt                                                     <  1:2819475/150‑1 (MQ=255)
                                                                                                  cTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGaa                                                  >  2:3156826/1‑150 (MQ=255)
                                                                                                   tCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAAt                                                 >  2:674834/1‑150 (MQ=255)
                                                                                                         aCCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAAttttttt                                           <  2:596812/150‑1 (MQ=255)
                                                                                                          ccAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCGCCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTtta                                          <  2:4352182/150‑1 (MQ=255)
                                                                                                          ccAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTtta                                          <  2:4352187/150‑1 (MQ=255)
                                                                                                              gtggtgCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTaa                                      <  1:5004895/150‑1 (MQ=255)
                                                                                                              gtggtgCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTaa                                      <  1:1437509/150‑1 (MQ=255)
                                                                                                                ggtgCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGt                                    <  1:2540341/150‑1 (MQ=255)
                                                                                                                 gtgCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTc                                   >  1:763515/1‑150 (MQ=255)
                                                                                                                  tgCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCa                                  >  1:2539074/1‑150 (MQ=255)
                                                                                                                    cggcggATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAg                                >  1:3448693/1‑150 (MQ=255)
                                                                                                                     ggcggATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCTCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGa                               <  1:357888/150‑1 (MQ=255)
                                                                                                                          aTGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCAc                          <  2:4313766/150‑1 (MQ=255)
                                                                                                                              cACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTg                      >  2:1583864/1‑150 (MQ=255)
                                                                                                                                ctgctgGTACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGtt                    <  2:4973402/150‑1 (MQ=255)
                                                                                                                                  gctgGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGTATTACTTCAATCTTGAATTTTTTTATTATGTCAATATGCACCCTGTTCa                  >  2:4816413/1‑150 (MQ=255)
                                                                                                                                    tgGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGTTCaaa                >  2:2688855/1‑150 (MQ=255)
                                                                                                                                        cttcGCCAGCCTGAAGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGTTCAAATCCt            >  1:170074/4‑150 (MQ=255)
                                                                                                                                         ttcGCCAGCCTGAAGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGTTCAAATCCtt           >  2:4672784/3‑150 (MQ=255)
                                                                                                                                             ccAGCCTGAAGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGTTCAAATCCTTCCCg       >  1:4204789/1‑150 (MQ=255)
                                                                                                                                               aGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGTTCAAATCCTTCCCGGc     <  2:312444/150‑1 (MQ=255)
                                                                                                                                               aGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGTTCAAATCCTTCCCGGc     >  2:3044412/1‑150 (MQ=255)
                                                                                                                                                  ttCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGTTCAAATCCTTCCCGGCGTg  >  1:1292035/1‑150 (MQ=255)
                                                                                                                                                  ttCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGTTCAAATCCTTCACGGCGTg  >  2:4070428/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGTTCAAATCCTTCCCGGCGTG  >  NC_002947/1981057‑1981352

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