Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,981,135 T→G 11.5% intergenic (+31/‑269) 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,1350TG11.5% 204.0 / 19.3 87intergenic (+31/‑269)rpsA/ihfB30S ribosomal protein S1/integration host factor subunit beta
Reads supporting (aligned to +/- strand):  ref base T (36/39);  new base G (4/6);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

AATCGAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCT  >  NC_002947/1980986‑1981284
                                                                                                                                                     |                                                                                                                                                       
atcggtaaGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGCAGCAATGGCCAAGCAGCACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCt                                                                                                                                                         >  1:1264389/7‑150 (MQ=255)
   cGAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGNCAAGCANAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTNgg                                                                                                                                                      <  1:1658412/150‑1 (MQ=255)
   cGAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCNCgg                                                                                                                                                      <  2:1510990/150‑1 (MQ=255)
    gAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCggg                                                                                                                                                     >  1:2023000/1‑150 (MQ=255)
     aaaGACGACGCTGGAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGt                                                                                                                                                    >  1:1194917/1‑150 (MQ=255)
      aaGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTc                                                                                                                                                   <  1:707925/150‑1 (MQ=255)
         acgacgCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGcc                                                                                                                                                <  2:2023000/150‑1 (MQ=255)
          cgacgCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCGCCGGGTCGccc                                                                                                                                               <  2:409309/150‑1 (MQ=255)
          cgacgCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCGAGGGTCGccc                                                                                                                                               <  2:2208254/150‑1 (MQ=255)
           gacgCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTACGGAAGCTGCTGCCGACAACACCCTGGCTGCGCTGTTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCGAGGGTCGCCCt                                                                                                                                             <  1:696959/150‑1 (MQ=255)
                 gAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCGAGGGTCGCCCtttttac                                                                                                                                       <  1:821787/150‑3 (MQ=255)
                         cgAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCcgc                                                                                                                               <  2:1122512/150‑1 (MQ=255)
                             gCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCcgc                                                                                                                               <  2:1612090/146‑1 (MQ=255)
                             gCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCcgc                                                                                                                               >  1:1612090/1‑146 (MQ=255)
                             gCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGAcc                                                                                                                           >  1:2216440/1‑150 (MQ=255)
                             gCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGAcc                                                                                                                           >  2:1254954/1‑150 (MQ=255)
                               tATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAACCGCCAAGCAGAACTGAGTTCTCTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAg                                                                                                                         <  2:1261866/150‑1 (MQ=255)
                                 tCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGccc                                                                                                                                               >  1:676034/1‑127 (MQ=255)
                                 tCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGccc                                                                                                                                               <  2:676034/127‑1 (MQ=255)
                                    agagCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAggtggt                                                                                                                    >  1:1538374/1‑150 (MQ=255)
                                    agagCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAggtggt                                                                                                                    >  1:2070155/1‑150 (MQ=255)
                                       gCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATTGGTAAAAAGGGCGACCCTCGGGTCGCCCtttttt                                                                                                                                        >  1:958506/1‑127 (MQ=255)
                                       gCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATTGGTAAAAAGGGCGACCCTCGGGTCGCCCtttttt                                                                                                                                        <  2:958506/127‑1 (MQ=255)
                                         cTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTgcggc                                                                                                               <  1:1254954/150‑1 (MQ=255)
                                          tGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGCCCAGGTGGTgcggcg                                                                                                              >  2:1359319/1‑150 (MQ=255)
                                            aaaaaCGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCCGGTGGTGCGGCGGa                                                                                                            <  2:1913503/150‑1 (MQ=255)
                                              aaaCGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATg                                                                                                          >  1:1094640/1‑150 (MQ=255)
                                                    tCCGGAAGCTGCTGCCGACACCACTATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCActg                                                                                                   >  2:1622232/1‑150 (MQ=255)
                                                     ccGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCActgc                                                                                                  <  1:1938605/150‑1 (MQ=255)
                                                       ggAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCActgctg                                                                                                <  2:660203/150‑1 (MQ=255)
                                                               ctgcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGc                                                                                        >  1:1017657/1‑150 (MQ=255)
                                                                tgcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCGAGGGTCGCCCt                                                                                                                                             <  2:1199178/97‑1 (MQ=255)
                                                                tgcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCGAGGGTCGCCCt                                                                                                                                             >  1:1199178/1‑97 (MQ=255)
                                                                 gcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCa                                                                                      <  2:2070155/150‑1 (MQ=255)
                                                                    gACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTc                                                                                                                                                   <  1:2314096/88‑1 (MQ=255)
                                                                    gACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTc                                                                                                                                                   <  1:2314688/88‑1 (MQ=255)
                                                                    gACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTc                                                                                                                                                   >  2:2314096/1‑88 (MQ=255)
                                                                    gACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTc                                                                                                                                                   >  2:2314688/1‑88 (MQ=255)
                                                                    gACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCt                                                                                   <  1:718778/150‑1 (MQ=255)
                                                                      caccacCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTc                                                                                 <  2:2216440/150‑1 (MQ=255)
                                                                       accaccATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCGAGGGTCGccc                                                                                                                                               <  1:2166043/89‑1 (MQ=255)
                                                                       accaccATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCGAGGGTCGccc                                                                                                                                               >  2:2166043/1‑89 (MQ=255)
                                                                        ccaccaGGGCTGCGCCGCTGAGCGAAGCAATGGCCAAGCACAACTGAGTTCTGTTTGACCGGTAAAAACGGCGACCCTCAGCTAGGCTTTCTTTTTGCTCCGCGACCAGGTGGTGCGAGGATGCGCACTGCTGGAACACGCCAGCTTCAg                                                                               <  1:663265/150‑1 (MQ=255)
                                                                         caccaTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGACCACGCCAGCTTCAgg                                                                              >  2:140168/1‑150 (MQ=255)
                                                                          accaTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAac                                                                                            <  2:1929030/135‑1 (MQ=255)
                                                                          accaTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAac                                                                                            >  1:1929030/1‑135 (MQ=255)
                                                                          accaTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGc                                                                             >  1:289207/1‑150 (MQ=255)
                                                                                 ctgcgctgcTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGaa                                                                      <  1:1258029/150‑1 (MQ=255)
                                                                                  tgcgctgcTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCcgc                                                                                                                               >  2:1548267/1‑93 (MQ=255)
                                                                                  tgcgctgcTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCcgc                                                                                                                               <  1:1548267/93‑1 (MQ=255)
                                                                                     gctgctGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGcc                                                                   >  2:2020736/1‑150 (MQ=255)
                                                                                      ctgctgCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCg                                                                 >  2:1444240/1‑150 (MQ=255)
                                                                                        gctgcGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGtt                                                               <  1:1121963/150‑1 (MQ=255)
                                                                                        gctgcGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGtt                                                               <  2:1445785/150‑1 (MQ=255)
                                                                                        gctgcGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGtt                                                               <  1:1359319/150‑1 (MQ=255)
                                                                                             gcgAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAgg                                                          >  1:409853/1‑150 (MQ=255)
                                                                                             gcgAAGCAATGGCCAAGCAGAACCGAGTCCTGGTGGATTGATCACTCTACCTTTCCCCTTTTCGCCCTTTTTTTTTGCCCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAg                                                           <  2:62189/150‑1 (MQ=255)
                                                                                                aaGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGTCGACCCGAGGGTCGCCCttt                                                                                                                                           >  1:1892921/1‑67 (MQ=37)
                                                                                                aaGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGTCGACCCGAGGGTCGCCCttt                                                                                                                                           <  2:1892921/67‑1 (MQ=37)
                                                                                                  gCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCgg                                                     <  1:933956/150‑1 (MQ=255)
                                                                                                    aaTGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGaa                                                                      <  2:2222060/131‑1 (MQ=255)
                                                                                                    aaTGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGaa                                                                      >  1:2222060/1‑131 (MQ=255)
                                                                                                    aaTGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGt                                                   >  1:2298768/1‑150 (MQ=255)
                                                                                                        gCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAc                                               <  1:1622232/150‑1 (MQ=255)
                                                                                                          cAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAccc                                             >  1:1189115/1‑150 (MQ=255)
                                                                                                            aGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCgg                                                     <  1:2326855/140‑1 (MQ=255)
                                                                                                            aGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCgg                                                     >  2:2326855/1‑140 (MQ=255)
                                                                                                               aGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGaa                                        <  2:1189115/150‑1 (MQ=255)
                                                                                                               aGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGaa                                        <  2:850181/150‑1 (MQ=255)
                                                                                                                gAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCgg                                                     >  2:1495962/1‑136 (MQ=255)
                                                                                                                gAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCgg                                                     <  1:1495962/136‑1 (MQ=255)
                                                                                                                    tGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAgg                                   >  2:2063472/1‑150 (MQ=255)
                                                                                                                           tGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACt                            >  1:1349280/1‑150 (MQ=255)
                                                                                                                             tttGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCCCAACAACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGCCGAAGCCCGGTGCAGGAGCGGGTCCACCCGCGAGGCAGGCGCCagcgg                          >  2:930580/1‑146 (MQ=255)
                                                                                                                             tttGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCCTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTgg                          >  2:1276593/1‑150 (MQ=255)
                                                                                                                              ttGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGCGCGGGTTCACCCGCGAAGCAGGCGCCACTGGc                         <  2:1349280/150‑1 (MQ=255)
                                                                                                                               tGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGcc                        <  2:1538374/150‑1 (MQ=255)
                                                                                                                                gATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGc                                                                             <  1:1366593/96‑1 (MQ=255)
                                                                                                                                gATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGc                                                                             >  2:1366593/1‑96 (MQ=255)
                                                                                                                                     gTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTa                  >  2:2126557/1‑150 (MQ=255)
                                                                                                                                        aaaaGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGc               <  1:1444240/150‑1 (MQ=255)
                                                                                                                                         aaaGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCggg                                                    <  2:98621/112‑1 (MQ=255)
                                                                                                                                         aaaGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCggg                                                    >  1:98621/1‑112 (MQ=255)
                                                                                                                                          aaGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCgg             <  2:1017657/150‑1 (MQ=255)
                                                                                                                                           aGGGCGCCCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGt            <  1:671968/150‑1 (MQ=255)
                                                                                                                                            gggCGACCCTCGGGTCGCCCTTTTTTTTTTTTCTCTTCCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTa           >  2:410743/1‑150 (MQ=255)
                                                                                                                                            gggCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTa           >  1:1577044/1‑150 (MQ=255)
                                                                                                                                             ggCGACCCGAGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGCACGGGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAg          >  1:1329292/1‑150 (MQ=255)
                                                                                                                                               cGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGAc        <  1:1081706/150‑1 (MQ=255)
                                                                                                                                                     tCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCt  >  2:1473586/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
AATCGAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCT  >  NC_002947/1980986‑1981284

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