Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 109,644 A→G T456A (ACT→GCT)  secA → preprotein translocase subunit, ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913109,6440AG98.8% 298.8 / ‑4.2 83T456A (ACT→GCT) secApreprotein translocase subunit, ATPase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (39/43);  total (40/43)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01

GCTGGATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTACTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACT  >  NC_000913/109496‑109781
                                                                                                                                                    |                                                                                                                                         
gCTGGATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCt                                                                                                                                         >  1:261783/1‑151 (MQ=255)
   ggATACCGTGGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTAtc                                                                                                                                      >  2:1272635/1‑151 (MQ=255)
   ggATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTa                                                                                                                                        <  1:66822/149‑1 (MQ=255)
     aTACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTAtctc                                                                                                                                    <  1:1065547/151‑1 (MQ=255)
         cgtcgtTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATc                                                                                                                                <  1:601597/151‑1 (MQ=255)
          gtcgtTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCg                                                                                                                               >  2:1157349/1‑151 (MQ=255)
                 ttCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATc                                                                                                                        >  2:127130/1‑151 (MQ=255)
                   ccGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCg                                                                                                                       <  2:1312150/150‑1 (MQ=255)
                    cGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCgg                                                                                                                      >  1:945638/1‑150 (MQ=255)
                    cGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGa                                                                                                                     <  2:715331/151‑1 (MQ=255)
                      aCCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGc                                                                                                                   >  1:615470/1‑151 (MQ=255)
                       ccAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCt                                                                                                                  >  1:539885/1‑151 (MQ=255)
                        cAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTg                                                                                                                 <  1:191383/151‑1 (MQ=255)
                         aaCCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTgg                                                                                                                <  1:167375/151‑1 (MQ=255)
                           ccGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGgtg                                                                                                              >  2:1209364/1‑151 (MQ=255)
                           ccGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGgtg                                                                                                              <  2:280311/151‑1 (MQ=255)
                           ccGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGgtg                                                                                                              <  2:278332/151‑1 (MQ=255)
                                  aTGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACg                                                                                                       >  1:883587/1‑151 (MQ=255)
                                     aTTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAAc                                                                                                    <  1:13125/151‑1 (MQ=255)
                                     aTTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAAc                                                                                                    <  1:356976/151‑1 (MQ=255)
                                        cGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGa                                                                                                 <  1:1122660/151‑1 (MQ=255)
                                         gTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGAGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGAc                                                                                                >  2:1073726/1‑151 (MQ=255)
                                             aGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCaa                                                                                             >  1:620516/1‑150 (MQ=255)
                                                   gCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGcc                                                                                         <  1:545118/148‑1 (MQ=255)
                                                   gCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCgg                                                                                       <  1:1172233/150‑1 (MQ=255)
                                                   gCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCgg                                                                                       <  1:1256542/150‑1 (MQ=255)
                                                       gACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTa                                                                                  <  2:880551/151‑1 (MQ=255)
                                                       gACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTa                                                                                  <  1:1011092/151‑1 (MQ=255)
                                                         ccTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAg                                                                                <  2:363449/151‑1 (MQ=255)
                                                            ggTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGcac                                                                             >  1:1033965/1‑151 (MQ=255)
                                                             gTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGcaca                                                                            <  2:458550/151‑1 (MQ=255)
                                                                 aCATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAAc                                                                          <  2:57708/149‑1 (MQ=255)
                                                                 aCATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACg                                                                         <  1:316906/150‑1 (MQ=255)
                                                                 aCATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGt                                                                        >  2:129544/1‑151 (MQ=255)
                                                                  cATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTc                                                                       <  2:332929/151‑1 (MQ=255)
                                                                   aTGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTcc                                                                      >  1:1468779/1‑151 (MQ=255)
                                                                    tGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCt                                                                     >  2:1335429/1‑151 (MQ=255)
                                                                     gACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCt                                                                     <  2:141974/150‑1 (MQ=255)
                                                                     gACTGAAGAGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGGGCTATCGCCATCGAAAAATCGGTGCTGGTGTCAAACGTACTGACCAAAGCCGGTATTAAGCACAACGTc                                                                       >  2:241230/1‑148 (MQ=255)
                                                                      aCTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGa                                                                   <  2:1225083/151‑1 (MQ=255)
                                                                         gAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGaa                                                                  <  1:341934/149‑1 (MQ=255)
                                                                         gAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGaa                                                                  <  2:909676/149‑1 (MQ=255)
                                                                         gAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACg                                                                <  2:1322205/151‑1 (MQ=255)
                                                                         gAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACg                                                                >  1:255950/1‑151 (MQ=255)
                                                                          aaGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGc                                                               <  2:239354/151‑1 (MQ=255)
                                                                           aGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGcc                                                              <  1:809896/151‑1 (MQ=255)
                                                                            gCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCa                                                             <  2:1252699/151‑1 (MQ=255)
                                                                              ggAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCaaa                                                           <  2:261755/151‑1 (MQ=255)
                                                                                aaaaaaTTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAAtt                                                         <  1:1124441/151‑1 (MQ=255)
                                                                                   aaaTTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAACCGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTc                                                        >  2:1157703/1‑149 (MQ=255)
                                                                                    aaTTCAGGCGATCATTGAAGATATCAAAGAACGTACTTCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGGGTCAAACGAACTGATCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCAc                                                     >  2:183540/1‑151 (MQ=255)
                                                                                       tCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGcc                                                  >  1:638064/1‑151 (MQ=255)
                                                                                          ggCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAAc                                               <  1:1320616/151‑1 (MQ=255)
                                                                                           gCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACg                                              >  2:433294/1‑151 (MQ=255)
                                                                                            cGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACg                                              >  2:834249/1‑150 (MQ=255)
                                                                                            cGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTACTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACg                                              >  2:1358585/1‑150 (MQ=255)
                                                                                             gATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGaa                                            <  2:501824/151‑1 (MQ=255)
                                                                                             gATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGaa                                            <  1:61577/151‑1 (MQ=255)
                                                                                                cATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCCTCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAgcg                                         >  2:42186/1‑151 (MQ=255)
                                                                                                cATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGAGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAgcg                                         >  1:456196/1‑151 (MQ=255)
                                                                                                 aTTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAgcgg                                        >  2:1130785/1‑151 (MQ=255)
                                                                                                   tGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAgcggcg                                      >  2:604951/1‑151 (MQ=255)
                                                                                                    gAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAgcggcg                                      <  1:1247439/150‑1 (MQ=255)
                                                                                                     aaGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAGCGGCGAt                                    >  2:366569/1‑151 (MQ=255)
                                                                                                       gATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAGCGGCGAttg                                  >  2:698592/1‑151 (MQ=255)
                                                                                                       gATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAGCGGCGAttg                                  >  2:1117958/1‑151 (MQ=255)
                                                                                                        atatCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAGCGGCGAttgt                                 >  1:547279/1‑151 (MQ=255)
                                                                                                         tatCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAGCGGCGAttgtt                                >  2:154694/1‑151 (MQ=255)
                                                                                                         tatCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAGCGGCGAttgtt                                >  2:182405/1‑151 (MQ=255)
                                                                                                          atCAAAGATCGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAGCGGCGAttgtt                                <  1:163128/150‑1 (MQ=255)
                                                                                                          atCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAGCGGCGAttgttg                               >  2:1169662/1‑151 (MQ=255)
                                                                                                            cAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAGCGGCGATTGTTGc                              <  2:447282/150‑1 (MQ=255)
                                                                                                               aGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAGCGGCGATTGTTGCTcag                          <  2:1043726/151‑1 (MQ=255)
                                                                                                                gAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAGCGGCGATTGTTGCTcagg                         <  1:1258545/151‑1 (MQ=255)
                                                                                                                  aCGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAGCGGCGATTGTTGCTcaggc                        >  1:829084/1‑150 (MQ=255)
                                                                                                                      aCTGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGtt                   >  2:111662/1‑151 (MQ=255)
                                                                                                                        tGCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTAt                 >  2:47909/1‑151 (MQ=255)
                                                                                                                         gCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTa                  <  1:1073827/149‑1 (MQ=255)
                                                                                                                         gCGAAAGGCCAGCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATc                >  1:208248/1‑151 (MQ=255)
                                                                                                                                    gCCGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTg     <  1:794068/151‑1 (MQ=255)
                                                                                                                                     ccGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTg     >  1:1012841/1‑150 (MQ=255)
                                                                                                                                      cGGTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGAc   >  2:92568/1‑151 (MQ=255)
                                                                                                                                       ggTGCTGGTGGGTGCTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACt  <  1:1354464/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                         
GCTGGATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGTGGGTACTATCTCCATCGAAAAATCGGAGCTGGTGTCAAACGAACTGACCAAAGCCGGTATTAAGCACAACGTCCTGAACGCCAAATTCCACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACT  >  NC_000913/109496‑109781

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