Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,400,313 A→C 100% D9A (GAT→GCT)  hfq → global sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,400,3130AC97.2% 363.9 / ‑0.4 102D9A (GAT→GCT) hfqglobal sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base C (43/56);  total (45/57)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 102 total aligned reads displayed.

AGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATC  >  NC_000913/4400216‑4400408
                                                                                                 |                                                                                               
agaCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTc                                                                                               <  2:492042/100‑1 (MQ=255)
agaCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTc                                                                                               <  2:1331457/100‑1 (MQ=255)
 naCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTc                                                                                               >  2:1172707/2‑99 (MQ=255)
 gaCGTATCGTGCGCAATTTTTTTTTTTTTGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTcc                                                                                              >  2:1708360/1‑100 (MQ=255)
 gaCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTcc                                                                                              >  1:1983358/1‑100 (MQ=255)
 gaCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTcc                                                                                              <  1:1844121/100‑1 (MQ=255)
  aCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCg                                                                                             <  2:959572/100‑1 (MQ=255)
   cGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGt                                                                                            >  1:1029603/1‑100 (MQ=255)
     tATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTc                                                                                          <  2:76550/100‑1 (MQ=255)
      aTCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTcc                                                                                         >  1:122189/1‑100 (MQ=255)
      aTCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTcc                                                                                         <  2:1658445/100‑1 (MQ=255)
         gTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTg                                                                                       <  2:1241385/99‑1 (MQ=255)
         gTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGa                                                                                      <  1:69788/100‑1 (MQ=255)
          tGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCGATCTTTACAATATCCGTTCTTGaa                                                                                     >  1:665377/1‑100 (MQ=255)
           gcgcAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAAc                                                                                    >  2:715999/1‑100 (MQ=255)
           gcgcAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAAc                                                                                    <  2:1048619/100‑1 (MQ=255)
           gcgcAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAAc                                                                                    <  2:367949/100‑1 (MQ=255)
            cgcAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACg                                                                                   >  2:1661881/1‑100 (MQ=255)
            cgcAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACg                                                                                   <  1:573822/100‑1 (MQ=255)
              cAATTTTTTTATTATCTTAAGTTTCAAAGTACAAATAAGCATATAAGTAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCa                                                                                 >  2:613601/1‑100 (MQ=255)
              cAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATATTTACAAGCTCAGTTCCTGAACGCa                                                                                 >  2:1958269/1‑100 (MQ=255)
                aTTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACt                                                                               <  1:235927/100‑1 (MQ=255)
                 ttttttCAGAATCGAAAAGTTCAAAGAACAAATAAGCATACAAGGACCAGAGAGAATGGCCTAGAAGCAATAATTACAAGCTCAGTTCCTCAACGCAc                                                                                >  2:194915/1‑98 (MQ=255)
                  tttttCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGc                                                                             <  1:387642/100‑1 (MQ=255)
                       cAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCTATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCgg                                                                        <  2:1110015/100‑1 (MQ=255)
                        aGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCggg                                                                       >  1:1135553/1‑100 (MQ=255)
                        aGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCggg                                                                       >  1:379341/1‑100 (MQ=255)
                           aTCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCCGACGCAGTGCGttcaggaac                                                                   <  1:1988381/100‑9 (MQ=37)
                             cGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACg                                                                   <  1:107783/99‑1 (MQ=255)
                              gAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACgtg                                                                 >  2:167874/1‑100 (MQ=255)
                              gAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACgtg                                                                 >  2:667659/1‑100 (MQ=255)
                              gAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACg                                                                   >  2:1832122/1‑98 (MQ=255)
                               aaaGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACgtgt                                                                <  1:679674/100‑1 (MQ=255)
                                aaGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGtt                                                               <  2:1668220/100‑1 (MQ=255)
                                 aGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTc                                                              <  1:1319818/100‑1 (MQ=255)
                                  ggTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTcc                                                             <  1:1688926/100‑1 (MQ=255)
                                  ggTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTcc                                                             >  2:2032159/1‑100 (MQ=255)
                                  ggTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTcc                                                             <  2:1320481/100‑1 (MQ=255)
                                     tCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGt                                                          >  2:715597/1‑100 (MQ=255)
                                     tCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGt                                                          >  1:1816861/1‑100 (MQ=255)
                                      cAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGtt                                                         >  2:1596081/1‑100 (MQ=255)
                                        aaGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTc                                                       <  1:1382702/100‑1 (MQ=255)
                                        aaGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTc                                                       <  2:1054751/100‑1 (MQ=255)
                                        aaGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTc                                                       <  2:373500/100‑1 (MQ=255)
                                           tACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCtat                                                    <  1:481603/100‑1 (MQ=255)
                                           tACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCtat                                                    <  2:172050/100‑1 (MQ=255)
                                            aCAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCtatt                                                   <  1:446202/100‑1 (MQ=255)
                                             cAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCtattt                                                  >  2:1659162/1‑100 (MQ=255)
                                             cAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCtattt                                                  <  2:265775/100‑1 (MQ=255)
                                             cAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCtattt                                                  <  2:1567799/100‑1 (MQ=255)
                                             cAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCtattt                                                  >  1:1930307/1‑100 (MQ=255)
                                              aaaTAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCtattta                                                 <  1:1920075/100‑1 (MQ=255)
                                               aaTAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCtatttat                                                >  1:1003447/1‑100 (MQ=255)
                                                aTAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCtatttatt                                               <  1:1661893/100‑1 (MQ=255)
                                                 tAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTatttattt                                              >  2:537646/1‑100 (MQ=255)
                                                  aaGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTg                                             >  2:250729/1‑100 (MQ=255)
                                                  aaGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTg                                             <  2:1481429/100‑1 (MQ=255)
                                                  aaGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTg                                             <  1:383475/100‑1 (MQ=255)
                                                   aGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTg                                             >  1:905692/1‑99 (MQ=255)
                                                    gCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGt                                           >  1:412141/1‑100 (MQ=255)
                                                     cATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTg                                          <  1:2061705/100‑1 (MQ=255)
                                                     cATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTg                                          <  1:2061838/100‑1 (MQ=255)
                                                       tataAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGaa                                        >  1:360054/1‑100 (MQ=255)
                                                       tataAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGaa                                        >  1:402900/1‑100 (MQ=255)
                                                        ataAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAAt                                       >  2:578963/1‑100 (MQ=255)
                                                         taAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATg                                      <  2:179251/100‑1 (MQ=255)
                                                           aGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGt                                    <  2:122187/100‑1 (MQ=255)
                                                             gAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTAt                                  <  2:412138/100‑1 (MQ=255)
                                                             gAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTAt                                  >  2:1253502/1‑100 (MQ=255)
                                                              aaaaGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTAtt                                 <  1:796663/100‑1 (MQ=255)
                                                               aaaGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTa                                <  2:333985/100‑1 (MQ=255)
                                                               aaaGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTa                                <  2:955982/100‑1 (MQ=255)
                                                               aaaGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTa                                <  2:1472864/100‑1 (MQ=255)
                                                                 agagagAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAg                              >  2:614179/1‑100 (MQ=255)
                                                                  gagagaATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGc                             >  1:1773351/1‑100 (MQ=255)
                                                                  gagagaATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGc                             >  2:306765/1‑100 (MQ=255)
                                                                     agaATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTg                           <  1:716002/99‑1 (MQ=255)
                                                                      gaATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCa                         <  1:848995/100‑1 (MQ=255)
                                                                      gaATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCa                         >  2:250630/1‑100 (MQ=255)
                                                                           gCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGc                    >  1:530853/1‑100 (MQ=255)
                                                                           gCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGc                    >  1:569089/1‑100 (MQ=255)
                                                                            cTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCa                   <  1:1703351/100‑1 (MQ=255)
                                                                            cTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCa                   <  1:1585375/100‑1 (MQ=255)
                                                                            cTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCa                   <  2:1135548/100‑1 (MQ=255)
                                                                            cTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCa                   <  2:1019308/100‑1 (MQ=255)
                                                                             tAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCa                   >  2:635767/1‑99 (MQ=255)
                                                                              aaGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCaaa                 <  1:1675028/100‑1 (MQ=255)
                                                                              aaGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCaaa                 <  2:1553167/100‑1 (MQ=255)
                                                                                ggggCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATc               >  1:1897270/1‑100 (MQ=255)
                                                                                ggggCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTACTTATTTGGTGAATGGTATTAAGCTGCAAGGGCaact                >  2:815290/1‑97 (MQ=255)
                                                                                 gggCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATca              >  2:273805/1‑99 (MQ=255)
                                                                                 gggCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATca              >  1:311340/1‑99 (MQ=255)
                                                                                 gggCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATca              <  2:1029600/100‑2 (MQ=255)
                                                                                  ggCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATcaa             >  1:1881209/1‑98 (MQ=255)
                                                                                   gCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAg            <  1:1247600/100‑1 (MQ=255)
                                                                                    cAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGt           >  1:848277/1‑100 (MQ=255)
                                                                                     aaTCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTc          <  2:760050/100‑1 (MQ=255)
                                                                                       tCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTTTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCt         <  1:114329/99‑1 (MQ=255)
                                                                                           taaaaGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCtttt      <  2:168585/95‑1 (MQ=255)
                                                                                             cAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATc  >  2:1979006/1‑100 (MQ=255)
                                                                                                 |                                                                                               
AGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATC  >  NC_000913/4400216‑4400408

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