Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,400,396 G→A 100% E37K (GAG→AAG)  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,3960GA91.1% 256.3 / 11.4 88E37K (GAG→AAG) hfqglobal sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication
Reads supporting (aligned to +/- strand):  ref base G (4/4);  new base A (36/44);  total (40/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTT  >  NC_000913/4400299‑4400491
                                                                                                 |                                                                                               
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)
            agctCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAg                                                                                   <  2:863468/97‑1 (MQ=255)
            agctCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAg                                                                                   >  2:85476/4‑100 (MQ=255)
             gctCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCCGTTTCTATTTATTTGGTGATTGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGt                                                                                  <  2:1134309/98‑1 (MQ=255)
              ctCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGtt                                                                                 <  1:667662/99‑1 (MQ=255)
                 cGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGt                                                                              <  1:273808/100‑1 (MQ=255)
                 cGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGt                                                                              >  1:960815/1‑100 (MQ=255)
                  gTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTg                                                                             >  2:230370/1‑100 (MQ=255)
                    tCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGAt                                                                           >  1:617609/1‑100 (MQ=255)
                     ccTGAACGCACTGCGTCGGGCACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATc                                                                          >  2:115680/1‑100 (MQ=255)
                     ccTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATc                                                                          >  2:1044818/1‑100 (MQ=255)
                      cTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATcc                                                                         <  1:306768/100‑1 (MQ=255)
                        gAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTg                                                                       >  1:1630299/1‑100 (MQ=255)
                        gAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTg                                                                       >  2:225747/1‑100 (MQ=255)
                        gAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTg                                                                       <  2:200932/100‑1 (MQ=255)
                        gAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATAAGTTCGTGATCCTg                                                                       >  1:485068/1‑100 (MQ=255)
                          aCGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGtt                                                                     >  1:1564779/1‑100 (MQ=255)
                           cGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTg                                                                    <  2:1051262/100‑1 (MQ=255)
                           cGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTg                                                                    <  1:1985144/100‑1 (MQ=255)
                            gCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGa                                                                   <  1:1313020/100‑1 (MQ=255)
                             cACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGaa                                                                  >  1:1802055/1‑100 (MQ=255)
                             cACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGaa                                                                  >  2:1647497/1‑100 (MQ=255)
                              aCTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGaaa                                                                 <  1:220423/100‑1 (MQ=255)
                              aCTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGaaa                                                                 <  2:379338/100‑1 (MQ=255)
                               cTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGaaaa                                                                >  1:1069769/1‑100 (MQ=255)
                                     cGGGAACGTGTTCCAGTTTCTATTTATTTGGTGTATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACgg                                                          <  2:168038/100‑1 (MQ=255)
                                     cGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAacac                                                            <  2:1181651/98‑1 (MQ=255)
                                     cGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACgg                                                          <  2:311337/100‑1 (MQ=255)
                                      gggAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGt                                                         <  2:386656/100‑1 (MQ=255)
                                        gAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCa                                                       <  2:402897/100‑1 (MQ=255)
                                         aaCGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAg                                                      <  1:260896/100‑1 (MQ=255)
                                           cGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGcc                                                    <  1:115682/100‑1 (MQ=255)
                                              gtTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGa                                                 >  1:604422/1‑100 (MQ=255)
                                              gtTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGa                                                 >  1:1843637/1‑100 (MQ=255)
                                              gtTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGa                                                 >  2:652014/1‑100 (MQ=255)
                                               ttCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGAt                                                <  1:1958285/100‑1 (MQ=255)
                                                tCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATg                                               >  1:371629/1‑100 (MQ=255)
                                                 ccAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATgg                                              >  1:1228888/1‑100 (MQ=255)
                                                   aGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGtt                                            >  2:1302886/1‑100 (MQ=255)
                                                    gTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGttt                                           >  2:441138/1‑100 (MQ=255)
                                                      ttttatttattGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACa                                        <  1:1596092/97‑1 (MQ=255)
                                                       tCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACa                                        <  2:848274/100‑1 (MQ=255)
                                                       tCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACa                                        <  2:1816848/100‑1 (MQ=255)
                                                       tCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACa                                        <  1:250732/100‑1 (MQ=255)
                                                       tCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACa                                        <  1:1979022/100‑1 (MQ=255)
                                                       tCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACa                                        <  1:1751064/100‑1 (MQ=255)
                                                         tatttattTGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACATGGTCAGCCAGATGGTTTACa                                        <  2:199806/98‑1 (MQ=255)
                                                          atttatttGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCAAGATGGTTTACAAGc                                     >  1:1929096/1‑100 (MQ=255)
                                                              atttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAcgc                                 >  1:2026786/1‑100 (MQ=255)
                                                                 tGGTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGa                               >  2:1035551/1‑99 (MQ=255)
                                                                  ggTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGAtt                             <  1:1445679/100‑1 (MQ=255)
                                                                  ggTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGAtt                             >  2:947255/1‑100 (MQ=255)
                                                                  ggTGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGAtt                             >  2:1489769/1‑100 (MQ=255)
                                                                    tGAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTc                           <  1:1789711/100‑1 (MQ=255)
                                                                     gAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCt                          <  1:85478/100‑1 (MQ=255)
                                                                     gAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCt                          <  2:1983342/100‑1 (MQ=255)
                                                                     gAATGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCt                          >  1:593629/1‑100 (MQ=255)
                                                                       aTGGTATTAAGCTGCACGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTa                         >  1:734344/1‑99 (MQ=255)
                                                                       aTGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAc                        <  2:1881194/100‑1 (MQ=255)
                                                                       aTGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAc                        <  1:1659173/100‑1 (MQ=255)
                                                                       aTGGTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATAAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAc                        <  1:613605/100‑1 (MQ=255)
                                                                         ggTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtg                      <  2:569085/100‑1 (MQ=255)
                                                                         ggTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtg                      <  2:360051/100‑1 (MQ=255)
                                                                         ggTATTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTCAAAAACACTGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtg                      <  2:928222/100‑1 (MQ=255)
                                                                            aTTAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtgttg                   <  1:537650/100‑1 (MQ=255)
                                                                              tAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTc                 >  1:1468629/1‑100 (MQ=255)
                                                                              tAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTc                 >  2:282430/1‑100 (MQ=255)
                                                                              tAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTc                 <  1:71173/100‑1 (MQ=255)
                                                                              tAAGCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTc                 <  2:1228883/100‑1 (MQ=255)
                                                                                 gCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCg              <  2:1930291/100‑1 (MQ=255)
                                                                                 gCTGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCg              <  1:2030331/100‑1 (MQ=255)
                                                                                   tGCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGtc            >  2:1286063/1‑100 (MQ=255)
                                                                                    gCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGtct           >  1:481242/1‑100 (MQ=255)
                                                                                    gCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGtct           <  1:167877/100‑1 (MQ=255)
                                                                                    gCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGtct           >  1:645267/1‑100 (MQ=255)
                                                                                    gCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGtct           >  2:1430333/1‑100 (MQ=255)
                                                                                    gCAAGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGtct           >  1:1176507/1‑100 (MQ=255)
                                                                                       aGGGCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACACGCCCGCGATTTCTACTGTTGTCCCGTTTCGc        >  2:741290/1‑100 (MQ=255)
                                                                                        gggCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGcc       >  1:445453/1‑100 (MQ=255)
                                                                                        gggCAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGcc       >  1:1230686/1‑100 (MQ=255)
                                                                                           cAAATCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCgg    >  1:588399/1‑100 (MQ=255)
                                                                                             aaTCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGtt  <  1:223888/100‑1 (MQ=255)
                                                                                             aaTCAAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGtt  <  2:371626/100‑1 (MQ=255)
                                                                                                 |                                                                                               
GCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTT  >  NC_000913/4400299‑4400491

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