Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,483,912 T→A 5.9% V53D (GTC→GAC)  evgA → DNA‑binding transcriptional activator EvgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,483,9120TA 5.9% 316.7 / 7.8 101V53D (GTC→GAC) evgADNA‑binding transcriptional activator EvgA
Reads supporting (aligned to +/- strand):  ref base T (40/55);  new base A (2/4);  total (42/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 104 total aligned reads displayed.

AATAATTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTAATGGTTTCGTGAGTAA  >  NC_000913/2483763‑2484059
                                                                                                                                                     |                                                                                                                                                   
aataatTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGt                                                                                                                                                     >  1:752381/1‑150 (MQ=255)
 ataatTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTc                                                                                                                                                    <  2:716039/150‑1 (MQ=255)
   aatTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAGACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGa                                                                                                                                                  >  2:1657378/1‑150 (MQ=255)
   aatTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGa                                                                                                                                                  <  1:1874646/150‑1 (MQ=255)
   aatTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGa                                                                                                                                                  <  2:1009989/150‑1 (MQ=255)
    attattGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGat                                                                                                                                                 >  1:1508063/1‑150 (MQ=255)
    attattGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGGGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGat                                                                                                                                                 <  2:676206/150‑1 (MQ=255)
     ttattGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTCAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACAATTAAGCCTGATATCGTCATCATTGATGTCGaaa                                                                                                                                                >  2:1003749/1‑148 (MQ=255)
          gatgaCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATcccc                                                                                                                                           >  2:730149/1‑150 (MQ=255)
          gatgaCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcagg                                                                                                                                           <  1:272517/150‑4 (MQ=255)
          gatgaCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcagg                                                                                                                                           <  1:1830608/150‑4 (MQ=255)
          gatgaCCATCCTCATGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcagg                                                                                                                                           <  1:1370123/150‑4 (MQ=255)
           atgaCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCg                                                                                                                                          <  1:1640898/150‑1 (MQ=255)
           atgaCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcaggc                                                                                                                                          <  1:1517967/150‑5 (MQ=255)
            tgaCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCgg                                                                                                                                         <  2:1508063/150‑1 (MQ=255)
            tgaCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCgg                                                                                                                                         <  2:1283046/150‑1 (MQ=255)
            tgaCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCgg                                                                                                                                         >  2:1589633/1‑150 (MQ=255)
             gaCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGa                                                                                                                                        <  2:1413944/150‑1 (MQ=255)
              aCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACGCTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAg                                                                                                                                       >  1:1957342/1‑150 (MQ=255)
              aCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAg                                                                                                                                       <  1:882843/150‑1 (MQ=255)
              aCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGAAATCGTCATCATTGATGTCGATATCCCCGGAg                                                                                                                                       >  1:549279/1‑150 (MQ=255)
                 aTCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTa                                                                                                                                    <  1:1408516/150‑1 (MQ=255)
                   ccTCTTGCTATCGCAGCAATTCGTAATTTATTGATCGAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATTTCGATATCCCCGGAGTTAAc                                                                                                                                  >  2:1509365/1‑150 (MQ=255)
                   ccTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAAc                                                                                                                                  <  2:507948/150‑1 (MQ=255)
                    ctctTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACg                                                                                                                                 <  1:1848365/150‑1 (MQ=255)
                    ctctTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACg                                                                                                                                 <  2:612155/150‑1 (MQ=255)
                     tctTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACgg                                                                                                                                <  1:1207316/150‑1 (MQ=255)
                     tctTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACgg                                                                                                                                <  1:576279/150‑1 (MQ=255)
                        tGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTAt                                                                                                                             >  1:1076260/1‑150 (MQ=255)
                          cTATCGCAGCAATTCGTAATTTATTGATCGAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATcc                                                                                                                           <  1:1509365/150‑1 (MQ=255)
                             tCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAgg                                                                                                                        <  1:1849401/150‑1 (MQ=255)
                             tCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAgg                                                                                                                        <  1:708478/150‑1 (MQ=255)
                               gcagcaATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGgtg                                                                                                                      <  2:1220255/150‑1 (MQ=255)
                               gcagcaATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGgtg                                                                                                                      <  2:423145/150‑1 (MQ=255)
                               gcagcaATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGgtg                                                                                                                      <  2:1224309/150‑1 (MQ=255)
                                 agcaATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGAAATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGtt                                                                                                                    <  2:549279/150‑1 (MQ=255)
                                  gcaATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGt                                                                                                                                                     >  1:1007776/1‑116 (MQ=255)
                                  gcaATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTa                                                                                                                   >  1:1637804/1‑150 (MQ=255)
                                    aaTTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGa                                                                                                                 <  2:1285389/150‑1 (MQ=255)
                                    aaTTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGa                                                                                                                 >  2:1484020/1‑150 (MQ=255)
                                    aaTTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGa                                                                                                                 <  2:1948540/150‑1 (MQ=255)
                                         gTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGc                                                                                                            <  2:1539934/150‑1 (MQ=255)
                                           aaTTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTg                                                                                                          >  1:1629303/1‑150 (MQ=255)
                                             tttATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAg                                                                                                        <  2:1281236/150‑1 (MQ=255)
                                             tttATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAg                                                                                                        <  1:1758224/150‑1 (MQ=255)
                                               tATTGATCAAAAAGGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGa                                                                                                      <  1:315173/150‑1 (MQ=255)
                                                   gATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAgcg                                                                                                  <  2:630169/150‑1 (MQ=255)
                                                     tCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGcc                                                                                                >  2:238947/1‑150 (MQ=255)
                                                       aaaaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCaa                                                                                              >  1:1288469/1‑150 (MQ=255)
                                                       aaaaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCa                                                                                               <  2:453300/149‑1 (MQ=255)
                                                       aaaaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCa                                                                                               >  1:453300/1‑149 (MQ=255)
                                                         aaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGATCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAata                                                                                            >  2:1300912/1‑150 (MQ=255)
                                                           aCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAAtata                                                                                          <  1:1228708/150‑1 (MQ=255)
                                                            cGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAg                                                                                         <  2:1629303/150‑1 (MQ=255)
                                                            cGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAg                                                                                         >  1:1558651/1‑150 (MQ=255)
                                                             catatTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGc                                                                                        <  2:1964944/149‑1 (MQ=255)
                                                                   gAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAatt                                                                                  >  1:11021/1‑150 (MQ=255)
                                                                        cTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattattat                                                                             <  2:887914/150‑1 (MQ=255)
                                                                            gCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATc                                                                         >  2:1960971/1‑150 (MQ=255)
                                                                            gCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATc                                                                         >  1:988894/1‑150 (MQ=255)
                                                                             cAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCg                                                                        <  2:1637804/150‑1 (MQ=255)
                                                                             cAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCg                                                                        <  1:1418086/150‑1 (MQ=255)
                                                                                  ttGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTcc                                                                   <  1:1960971/150‑1 (MQ=255)
                                                                                   tGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCACCAATATAGCGGAATTATTATTATCGTCTCCg                                                                  >  1:512518/1‑150 (MQ=255)
                                                                                      cTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGGCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGCAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTa                                                               >  2:1087371/1‑150 (MQ=255)
                                                                                       tGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTaa                                                              >  1:1966191/1‑150 (MQ=255)
                                                                                           ggCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAAt                                                          >  1:1270219/1‑150 (MQ=255)
                                                                                            gCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATg                                                         <  1:238947/150‑1 (MQ=255)
                                                                                             cGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATgg                                                        >  1:1110913/1‑149 (MQ=255)
                                                                                                  gTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCAtt                                                   <  1:1300912/150‑1 (MQ=255)
                                                                                                   tGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCAttt                                                  >  2:1263361/1‑150 (MQ=255)
                                                                                                    gCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAACGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCAtttt                                                 >  2:137886/1‑150 (MQ=255)
                                                                                                       gTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTAc                                              <  2:116614/150‑1 (MQ=255)
                                                                                                       gTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCCATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTAc                                              <  2:589219/150‑1 (MQ=255)
                                                                                                        ttCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGGTCATTTTTagc                                             <  1:1657378/150‑3 (MQ=255)
                                                                                                          cAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAATATGCCCATTTTTACggg                                           >  2:1896086/1‑150 (MQ=255)
                                                                                                          cAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGGTCATTTTTagcgg                                           <  2:175011/150‑5 (MQ=255)
                                                                                                          cAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACggg                                           <  2:1076260/150‑1 (MQ=255)
                                                                                                           aGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGa                                          >  2:1706516/1‑150 (MQ=255)
                                                                                                                 tGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACAtt                                    >  2:1489030/1‑150 (MQ=255)
                                                                                                                  ggAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATtg                                   >  1:311406/1‑150 (MQ=255)
                                                                                                                  ggAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCACCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATtg                                   <  2:512518/150‑1 (MQ=255)
                                                                                                                   gAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATtgt                                  >  2:1243900/1‑150 (MQ=255)
                                                                                                                     aaCACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGc                                >  1:377800/1‑150 (MQ=255)
                                                                                                                      acacTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCt                               <  2:329365/150‑1 (MQ=255)
                                                                                                                           tAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGc                          >  1:563049/1‑150 (MQ=255)
                                                                                                                              tccTGATATCGTCATCATTTATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGTAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTgg                       <  2:861567/149‑1 (MQ=255)
                                                                                                                               ggggATATCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGgc                      >  2:508867/4‑150 (MQ=255)
                                                                                                                                gggATATCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGgcg                     >  1:742540/3‑150 (MQ=255)
                                                                                                                                cTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGgcg                     >  1:1760404/1‑150 (MQ=255)
                                                                                                                                    tatCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTaa                 >  1:1693470/1‑150 (MQ=255)
                                                                                                                                     atCGGCATCATTGATGTCTATATCACCGGAGTTAACGGTATCCAGGTGTTATAAACGCTGAGTAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTTCTTATGCTGGCGCTAAt                <  2:829883/150‑1 (MQ=255)
                                                                                                                                      tCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTAATg               <  1:730149/150‑1 (MQ=255)
                                                                                                                                      tCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTAATg               <  1:1003749/150‑1 (MQ=255)
                                                                                                                                       cGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTAATgg              >  2:550666/1‑150 (MQ=255)
                                                                                                                                        gTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTAATGGt             >  2:416646/1‑150 (MQ=255)
                                                                                                                                           atcatTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTAATGGTTTc          <  2:1947507/150‑1 (MQ=255)
                                                                                                                                             tatTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGTAAGCGCAAATATAGCGTAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGGTGATGCTGGCGCTAATGGTTTCGt        <  2:223610/149‑1 (MQ=255)
                                                                                                                                             catTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTAATGGTTTCGt        <  1:726386/150‑1 (MQ=255)
                                                                                                                                                   tGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTAATGGTTTCGTGAGTaa  >  1:1183110/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AATAATTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTAATGGTTTCGTGAGTAA  >  NC_000913/2483763‑2484059

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