Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,411,413 A→T 5.2% intergenic (‑61/+26) yfbS ← / ← hxpA putative transporter YfbS/hexitol phosphatase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,411,4130AT 5.2% 287.8 / 4.1 96intergenic (‑61/+26)yfbS/hxpAputative transporter YfbS/hexitol phosphatase A
Reads supporting (aligned to +/- strand):  ref base A (37/51);  new base T (3/2);  total (41/55)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

GACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCAATGACATG  >  NC_000913/2411266‑2411561
                                                                                                                                                   |                                                                                                                                                    
gACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACg                                                                                                                                                    >  2:1396676/1‑150 (MQ=255)
 aCCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACgg                                                                                                                                                   >  2:1377588/1‑150 (MQ=255)
    gCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGc                                                                                                                                                <  1:401228/150‑1 (MQ=255)
    gCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGc                                                                                                                                                >  1:1046603/1‑150 (MQ=255)
       tCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTaa                                                                                                                                             >  1:729715/1‑150 (MQ=255)
        ccATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTaaa                                                                                                                                            >  2:446220/1‑150 (MQ=255)
         cATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGTCTCGTCGGGGCttttt                                                                                                                                           <  1:1109766/150‑5 (MQ=255)
         cATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGTCTCGTCGGGGCttttt                                                                                                                                           <  1:467177/150‑5 (MQ=255)
         cATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAAt                                                                                                                                           >  1:1524390/1‑150 (MQ=255)
         cATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATAAACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAAt                                                                                                                                           >  1:1580347/1‑150 (MQ=255)
            aCGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCAt                                                                                                                                        >  2:37341/1‑150 (MQ=255)
             cGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAATTCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATg                                                                                                                                       <  2:859093/150‑1 (MQ=255)
             cGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATg                                                                                                                                       <  2:864821/150‑1 (MQ=255)
              gCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGa                                                                                                                                      <  2:1471924/150‑1 (MQ=255)
               cACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGAt                                                                                                                                     >  2:129877/1‑150 (MQ=255)
                aCTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATc                                                                                                                                    >  1:1339687/1‑150 (MQ=255)
                aCTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATc                                                                                                                                    <  1:1190894/150‑1 (MQ=255)
                 ctctGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCa                                                                                                                                   <  1:951273/150‑1 (MQ=255)
                 ctctGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCa                                                                                                                                   >  1:363640/1‑150 (MQ=255)
                    tGCCCGTCGAAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGgt                                                                                                                                <  1:331238/150‑1 (MQ=255)
                     gCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTATTTGAGTGTGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGTCTAAATCATGATCAGgtt                                                                                                                               <  2:12779/150‑2 (MQ=255)
                       ccGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGtt                                                                                                                             >  2:523832/1‑150 (MQ=255)
                         gTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGCGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGCCGGGGCTAAATCATGATCAGGTGTTTc                                                                                                                           <  2:368460/150‑1 (MQ=255)
                             cAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTg                                                                                                                       <  1:1377588/150‑1 (MQ=255)
                             cAAACAAGACGATGG‑AACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGa                                                                                                                      >  1:595485/1‑150 (MQ=255)
                               aaCAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGaa                                                                                                                     >  1:1417298/1‑150 (MQ=255)
                                aCAAGACTATGGCAACGGCCAGTAATGAAAGAACCCAATTCAAATCACGGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGaat                                                                                                                    <  2:171508/150‑1 (MQ=255)
                                aCAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGaat                                                                                                                    <  2:660136/150‑1 (MQ=255)
                                    gACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAataata                                                                                                                >  1:260132/1‑150 (MQ=255)
                                    gACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAataata                                                                                                                >  2:513057/1‑150 (MQ=255)
                                    gACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAataata                                                                                                                <  2:1545516/150‑1 (MQ=255)
                                      cGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCAc                                                                                                                         <  1:993200/139‑1 (MQ=255)
                                      cGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCAc                                                                                                                         >  2:993200/1‑139 (MQ=255)
                                       gATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACa                                                                                                             <  1:541765/150‑1 (MQ=255)
                                          ggCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCg                                                                                                          <  1:1395451/150‑1 (MQ=255)
                                          ggCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCg                                                                                                          >  1:504756/1‑150 (MQ=255)
                                           gCAACGGCCAGTAATTAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGACCCGACGAGACGGGGCTAAATCATTATCAGGTGTTTCACTGAATAATAACATCGc                                                                                                         <  2:216816/150‑1 (MQ=255)
                                           gCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAATTCATGATCAGGTGTTTCACTGAATAATAACATCGc                                                                                                         >  2:993584/1‑150 (MQ=255)
                                           gCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCACCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGc                                                                                                         >  1:1208737/1‑150 (MQ=255)
                                            cAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGTGGATGACTTGATTCTGCCATAAAAAATCCCCGACGAGACGGGGCTAAATCATTATCAGGTGTTTCACTGAATAATAACATCGcc                                                                                                        <  2:432924/150‑1 (MQ=255)
                                            cAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGcc                                                                                                        <  2:1660835/150‑1 (MQ=255)
                                             aaCGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCa                                                                                                       >  2:72306/1‑150 (MQ=255)
                                                   cAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGc                                                                                                 <  1:1305346/150‑1 (MQ=255)
                                                    aGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCt                                                                                                >  2:161319/1‑150 (MQ=255)
                                                       aaTGAAAGAACCCAAATCAATTCACCGTTCACAACTTATGCTTGTTAATTGAGGGGGATGACTTGATTAGGCCATAAAAAAGCCCCGACGAGCCGGGGCTAAATCATGATCAGGTGTTTCACTGAAGAATAACATCGCCAATTGGCTGtt                                                                                             >  2:1252656/1‑150 (MQ=255)
                                                       aaTGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGtt                                                                                             <  2:134519/150‑1 (MQ=255)
                                                       aaTGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTAGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGtt                                                                                             >  2:1501341/1‑150 (MQ=255)
                                                        aTGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGttt                                                                                            <  2:323936/150‑1 (MQ=255)
                                                         tGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTg                                                                                           <  2:1442889/150‑1 (MQ=255)
                                                         tGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTg                                                                                           <  1:542801/150‑1 (MQ=255)
                                                            aaGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCTTTTGGCTGTTTGGTc                                                                                        <  2:597094/150‑1 (MQ=255)
                                                            aaGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTc                                                                                        <  2:486682/150‑1 (MQ=255)
                                                             aGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTca                                                                                       <  2:457217/150‑1 (MQ=255)
                                                               aaCCCAAATTAATTCATCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTcaca                                                                                     <  2:709315/150‑1 (MQ=255)
                                                                aCCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAg                                                                                    >  1:528427/1‑150 (MQ=255)
                                                                aCCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAg                                                                                    <  2:1339687/150‑1 (MQ=255)
                                                                 cccAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGt                                                                                   >  1:1480381/1‑150 (MQ=255)
                                                                      aTCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAAttt                                                                              <  1:94800/150‑1 (MQ=255)
                                                                      aTCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAAttt                                                                              <  1:129877/150‑1 (MQ=255)
                                                                       tCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTg                                                                             <  2:490797/150‑1 (MQ=255)
                                                                       tCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTg                                                                             <  2:1046603/150‑1 (MQ=255)
                                                                          aTTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTc                                                                          <  2:154063/150‑1 (MQ=255)
                                                                            tCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGTCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCa                                                                        >  2:1661343/1‑150 (MQ=255)
                                                                             cACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGt                                                                                         <  2:802214/132‑1 (MQ=255)
                                                                             cACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGt                                                                                         >  1:802214/1‑132 (MQ=255)
                                                                                   tCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACtgtg                                                                 >  2:602670/1‑150 (MQ=255)
                                                                                    cacaACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACCTCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTgg                                                                >  2:560670/1‑150 (MQ=255)
                                                                                    cacaACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTgg                                                                >  2:876056/1‑150 (MQ=255)
                                                                                      caACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGGTCCAGCCTGTgggg                                                              >  2:312530/1‑148 (MQ=255)
                                                                                      caACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTggag                                                              >  2:1200343/1‑150 (MQ=255)
                                                                                       aaCTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTggagg                                                             <  2:197180/150‑1 (MQ=255)
                                                                                         cTTATCCTTGTTAATTGAGGGTGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGAc                                                           <  2:575603/150‑1 (MQ=255)
                                                                                         cTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGAc                                                           <  2:1399681/150‑1 (MQ=255)
                                                                                                ttGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCg                                                    <  2:1583960/150‑1 (MQ=255)
                                                                                                 tGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGa                                                   <  2:1315824/150‑1 (MQ=255)
                                                                                                   ttAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGAcc                                                 <  1:1501341/150‑1 (MQ=255)
                                                                                                   ttAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGAcc                                                 <  2:1580347/150‑1 (MQ=255)
                                                                                                       ttGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCAt                                             <  2:1417298/150‑1 (MQ=255)
                                                                                                       ttGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCAt                                             <  1:72306/150‑1 (MQ=255)
                                                                                                        tGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCAtt                                            <  2:528427/150‑1 (MQ=255)
                                                                                                         gcgggggTTGAGTTTATTCTTCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTGCACTGAATAATAACATCGCCTTTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCTTTc                                           <  2:1304782/148‑1 (MQ=255)
                                                                                                         gAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTc                                           <  1:513057/150‑1 (MQ=255)
                                                                                                          aGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCa                                          >  1:76308/1‑150 (MQ=255)
                                                                                                           gggggATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCCCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAg                                         <  1:446220/150‑1 (MQ=255)
                                                                                                             gggATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGgc                                       >  1:308586/1‑150 (MQ=255)
                                                                                                              ggATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGgcg                                      >  2:346866/1‑150 (MQ=255)
                                                                                                                 tGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCgg                                   <  1:921087/150‑1 (MQ=255)
                                                                                                                   aCTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCgggg                                 >  2:148993/1‑150 (MQ=255)
                                                                                                                    cTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGt                                >  1:515219/1‑150 (MQ=255)
                                                                                                                         ttCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTg                           <  2:1524390/150‑1 (MQ=255)
                                                                                                                                 tgatttaGCCCCGTCTCGTCGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCg                   >  1:610268/7‑150 (MQ=255)
                                                                                                                                  tatttaGCCCCGTCTCGTCGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGt                  >  1:1379508/6‑150 (MQ=255)
                                                                                                                                    tttaGCCCCGTCTCGTCGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTa                >  1:1393920/4‑150 (MQ=255)
                                                                                                                                    aaaaGCCCCGCCGAGCCGGCGCTAAATCATGATCAGGTGTTTCACTGAAAAATAACATCGCCTTTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTa                <  2:332818/150‑1 (MQ=255)
                                                                                                                                       aGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAAcc             <  2:504756/150‑1 (MQ=255)
                                                                                                                                       aGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAAcc             <  1:602670/150‑1 (MQ=255)
                                                                                                                                                gagaCGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCAatgaca    <  2:260132/150‑1 (MQ=255)
                                                                                                                                                 agaCGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCAatgacat   >  2:598213/1‑150 (MQ=255)
                                                                                                                                                  gaCGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCAatgacatg  <  1:1396676/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTTCACTGAATAATAACATCGCCATTTGGCTGTTTGGTCACAGTAATTTGCTCCAGACTGTGGAGGACCAAATCGACCTCATTCAGGCGCGGGGTATCTGCCGGAGCGTTAACCGCAATGACATG  >  NC_000913/2411266‑2411561

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