Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,710,708 T→C 16.2% intergenic (‑233/+91) eptB ← / ← yhjX KDO phosphoethanolamine transferase, Ca(2+)‑inducible/pyruvate‑inducible inner membrane protein, putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,710,7080TC16.2% 224.1 / 35.2 99intergenic (‑233/+91)eptB/yhjXKDO phosphoethanolamine transferase, Ca(2+)‑inducible/pyruvate‑inducible inner membrane protein, putative transporter
Reads supporting (aligned to +/- strand):  ref base T (37/46);  new base C (8/8);  total (45/54)
Fisher's exact test for biased strand distribution p-value = 7.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

Only 100 of 101 total aligned reads displayed.

AGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGA  >  NC_000913/3710610‑3710807
                                                                                                  |                                                                                                   
aGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAAtt                                                                                                    <  2:2541296/100‑1 (MQ=255)
     ggggAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGccc                                                                                               <  2:2155842/100‑1 (MQ=255)
     ggggAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAATTAACTGGccc                                                                                               <  2:304061/100‑1 (MQ=255)
       ggAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCaa                                                                                             >  2:352099/1‑100 (MQ=255)
       ggAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCaa                                                                                             >  2:1768364/1‑100 (MQ=255)
          aaGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAATTAACTGGcccttttt                                                                                          <  2:1567534/100‑6 (MQ=255)
          aaGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAATTAACTGGcccttttt                                                                                          <  2:238788/100‑6 (MQ=255)
           agagTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAATTAACTGGcccttttta                                                                                         <  2:2538925/100‑7 (MQ=255)
            gagTGTGCATGAAGCACAAATTTACAAAAATGCAAATTTTAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAATTAACTGGccctttttat                                                                                        <  1:1630472/100‑8 (MQ=255)
              gtgtGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTAc                                                                                      <  2:219072/100‑1 (MQ=255)
              gtgtGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTAc                                                                                      >  2:2118160/1‑100 (MQ=255)
              gtgtGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAATTAACTGGccctttttactg                                                                                      <  2:589791/100‑10 (MQ=38)
              gtgtGCATGAAGCACAAAAAAACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAATTAACTGGccctttttactg                                                                                      <  1:1636427/100‑10 (MQ=37)
               tgtgCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAATTAACTGGccctttttactgc                                                                                     <  2:82976/100‑11 (MQ=38)
               tgtgCATGAAGCAAAAAATAACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACg                                                                                     <  1:2491174/100‑1 (MQ=255)
                gtgCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGa                                                                                    >  1:2319014/1‑100 (MQ=255)
                 tgCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAg                                                                                   >  1:385193/1‑100 (MQ=255)
                  gCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGt                                                                                  <  2:573889/100‑1 (MQ=255)
                   cATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTc                                                                                 >  2:2308303/1‑100 (MQ=255)
                   cATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTc                                                                                 <  2:1281734/100‑1 (MQ=255)
                   cATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTc                                                                                 <  1:532167/100‑1 (MQ=255)
                    aTGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCa                                                                                >  2:1781422/1‑100 (MQ=255)
                     tGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCAt                                                                               >  1:2059958/1‑100 (MQ=255)
                      gAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGAAGTAAAAAGGGCCAGTTAATTGGCCAAAATTTTACGAGTCAtt                                                                              >  1:461510/1‑100 (MQ=255)
                       aaGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTg                                                                             <  2:1702548/100‑1 (MQ=255)
                       aaGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTg                                                                             <  1:160591/100‑1 (MQ=255)
                         gCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGtt                                                                           <  2:879916/100‑1 (MQ=255)
                           acaAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTta                                                                         <  1:1941732/100‑1 (MQ=255)
                            caAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTtat                                                                        >  1:1732153/1‑100 (MQ=255)
                                tttACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACAt                                                                    >  1:185867/1‑100 (MQ=255)
                                    cAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCg                                                                >  2:1783638/1‑100 (MQ=255)
                                    cAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCg                                                                >  1:1221401/1‑100 (MQ=255)
                                     aaaaaTGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGt                                                               <  2:2573463/100‑1 (MQ=255)
                                      aaaaTGTAAGTTTGAAAAGAAACTTAAGCGTTACAAAAATCAGTAAAAAGGTCCAGTTAGTTGGCCCAAATTGTACGAGTCATTGTTTATACATTACg                                                                <  1:2265883/98‑1 (MQ=255)
                                      aaaaTGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGtt                                                              >  2:253711/1‑100 (MQ=255)
                                      aaaaTGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGtt                                                              >  1:2067130/1‑100 (MQ=255)
                                       agaTGGAAGTTGAAAGAAAGACTTTAGCGTGACAAAAATCAGTAAAAATGGCCAGTTAATTGGCCCAATTTTTACGAGTCATTGTTTATACATTTCGTTg                                                             <  1:1183348/98‑1 (MQ=255)
                                       aaaTGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGttn                                                             <  2:1652785/100‑2 (MQ=255)
                                        aaTGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGa                                                            <  2:1608058/100‑1 (MQ=255)
                                          tGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCGAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATg                                                          >  2:392720/1‑100 (MQ=255)
                                            cAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCg                                                        >  1:1629361/1‑100 (MQ=255)
                                            cAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTAATTGTTTATACATTTCGTTGATGCg                                                        >  1:663350/1‑100 (MQ=255)
                                            cAAGTTTGAAAAGAGACCTTAGCGTTACAAAAAGCAATAAAAAAGGCGAGGTAATTGGGCCAACTTATACGAGTCAATGTTTATACATTTAGCTGATGCg                                                        >  1:790956/1‑100 (MQ=255)
                                             aaGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGt                                                       >  2:204967/1‑100 (MQ=255)
                                               gTTTGAACAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGttt                                                     <  1:1392473/100‑1 (MQ=255)
                                               gTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGttt                                                     <  2:1721973/100‑1 (MQ=255)
                                                tttGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTg                                                    <  2:325616/100‑1 (MQ=255)
                                                  tGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCGAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGAt                                                  <  2:855959/100‑1 (MQ=255)
                                                   gAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATg                                                 <  2:2141817/100‑1 (MQ=255)
                                                   gAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATg                                                 <  1:460146/100‑1 (MQ=255)
                                                     aaaGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGca                                               >  2:1890814/1‑100 (MQ=255)
                                                       agagACTTTATCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATGGTCTATACATTTCGTTGATGCGTTTGATGcaca                                             <  1:1702713/100‑1 (MQ=255)
                                                       agagACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGcaca                                             >  2:1413631/1‑100 (MQ=255)
                                                       agagACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGcaca                                             <  2:1916886/100‑1 (MQ=255)
                                                          gaCTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACgg                                          >  2:482498/1‑100 (MQ=255)
                                                           aCTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGa                                         <  1:343874/100‑1 (MQ=255)
                                                           aCTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGa                                         <  1:482496/100‑1 (MQ=255)
                                                             tttAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTGACGAGTCATTGTTTATACATGTCGTTGATGCGTTTGATGCACACGGAAg                                       <  2:1732857/100‑1 (MQ=255)
                                                              ttAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGc                                      <  2:573533/100‑1 (MQ=255)
                                                              ttAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGc                                      <  2:543662/100‑1 (MQ=255)
                                                                  cGTTACAAAAAGCAGTAAAAAGGGCCATTTAATTGGACCAAATTTTACGAGTCAGTGTTTATACATTTCGTTGATGCGTTTTATGCACACGGAAg                                       <  1:1447762/95‑1 (MQ=255)
                                                                   gTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAg                                 >  2:452074/1‑100 (MQ=255)
                                                                   gTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAg                                 <  2:1923688/100‑1 (MQ=255)
                                                                    ttACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGc                                <  2:1011168/100‑1 (MQ=255)
                                                                    ttACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGc                                >  1:730197/1‑100 (MQ=255)
                                                                     tACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGATGCTGAAGc                                >  1:646061/1‑99 (MQ=255)
                                                                      aCAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGccc                              >  1:548464/1‑100 (MQ=255)
                                                                      aCAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGccc                              >  1:2571985/1‑100 (MQ=255)
                                                                       cAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCa                             >  2:105513/1‑100 (MQ=255)
                                                                        aaaaaGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAg                            <  2:2367865/100‑1 (MQ=255)
                                                                        aaaaaGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAg                            <  1:766616/100‑1 (MQ=255)
                                                                         aaaaGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGt                           <  1:1138227/100‑1 (MQ=255)
                                                                            aGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGc                        <  2:878526/100‑1 (MQ=255)
                                                                            aGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATACACACGGAAGCTGAAGCCCAGTAGc                        <  1:60741/100‑1 (MQ=255)
                                                                             gCAGTAAAAAGGGCCAATTAATTGGCCCAACTTTTACGAGTCATTGTTTATACATTTCGTTGATGCTTTTGCTGCACACGGAAGCTGAAGGCCAg                            >  1:2007124/1‑95 (MQ=255)
                                                                              gtaaaatttgggCCAATTAACTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTc                      >  2:2011026/10‑100 (MQ=38)
                                                                              gtaaaatttgggCCAATTAACTGGCCCAAATTTTACGAGTCATTGAAAATACATTTCGTTGATGCGTTTGATGCACACGGAAGATGAAGCCCAGTAGCTc                      >  1:465446/10‑100 (MQ=37)
                                                                              gtaaaatttgggCCAATTAACTGGCACAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCCCACGGAAGCTGAAGCCCAGTAGCTc                      >  2:875337/10‑100 (MQ=255)
                                                                               ngTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGCTGCGTTGGATGCACACGGAAGCTGAAGCCCAGTAGCTcg                     >  1:1818925/2‑100 (MQ=255)
                                                                                gTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTcgc                    >  1:1371947/1‑100 (MQ=255)
                                                                                aaaatttgggCCAATTAACTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTcgc                    >  1:786042/8‑100 (MQ=255)
                                                                                 aaatttgggCCAATTAACTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTcgcg                   >  2:1398824/7‑100 (MQ=255)
                                                                                  aatttgggCCAATTAACTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCgg                  >  1:744088/6‑100 (MQ=255)
                                                                                   aaaaGGGCCAGTTTATTGGCCCAAATTTTACGAGCCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGc                 <  1:608974/100‑1 (MQ=255)
                                                                                   aaaaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGc                 >  2:1936309/1‑100 (MQ=255)
                                                                                   aaaaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGc                 <  2:1132954/100‑1 (MQ=255)
                                                                                     ttgggCCAATTAACTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTg               >  1:422587/3‑100 (MQ=255)
                                                                                     ttgggCCAATTAACTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTg               >  2:133212/3‑100 (MQ=255)
                                                                                     aaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTg               <  2:239920/100‑1 (MQ=255)
                                                                                     aaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTg               <  1:738815/100‑1 (MQ=255)
                                                                                       tccgatcttAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAg             <  1:625097/93‑1 (MQ=255)
                                                                                       gggCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAg             <  2:1134351/100‑1 (MQ=255)
                                                                                       gggCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAg             <  2:2067143/100‑1 (MQ=255)
                                                                                        ggCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGc            <  2:1044048/100‑1 (MQ=255)
                                                                                        ggCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGc            >  2:1493792/1‑100 (MQ=255)
                                                                                        ggCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGc            <  1:352097/100‑1 (MQ=255)
                                                                                             gTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTaa       <  2:2440472/100‑1 (MQ=255)
                                                                                               tAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAgagatcgga     >  2:625100/1‑92 (MQ=255)
                                                                                               tAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAg     >  2:2021602/1‑100 (MQ=255)
                                                                                                  ttGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGa  <  2:2158545/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
AGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGA  >  NC_000913/3710610‑3710807

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