Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,046,782 G→T 6.2% intergenic (‑180/‑184) 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,7820GT 6.2% 587.4 / 20.5 162intergenic (‑180/‑184)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 G (72/80);  new base T (2/8);  total (74/88)
Fisher's exact test for biased strand distribution p-value = 1.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

Only 100 of 162 total aligned reads displayed.

ATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCTGTGGCATC  >  NC_000913/4046634‑4046930
                                                                                                                                                    |                                                                                                                                                    
aTTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGt                                                                                                                                                     >  2:1659606/1‑150 (MQ=255)
   ttAGCATAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTaaa                                                                                                                                                  <  1:767275/150‑1 (MQ=255)
   ttAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTaaa                                                                                                                                                  <  2:2213095/150‑1 (MQ=255)
   ttAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTaaa                                                                                                                                                  <  2:2213099/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:2303666/150‑1 (MQ=255)
        aGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCa                                                                                                                                             <  2:2714904/150‑1 (MQ=255)
         gAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCAt                                                                                                                                            <  2:2313884/150‑1 (MQ=255)
            aaCGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGat                                                                                                                                         >  1:137811/1‑150 (MQ=255)
              cGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTTTTTACAGCATGatgt                                                                                                                                       <  2:2735848/150‑3 (MQ=255)
               gACCCAAATACCGGATGATCCTTAAGGAGAAAAATAATTCATATCTATACACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTTTTTTTTACTGCATGAAATACAATGACAAACATCATGCTTTTTACAGCATGatgtt                                                                                                                                      <  2:2418185/150‑4 (MQ=255)
                aCCCAACTAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTTTTTACAGCATGatgttt                                                                                                                                     <  1:166722/150‑5 (MQ=255)
                aCCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTTTTTACAGCATGatgttt                                                                                                                                     <  1:166807/150‑5 (MQ=255)
                 cccAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCATTTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCCGTAAAAAGCATGAtaataa                                                                                                                                    <  2:2030842/150‑1 (MQ=255)
                  ccAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATaaa                                                                                                                                   >  1:807672/1‑150 (MQ=255)
                  ccAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATaaa                                                                                                                                   <  1:1307943/150‑1 (MQ=255)
                  ccAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATaaa                                                                                                                                   <  1:1306768/150‑1 (MQ=255)
                  ccAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATaaa                                                                                                                                   >  1:497303/1‑150 (MQ=255)
                    aaaTAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAAtt                                                                                                                                 >  1:101852/1‑150 (MQ=255)
                            gatgatCCTTAATGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGa                                                                                                                         <  2:2654628/150‑1 (MQ=255)
                            gatgatCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGTGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGa                                                                                                                         <  2:1903996/150‑1 (MQ=255)
                            gatgatCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATAAAAAGCGa                                                                                                                         >  1:1657031/1‑150 (MQ=255)
                             atgatCCTTAAGGAGAAAAATAATTCTTATCTATCCACATTATAAAAAATCCCATTATCTCAATTATTAGGGTTGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGAt                                                                                                                        <  1:2736512/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)
                                      aaGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAaataat                                                                                                               <  1:403399/150‑1 (MQ=255)
                                         tagaAAAATAATTCATATCTATCCACATTAGAAAAAATCCAATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTa                                                                                                            <  1:2305374/149‑1 (MQ=255)
                                         gagaAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTa                                                                                                            >  2:305091/1‑150 (MQ=255)
                                         gagaAAAATAATTCATATCTATCCACATGAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTa                                                                                                            <  2:1639665/150‑1 (MQ=255)
                                          agaAAAATAATTCATATCTATCTACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTAt                                                                                                           <  2:251199/150‑1 (MQ=255)
                                          agaAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTAt                                                                                                           <  1:856886/150‑1 (MQ=255)
                                          agaAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTAt                                                                                                           <  1:2183473/150‑1 (MQ=255)
                                          agaAAAATAATTCATATATATCTACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACATAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTAt                                                                                                           <  2:25241/150‑1 (MQ=255)
                                                ataatTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACaa                                                                                                     <  2:694648/150‑1 (MQ=255)
                                                 taatTCATATCTATCCACATGAGAAAAAATCCCATTATCTCTATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACaaa                                                                                                    <  1:2360535/150‑1 (MQ=255)
                                                    ttCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGtt                                                                                                 >  2:532471/1‑150 (MQ=255)
                                                      cATATCTATCCACATTAGAAAAAATCCCATTATCTCACTTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTAt                                                                                               <  2:2266825/150‑1 (MQ=255)
                                                       atatCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATc                                                                                              >  1:1883180/1‑150 (MQ=255)
                                                         atctatcCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCca                                                                                            <  2:260996/150‑1 (MQ=255)
                                                          tctatcCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCcac                                                                                           <  2:581488/150‑1 (MQ=255)
                                                              tcCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGa                                                                                       >  1:2302345/1‑150 (MQ=255)
                                                                   aTTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGAtt                                                                                  >  2:2438453/1‑150 (MQ=255)
                                                                      aGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTgc                                                                               <  2:1680974/150‑1 (MQ=255)
                                                                      aGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTgc                                                                               <  1:2354561/150‑1 (MQ=255)
                                                                        aaaaaaTCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTgcga                                                                             >  2:1015066/1‑150 (MQ=255)
                                                                         aaaaaTCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTgcgag                                                                            >  1:1660296/1‑150 (MQ=255)
                                                                         aaaaaTCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTgcgag                                                                            >  2:699707/1‑150 (MQ=255)
                                                                         aaaaaTCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTgcgag                                                                            >  1:2484357/1‑150 (MQ=255)
                                                                          aaaaTCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTgcgagc                                                                           >  2:2197581/1‑150 (MQ=255)
                                                                            aaTCCCATTATCTCAATTTTTTGGGATGTTTTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTgcgagcga                                                                         <  1:1267592/150‑1 (MQ=255)
                                                                              tCCCAGTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGCCAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATcc                                                                      <  2:1094209/150‑1 (MQ=255)
                                                                                 cATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCag                                                                    >  2:1871702/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                                                             >  1:1981352/1‑150 (MQ=255)
                                                                                        tcAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCt                                                             >  2:2012084/1‑150 (MQ=255)
                                                                                        tcAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCt                                                             >  1:222953/1‑150 (MQ=255)
                                                                                        tcAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCt                                                             <  1:1803606/150‑1 (MQ=255)
                                                                                        tcAATTATTAGGGAAGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGGTAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCt                                                             >  2:2078545/1‑150 (MQ=255)
                                                                                          aaTTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTAc                                                           <  1:282605/150‑1 (MQ=255)
                                                                                          aaTTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTAc                                                           <  2:2268530/150‑1 (MQ=255)
                                                                                              attaGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACaaaa                                                       <  2:951587/150‑1 (MQ=255)
                                                                                               ttaGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAg                                                      >  1:2462997/1‑150 (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:2694993/1‑150 (MQ=255)
                                                                                                         atttatttTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGa                                            <  1:2402554/150‑1 (MQ=255)
                                                                                                         atttatatTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGGGAGCGATGCAGAAGATCTACAAAAGATTTTCACGa                                            <  2:2083427/150‑1 (MQ=255)
                                                                                                              tttttAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGc                                       >  1:1255731/1‑150 (MQ=255)
                                                                                                              tttttAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGc                                       <  1:1717950/150‑1 (MQ=255)
                                                                                                                tttAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCgg                                     <  1:2583680/150‑1 (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)
                                                                                                                     cTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGaaa                                >  1:587237/1‑150 (MQ=255)
                                                                                                                      tGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGaaaa                               >  2:240574/1‑150 (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                         >  1:31526/1‑150 (MQ=255)
                                                                                                                            aaaaaCAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCAt                         >  1:82262/1‑150 (MQ=255)
                                                                                                                            aaaaaCAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCAt                         >  2:373893/1‑150 (MQ=255)
                                                                                                                            aaaaaCAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCAt                         >  1:1329941/1‑150 (MQ=255)
                                                                                                                                aCAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGttt                     >  2:434216/1‑150 (MQ=255)
                                                                                                                                 cAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGtttt                    >  2:1154781/1‑150 (MQ=255)
                                                                                                                                   aaGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCa                  <  1:1410343/150‑1 (MQ=255)
                                                                                                                                       cAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCct              <  1:532471/150‑1 (MQ=255)
                                                                                                                                        aaaCATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCctg             >  2:1428107/1‑150 (MQ=255)
                                                                                                                                          ttatcatGCTTTTTACAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGAGTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTGTCATCctgtc           >  2:1695444/3‑150 (MQ=255)
                                                                                                                                              catGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCtgtg       <  1:955377/150‑1 (MQ=255)
                                                                                                                                              catGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCtgtg       <  1:1560232/150‑1 (MQ=255)
                                                                                                                                               atGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCTGTgg      >  1:860506/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCTGTGGCATC  >  NC_000913/4046634‑4046930

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