Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,046,784 A→T 6.6% intergenic (‑182/‑182) yihG ← / → polA inner membrane protein, inner membrane acyltransferase/5' to 3' DNA polymerase and 3' to 5'/5' to 3' exonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,046,7840AT 6.6% 531.2 / 19.2 150intergenic (‑182/‑182)yihG/polAinner membrane protein, inner membrane acyltransferase/5' to 3' DNA polymerase and 3' to 5'/5' to 3' exonuclease
Reads supporting (aligned to +/- strand):  ref base A (69/71);  new base T (2/8);  total (71/79)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 160 total aligned reads displayed.

TTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCTGTGGCATC  >  NC_000913/4046636‑4046930
                                                                                                                                                    |                                                                                                                                                  
tttAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTTCATGAAAAACAAATACAAACATAATGCTGTaa                                                                                                                                                   >  2:2659996/1‑150 (MQ=255)
 ttAGCATAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTaaa                                                                                                                                                  <  1:767275/150‑1 (MQ=255)
 ttAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTaaa                                                                                                                                                  <  2:2213095/150‑1 (MQ=255)
 ttAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTaaa                                                                                                                                                  <  2:2213268/150‑1 (MQ=255)
  tAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTaaaa                                                                                                                                                 <  2:108739/150‑1 (MQ=255)
   aGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTaaaaa                                                                                                                                                <  1:1096642/150‑1 (MQ=255)
   aGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTaaaaa                                                                                                                                                <  1:512141/150‑1 (MQ=255)
    gCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAg                                                                                                                                               <  1:2303638/150‑1 (MQ=255)
      aGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCa                                                                                                                                             <  2:2714904/150‑1 (MQ=255)
          aaCGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGat                                                                                                                                         >  2:1687221/1‑150 (MQ=255)
          aaCGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGat                                                                                                                                         >  1:137811/1‑150 (MQ=255)
            cGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTTTTTACAGCATGatgt                                                                                                                                       <  2:2213011/150‑3 (MQ=255)
             gACCCCCATAACGGATGATCCGGACGGAGAAACAGACTTCAGAACTATCCACCGTAGAAACAATCCCATGATCTCAAGGATTAGGTATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTTTTTACAGCATGatgtt                                                                                                                                      <  2:2418186/150‑4 (MQ=255)
              aCCCAACTAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTTTTTACAGCATGatgttt                                                                                                                                     <  1:166722/150‑5 (MQ=255)
              aCCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTTTTTACAGCATGatgttt                                                                                                                                     <  1:166807/150‑5 (MQ=255)
                ccAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATaaa                                                                                                                                   >  2:1879915/1‑150 (MQ=255)
                ccAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATaaa                                                                                                                                   <  1:1306768/150‑1 (MQ=255)
                ccAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATaaa                                                                                                                                   >  1:807672/1‑150 (MQ=255)
                ccAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATaaa                                                                                                                                   <  1:1307943/150‑1 (MQ=255)
                  aaaTAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAAtt                                                                                                                                 >  1:101852/1‑150 (MQ=255)
                  aaaTAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAAtt                                                                                                                                 >  1:440925/1‑150 (MQ=255)
                  aaaTAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAAtt                                                                                                                                 >  1:2337506/1‑150 (MQ=255)
                          gatgatCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGa                                                                                                                         >  2:2271369/1‑150 (MQ=255)
                          gatgatCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGa                                                                                                                         >  1:1657043/1‑150 (MQ=255)
                          gatgatCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATAAAAAGCGa                                                                                                                         >  1:1657031/1‑150 (MQ=255)
                           atgatCCTTAAGTAGAAAAATAATTCATATCTATCCAAATTAGAAAAAATCCCATTATCTCAATTATTTGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGAt                                                                                                                        <  1:2736496/150‑1 (MQ=255)
                           atgatCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGAt                                                                                                                        <  2:902079/150‑1 (MQ=255)
                           atgatCATTAAGTAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATATCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGAt                                                                                                                        <  1:2736486/150‑1 (MQ=255)
                               tCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTaa                                                                                                                    >  2:1769542/1‑150 (MQ=255)
                                  ttAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAaata                                                                                                                 >  2:154850/1‑150 (MQ=255)
                                       tagaAAAATAATTCATATCTATCCACATTAGAAAAAATCCAATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTa                                                                                                            <  1:2305374/149‑1 (MQ=255)
                                       gagaAAAATAATTCATATCTATCCACATTATAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTa                                                                                                            <  1:1925049/150‑1 (MQ=255)
                                       gagaAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTa                                                                                                            >  2:305091/1‑150 (MQ=255)
                                        agaAAAATAATTCATATCTATCTACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTAt                                                                                                           <  2:25244/150‑1 (MQ=255)
                                        agaAAAATAATTCATATCTATCTACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTAt                                                                                                           <  2:251199/150‑1 (MQ=255)
                                        agaAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTAt                                                                                                           <  1:856886/150‑1 (MQ=255)
                                             aataatTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGcaca                                                                                                      >  2:2402554/1‑150 (MQ=255)
                                              ataatTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACaa                                                                                                     <  2:694648/150‑1 (MQ=255)
                                               taatTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACaaa                                                                                                    <  1:574598/150‑1 (MQ=255)
                                               taatTCATATCTATCCACATGAGAAAAAATCCCATTATCTCTATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACaaa                                                                                                    <  1:2360535/150‑1 (MQ=255)
                                                  ttCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGtt                                                                                                 >  2:532471/1‑150 (MQ=255)
                                                   tCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATGATTTATGCACAAAGTTa                                                                                                >  2:1661175/1‑150 (MQ=255)
                                                    cATATCTATCCACATTAGAAAAAATCCCATTATCTCACTTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTAt                                                                                               <  2:2266825/150‑1 (MQ=255)
                                                     atatCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATc                                                                                              >  1:504644/1‑150 (MQ=255)
                                                     atatCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATc                                                                                              >  1:1883180/1‑150 (MQ=255)
                                                       atctatcCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCca                                                                                            <  2:260996/150‑1 (MQ=255)
                                                           atcCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATg                                                                                        >  2:1729663/1‑150 (MQ=255)
                                                            tcCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGa                                                                                       >  1:2302345/1‑150 (MQ=255)
                                                                 aTTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGAtt                                                                                  >  2:2438453/1‑150 (MQ=255)
                                                                 aTTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGAtt                                                                                  >  2:2061658/1‑150 (MQ=255)
                                                                      aaaaaaTCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTgcga                                                                             >  2:1015066/1‑150 (MQ=255)
                                                                       aaaaaTCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTgcgag                                                                            >  2:699707/1‑150 (MQ=255)
                                                                       aaaaaTCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTgcgag                                                                            >  1:212845/1‑150 (MQ=255)
                                                                       aaaaaTCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTgcgag                                                                            >  1:2484357/1‑150 (MQ=255)
                                                                       aaaaaTCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTgcgag                                                                            <  2:28055/150‑1 (MQ=255)
                                                                          aaTCCCATTATCTCAATTTTTTGGGATGTTTTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTgcgagcga                                                                         <  1:1267592/150‑1 (MQ=255)
                                                                            tCCCAGTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGCCAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATcc                                                                      <  2:1094209/150‑1 (MQ=255)
                                                                              ccATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCa                                                                     >  2:1942089/1‑150 (MQ=255)
                                                                                 ttATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCagaa                                                                  <  1:614240/150‑1 (MQ=255)
                                                                                     ctcAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATc                                                              >  1:943371/1‑150 (MQ=255)
                                                                                     ctcAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATc                                                              <  1:1160610/150‑1 (MQ=255)
                                                                                      tcAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCt                                                             >  2:2012084/1‑150 (MQ=255)
                                                                                      tcAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCt                                                             >  1:222953/1‑150 (MQ=255)
                                                                                      tcAATTATTAGGGAAGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGGTAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCt                                                             >  2:2078545/1‑150 (MQ=255)
                                                                                        aaTTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTAc                                                           <  1:282605/150‑1 (MQ=255)
                                                                                        aaTTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTAc                                                           <  1:282606/150‑1 (MQ=255)
                                                                                           aattaGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAAATAAAACCGATCTAAATAAAATATGCACAAAGTTAACCACATGACGATTTGCGAGCGAACCAGAAGAACTACaaa                                                        <  2:1595133/149‑1 (MQ=255)
                                                                                            attaGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACaaaa                                                       <  2:951587/150‑1 (MQ=255)
                                                                                               aGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGAt                                                    >  1:1816733/1‑150 (MQ=255)
                                                                                                   atggatTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTc                                                >  2:1006457/1‑150 (MQ=255)
                                                                                                    tggatTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCa                                               >  2:1211722/1‑150 (MQ=255)
                                                                                                       atttatttTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGa                                            <  1:2402554/150‑1 (MQ=255)
                                                                                                            tttttAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGc                                       >  1:1255731/1‑150 (MQ=255)
                                                                                                            tttttAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGc                                       <  1:1717950/150‑1 (MQ=255)
                                                                                                            tttttAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGc                                       >  1:2602568/1‑150 (MQ=255)
                                                                                                              tttAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCgg                                     <  1:2074039/150‑1 (MQ=255)
                                                                                                              tttAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCgg                                     <  1:890462/150‑1 (MQ=255)
                                                                                                               ttAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGt                                    <  1:40781/150‑1 (MQ=255)
                                                                                                                  aCTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGaa                                 >  1:1190939/1‑150 (MQ=255)
                                                                                                                    tGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGaaaa                               >  2:240574/1‑150 (MQ=255)
                                                                                                                     gCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGaaaaa                              <  1:1729663/150‑1 (MQ=255)
                                                                                                                      cATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAAc                             >  2:272418/1‑150 (MQ=255)
                                                                                                                        tGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTAc                                                           >  1:168196/1‑118 (MQ=255)
                                                                                                                        tGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTAc                                                           <  2:168196/118‑1 (MQ=255)
                                                                                                                          aaaaaCAAAGACAAACATCCTGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCAt                         <  1:875512/150‑1 (MQ=255)
                                                                                                                          aaaaaCAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCAt                         >  2:373893/1‑150 (MQ=255)
                                                                                                                          aaaaaCAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCAt                         >  1:82262/1‑150 (MQ=255)
                                                                                                                          aaaaaCAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCAt                         >  1:1329941/1‑150 (MQ=255)
                                                                                                                              aCAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGttt                     >  2:434216/1‑150 (MQ=255)
                                                                                                                              aCAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGttt                     >  1:352644/1‑150 (MQ=255)
                                                                                                                               cAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGtttt                    >  2:1154781/1‑150 (MQ=255)
                                                                                                                                 aaGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCa                  <  1:1410343/150‑1 (MQ=255)
                                                                                                                                      aaaCATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCctg             >  1:2438355/1‑150 (MQ=255)
                                                                                                                                      aaaCATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCctg             >  2:1428107/1‑150 (MQ=255)
                                                                                                                                      aaaCATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATATTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCctg             >  2:2175588/1‑150 (MQ=255)
                                                                                                                                        ttatcatGCTTTTTACAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCctgtc           >  1:2551791/3‑150 (MQ=255)
                                                                                                                                            catGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCtgtg       <  1:1560232/150‑1 (MQ=255)
                                                                                                                                            catGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCtgtg       <  2:587237/150‑1 (MQ=255)
                                                                                                                                            catGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCtgtg       <  1:955377/150‑1 (MQ=255)
                                                                                                                                             atGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCTGTgg      >  1:860506/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCTGTGGCATC  >  NC_000913/4046636‑4046930

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