Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,950,467 G→T 42.1% intergenic (+47/‑40) ilvL → / → ilvX ilvG operon leader peptide/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,950,4670GT42.1% 74.9 / 127.7 109intergenic (+47/‑40)ilvL/ilvXilvG operon leader peptide/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (27/36);  new base T (17/29);  total (44/65)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 111 total aligned reads displayed.

GGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATG  >  NC_000913/3950369‑3950562
                                                                                                   |                                                                                               
ggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCggg‑gt                                                                                                >  1:1862817/1‑99 (MQ=255)
    attattATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTtttt                                                                                            <  1:1286279/100‑1 (MQ=255)
      tattatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTtttttt                                                                                          >  2:2329400/1‑100 (MQ=255)
      tattatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCC‑‑G‑GGtttttttt                                                                                        <  1:1334372/100‑1 (MQ=255)
       attatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑Gtttttttt                                                                                         <  2:1325464/100‑1 (MQ=255)
         tatCCCGCCGTTCGGGGCTACCCTTCGACGAGGAAAGGCTTAGAGATGAAGCCGTTACTTAATAAGACCCCCTCACCGAAATGTCCGGG‑TTTTTttgttg                                                                                       <  2:2116350/100‑1 (MQ=255)
         tatCCCACCGTGCGGGGTTGCACTTGGACGAGGAAAGGCTTAGGGATCAAGCCTTAACTAACTAATACCCCCGCATTGAAATGTCCGTG‑GTGTTTTTTTg                                                                                       <  2:1154802/100‑1 (MQ=255)
         tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTg                                                                                       <  1:1264336/100‑1 (MQ=255)
         tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTg                                                                                       <  1:1308205/100‑1 (MQ=255)
         tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTg                                                                                       <  2:2289185/100‑1 (MQ=255)
         tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GGTTTTTTTTg                                                                                       <  2:849301/100‑1 (MQ=255)
          atCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGa                                                                                      >  1:375311/1‑100 (MQ=255)
          atCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGa                                                                                      >  1:1060948/1‑100 (MQ=255)
            cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGAcc                                                                                    <  1:1696245/100‑1 (MQ=255)
            cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGAcc                                                                                    <  1:203734/100‑1 (MQ=255)
            cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGAc                                                                                     >  1:767340/1‑99 (MQ=255)
            cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCC‑GG‑GTTTTTTTTTGACCt                                                                                   >  1:348861/1‑100 (MQ=255)
             ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCt                                                                                   <  1:1002564/100‑1 (MQ=255)
             ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCt                                                                                   >  1:384040/1‑100 (MQ=255)
             ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCt                                                                                   >  1:2377644/1‑100 (MQ=255)
              cACGGTGGGGGGTTTCAATTTGACGAGGAAGGGGTTAGACGTGAAGCCGTAAAGAACTTAGGCCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCtt                                                                                  <  2:943571/100‑1 (MQ=255)
              cACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCt                                                                                   >  1:1583606/1‑99 (MQ=255)
              cACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCCCACCGAAAGGTCC‑GG‑GTTTTTTTTTGACCTTa                                                                                 >  1:607113/1‑100 (MQ=255)
                  gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTaaaa                                                                              <  1:1206867/100‑1 (MQ=255)
                     cGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACa                                                                           >  1:1478020/1‑100 (MQ=255)
                     cGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACa                                                                           <  1:910097/100‑1 (MQ=255)
                      ggggTTACACTTAGCCGAGGGAAGGGTTAGAGATCAAGCCTAAACAAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTTGACCTTAAAAACa                                                                           <  1:1762429/100‑1 (MQ=255)
                      ggggCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACAt                                                                          >  1:2121898/1‑100 (MQ=255)
                      ggggATGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCGTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACa                                                                           <  1:494849/99‑1 (MQ=255)
                       gggCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACa                                                                           <  1:2011724/98‑1 (MQ=255)
                       gggCTGCACTCGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATa                                                                         <  2:1574913/100‑1 (MQ=255)
                         gCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTTAACTTaaaaaaaaaaa                                                                       >  1:212847/1‑94 (MQ=255)
                         gCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAAc                                                                       >  2:692860/1‑100 (MQ=255)
                           tGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAAcc                                                                      <  1:829579/99‑1 (MQ=255)
                            gCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCga                                                                    >  1:1196522/1‑100 (MQ=255)
                             cACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCC‑GG‑GGTTTTTTTTGACCTTAAAAACATAACCgagg                                                                  <  1:2327318/100‑1 (MQ=255)
                               cTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCgagga                                                                 >  1:553976/1‑100 (MQ=255)
                                 tGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                               <  2:2385266/100‑1 (MQ=255)
                                  ggACGAGGAAAGGCTCAGAGATCGAGCCTTAACGAACTAAGACACCCGCACCGAAAGGTCC‑GG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                              <  1:1062418/99‑1 (MQ=255)
                                   gACGAGGAAAGGTTTAAAGATCACGCGTTAGGGAAATAAAATGCCCGGACCGAAAGGTCAGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                              <  2:722779/99‑1 (MQ=255)
                                   gACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACGAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTGTGACCTTAAAAACATAACCGAGAAGCa                                                              >  2:344404/1‑99 (MQ=255)
                                       aGGGAAGGCCTAGTGATAAAGCCTGAGTGAAGCAAGCCACACGAAACGACAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGAc                                                           <  2:703334/98‑1 (MQ=255)
                                       aGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACaa                                                         >  1:1462545/1‑100 (MQ=255)
                                       aGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCC‑GG‑GGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                        >  2:150367/1‑100 (MQ=255)
                                       aGGAAAGGCTGAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                              <  1:241802/95‑1 (MQ=255)
                                        ggAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                        >  1:1593422/1‑100 (MQ=255)
                                         gagaGGATTAGAGACCAAGCGTGAGGGGACGAAGGCGCCCGCACAGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                       <  1:815292/97‑1 (MQ=255)
                                         gAAATGCTTAGAGATCAAGGCTTAACGAAATAAGAACCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                       <  1:792804/100‑1 (MQ=255)
                                         gAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCTAGGAGCAGACAATg                                                       >  2:1361065/1‑100 (MQ=255)
                                         gAAAGGCTTAGAGACCAAGCCGTAGCGACCGAAGACCCCCGCACAGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                       <  1:2099518/100‑1 (MQ=255)
                                          aaaGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGa                                                      >  1:2308521/1‑100 (MQ=255)
                                             gggttAGAGATCAAGCCTTAACGAACTAAGTCCCCCGCACCGAAAGGACCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                       <  2:2372746/93‑1 (MQ=255)
                                             ggCTTAGAGATCAATCCTTACCGAACTGAGGCCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAAt                                                    <  2:714037/99‑1 (MQ=255)
                                             ggCTTAAAGATCAAGCCTGAAGGAAGGAAGACCCCCGCACCGAAAGGGCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATa                                                   <  1:466358/100‑1 (MQ=255)
                                               tttAGAGACCAGGCCTGAACGAACTAGGACCCCGGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                        <  1:1366162/92‑1 (MQ=255)
                                               cTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                 >  1:3535/1‑100 (MQ=255)
                                                ttAGAGATCAAGCCTTAACGAACTAAGACCACCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                       <  1:219805/93‑1 (MQ=255)
                                                 tAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GGTTTTTTTTAAAATTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                 <  2:357030/98‑1 (MQ=255)
                                                 tAGAGATCAAGCCGTAACGGACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                 <  1:816537/98‑1 (MQ=255)
                                                   gagaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGca                                             >  2:1624090/1‑100 (MQ=255)
                                                     gaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGcaca                                           >  1:1285731/1‑100 (MQ=255)
                                                     gaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTTACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGcaca                                           >  2:899103/1‑100 (MQ=255)
                                                      aTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaa                                          >  1:1393021/1‑100 (MQ=255)
                                                       tCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTTACTTTAAAAACTTAACCGAGGAGCATACAATGAATAACAGcaca                                           >  2:182201/1‑98 (MQ=255)
                                                       tCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaaa                                         >  2:1759965/1‑100 (MQ=255)
                                                       tCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGGCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGGAAACAATGAATAACAGAACcaa                                         >  1:1917329/1‑100 (MQ=255)
                                                        cAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaaaa                                        >  1:1715419/1‑100 (MQ=255)
                                                        cAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTCCTTTTTTGACCTCAAAACCATAACCGAGCAGCAGACCATGAATAACAGCACaaa                                         >  1:57292/1‑99 (MQ=255)
                                                           gCCTTAACGAATTTATACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAt                                       <  2:1057451/98‑1 (MQ=255)
                                                            ccTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCt                                    >  1:2372932/1‑100 (MQ=255)
                                                             cTTAACGAACTAAGACCCCCGCGCCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTg                                   <  1:1659619/100‑1 (MQ=255)
                                                              ttAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATGACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGt                                  >  1:1171729/1‑100 (MQ=255)
                                                              ttAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGt                                  >  1:1245761/1‑100 (MQ=255)
                                                               tAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGtt                                 >  1:1937079/1‑100 (MQ=255)
                                                                  cGAACTTAGCCCCTCGCACAGAAAGGTCCGGG‑GTATTTTTTTGACCTTAAAAACATAACCGAGGAGCAAACAATGAATAACAGCACAAAATTCTGTTtc                               <  1:2386109/99‑1 (MQ=255)
                                                                   gaacatagtcccccGCACGGGAGGGGCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAACAACAGCACAAAATTCTGTTtc                               <  1:526491/89‑1 (MQ=255)
                                                                     aCTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCaa                           <  1:1140405/100‑1 (MQ=255)
                                                                      cTAAGTCCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCaa                           <  1:293853/99‑1 (MQ=255)
                                                                      cTAAGACCCCCGCACCGAAGGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAg                          <  1:2314322/100‑1 (MQ=255)
                                                                      cTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAg                          >  1:9776/1‑100 (MQ=255)
                                                                       ttggACGCCTGCAACGAAAGGGCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtctc                             <  2:820760/93‑1 (MQ=255)
                                                                            cccccGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGAtt                       <  2:1715399/97‑1 (MQ=255)
                                                                              cccGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGa                  >  1:410671/1‑100 (MQ=255)
                                                                               ccGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGAc                 <  1:913899/100‑1 (MQ=255)
                                                                                cGCAGCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAgg                   <  2:2147825/97‑1 (MQ=255)
                                                                                 gcgccgattggTCCGGG‑TTTTTTTGTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgg               <  1:420760/91‑1 (MQ=255)
                                                                                   agcGAATGGTCCGGT‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgggg             <  2:26527/98‑1 (MQ=255)
                                                                                   aCCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgggg             <  2:1591872/100‑1 (MQ=255)
                                                                                    ccGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGa            >  1:76411/1‑100 (MQ=255)
                                                                                      gacaGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACggg              <  2:1285422/93‑1 (MQ=255)
                                                                                      gAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac          >  1:1709937/1‑100 (MQ=255)
                                                                                      gAAAGGTCCGGG‑GTTTTTTTTTTACCTTAAAAACATAACCGAGGAGAAGACAATAAATAACAGCACAAAATTCTGTTTCTCAAGATTAAGAACGGGaaa           >  2:1315955/1‑99 (MQ=255)
                                                                                      gAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac          <  1:1598773/100‑1 (MQ=255)
                                                                                       aaaGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaact         <  1:1529695/100‑1 (MQ=255)
                                                                                         aGGTCCGGG‑GTTTTTTTTTTTTTTTTTAAACATAACCGAGGACCAGACAATGAATAACAGCACAAAATTCTTTTTCTCAAGATTCACGACGGGGaactaa       >  1:292940/1‑100 (MQ=255)
                                                                                          ggTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCCGTTTCTCAAGATTCAGGACGGGGaactaac      <  2:2380321/100‑1 (MQ=255)
                                                                                          ggTCCGGG‑GTTTTTTTTTGACCTGAAAAACATAACCGAGGAGCAGACAATGAATACCAGCACAAAACTCTGTTTCTCAAGATTCAGGACGGGGaacga        <  1:1620459/98‑3 (MQ=255)
                                                                                            tCCGGG‑TTTTTTTTTTGACCTTAACAACATAACCGAGGAGCAGGCAATGAATAACTGCACATAATTCTGTTTCTCAAGATTCAGGCCGGGGaactaac      <  2:866245/98‑1 (MQ=255)
                                                                                              cGGG‑GTTTTTTTTTGACCTGAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATg  <  1:1238302/100‑1 (MQ=255)
                                                                                               ggg‑gtTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGAAGGGGAACTAACTATg  <  2:384037/99‑1 (MQ=255)
                                                                                                   |                                                                                               
GGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATG  >  NC_000913/3950369‑3950562

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