Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 584,262 (A)7→8 coding (582/750 nt) appY → global transcriptional activator; DLP12 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913584,2551.A89.2% 239.9 / 22.1 75coding (575/750 nt)appYglobal transcriptional activator; DLP12 prophage
Reads supporting (aligned to +/- strand):  ref base . (3/5);  new base A (41/25);  total (45/30)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGAATTGATTTATACGAGTGAAAGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑C‑AAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTA  >  NC_000913/584156‑584347
                                                                                                      |                                                                                             
cGGAATTGATTTATACGAGTGAAAGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑c                                                                                                >  2:846530/1‑100 (MQ=255)
cGGAATTGATTTATACGAGTGAAAGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑c                                                                                                <  1:986738/100‑1 (MQ=255)
 ggAATTGATTTATACGAGTGAAAGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑C‑a                                                                                              <  2:806518/100‑1 (MQ=255)
     ttGATTTATACGAGTGAAAGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑C‑aaaaa                                                                                          <  2:1632474/100‑1 (MQ=255)
      tGATTTATACGAGTGAAAGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑C‑aaaaa‑a                                                                                        >  2:179403/1‑100 (MQ=255)
       gATTTATACGAGTGAAAGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑C‑aaaaa‑aa                                                                                       >  1:1486692/1‑100 (MQ=255)
        aTTTATACGAGTGAAAGTTTAAT‑AAAAAAAAGCTTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑C‑aaaaa‑aa                                                                                       >  1:228159/1‑99 (MQ=255)
        aTTTATACGAGTGAAAGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑C‑CAAAA‑AAc                                                                                      <  2:1025106/100‑1 (MQ=255)
          ttATACGAGTGAAAGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑C‑aaaaa‑aaact                                                                                    <  1:791840/100‑4 (MQ=255)
           tataCGAGTGAAAGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGAAACGTCATTTACTGAAATATTGAGAGATACTAGGATTAGGTATA‑C‑aaaaa‑aaac                                                                                     >  1:189133/1‑96 (MQ=255)
            ataCGAGTGAAAGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑C‑AAAAA‑AACTCAt                                                                                  >  2:1764141/1‑100 (MQ=255)
              aCGAGTGAAAGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATa                                                                                 >  2:501408/1‑100 (MQ=255)
              aCGAGTGAAAGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATa                                                                                 <  2:1091852/100‑1 (MQ=255)
                gAGTGAAAGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAAc                                                                               >  2:369536/1‑100 (MQ=255)
                 aGTGAAAGTTTAATAAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTAAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAAc                                                                               >  1:230420/1‑100 (MQ=255)
                 aGTGAAAGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACt                                                                              >  2:1592712/1‑100 (MQ=255)
                  gTGAAAGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACtt                                                                             >  2:1257018/1‑100 (MQ=255)
                  gTGAAAGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACtt                                                                             >  1:1749369/1‑100 (MQ=255)
                  gTGAAAGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACtt                                                                             >  1:1390392/1‑100 (MQ=255)
                  gTGAAAGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑C‑AAAAA‑AACTCATAACTTc                                                                            <  2:800103/100‑1 (MQ=255)
                   tGAAAGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTc                                                                            >  1:793869/1‑100 (MQ=255)
                     aaaaTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAAAGTCATTTACTGAAATATTGAGAGATACTAGGATGAAAAATA‑AAAAAAA‑AACTCATAACTTCaa                                                                          <  1:1592701/100‑1 (MQ=255)
                     aaaGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCa                                                                           >  2:1191614/1‑99 (MQ=255)
                       aGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAAc                                                                        <  2:179708/100‑1 (MQ=255)
                             at‑aaAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTAt                                                                  >  1:1789198/1‑100 (MQ=255)
                             at‑aaAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTAt                                                                  >  2:100306/1‑100 (MQ=255)
                             at‑aaAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑C‑‑AAAA‑AACTCATAACTTCAAACTCTTaat                                                                 >  2:1724636/1‑97 (MQ=255)
                                aaaaaaaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTc                                                                >  1:1769070/1‑100 (MQ=255)
                                 aaaaaaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCt                                                               <  1:880753/100‑1 (MQ=255)
                                  aaaaaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTa                                                              <  2:1033969/100‑1 (MQ=255)
                                  aaaaaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑C‑AAAAATAACTCATAGCTTCAAACTCTTATTc                                                                <  1:850676/98‑1 (MQ=255)
                                   aaaaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTAt                                                             >  1:1904689/1‑100 (MQ=255)
                                      aaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATAAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCaa                                                          <  2:1485630/100‑1 (MQ=255)
                                        gATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATg                                                        >  1:506822/1‑100 (MQ=255)
                                         aTTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑C‑AAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTc                                                      >  2:1776144/1‑100 (MQ=255)
                                          ttAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTc                                                      >  2:1595581/1‑100 (MQ=255)
                                          ttAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTc                                                      >  1:1942630/1‑100 (MQ=255)
                                           tAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCg                                                     <  1:193264/100‑1 (MQ=255)
                                              gggATGAAGGAACGTCATTTACTGGAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAg                                                  >  2:986520/1‑100 (MQ=255)
                                                gATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAgac                                                >  1:1181635/1‑98 (MQ=255)
                                                gATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGcc                                                >  1:617360/1‑100 (MQ=255)
                                                gATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGcc                                                <  2:1670184/100‑1 (MQ=255)
                                                gATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGcc                                                <  2:1579722/100‑1 (MQ=255)
                                                 aTGAAGGAACATCATTTACTTAAATATTGAGATATTCTAGGATGAAAAAAA‑AAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGccc                                               <  1:407669/100‑1 (MQ=255)
                                                  tGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCa                                              <  2:1594810/100‑1 (MQ=255)
                                                   gAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAg                                             <  2:295198/100‑1 (MQ=255)
                                                    aaGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGa                                            <  1:519867/100‑1 (MQ=255)
                                                       gAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAAt                                         >  2:63408/1‑100 (MQ=255)
                                                        aaCGTCATTTACTGAAATATTGAGAGATACTAGGATGAAGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAAtg                                        <  1:1791536/100‑1 (MQ=255)
                                                         acatCATTTACTGAAATATTGAGAGATACTAGGAAAAAAAAAAAA‑AAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGaaa                                          <  1:425863/94‑1 (MQ=255)
                                                          cGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAAtgtg                                      <  1:716675/100‑1 (MQ=255)
                                                          cGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAAtgtg                                      >  2:1824990/1‑100 (MQ=255)
                                                          cGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑C‑AAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTgg                                     <  1:1055338/100‑1 (MQ=255)
                                                           gTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTgg                                     >  2:1238693/1‑100 (MQ=255)
                                                            tCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGc                                    <  2:1256/100‑1 (MQ=255)
                                                              aTTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCta                                  >  2:1659943/1‑100 (MQ=255)
                                                                 tACTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATaa                               >  2:1821843/1‑100 (MQ=255)
                                                                 tACTGAAATATTGAGAGATACTAGGATGAAAAATA‑AAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCtata                                <  1:216027/99‑1 (MQ=255)
                                                                  aCTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAAc                              <  2:1749380/100‑1 (MQ=255)
                                                                   cTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACa                             >  2:277733/1‑100 (MQ=255)
                                                                   cTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACa                             <  2:91001/100‑1 (MQ=255)
                                                                   cTGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCACAGAAATGTGGCTATAACa                             >  1:1058145/1‑100 (MQ=255)
                                                                    tGAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAg                            <  2:591199/100‑1 (MQ=255)
                                                                     gAAATATTGAGAGATACTAGGATGAGGTATG‑CGAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGt                           >  1:564530/1‑100 (MQ=255)
                                                                     gAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGGAATGTGGCTATAACAGt                           >  2:1899066/1‑100 (MQ=255)
                                                                     gAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGt                           <  2:55084/100‑1 (MQ=255)
                                                                     gAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGt                           >  2:1682325/1‑100 (MQ=255)
                                                                     gAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGt                           >  2:178994/1‑100 (MQ=255)
                                                                     gAAATATTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCGTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGt                           >  1:1364871/1‑100 (MQ=255)
                                                                         tatTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACtt                       >  1:38072/1‑100 (MQ=255)
                                                                          atTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttc                      >  1:751394/1‑100 (MQ=255)
                                                                          atTGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttc                      <  1:1606097/100‑1 (MQ=255)
                                                                           ttGAGAGATACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttca                     >  1:257402/1‑100 (MQ=255)
                                                                             gagagaTACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcata                   >  1:878619/1‑100 (MQ=255)
                                                                             gagagaTACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcata                   >  1:1541166/1‑100 (MQ=255)
                                                                                   tACTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttcat             >  1:713134/1‑100 (MQ=255)
                                                                                    aCTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttcata            >  1:1467569/1‑100 (MQ=255)
                                                                                    aCTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttcata            <  2:564766/100‑1 (MQ=255)
                                                                                    aCTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttcata            <  2:1876485/100‑1 (MQ=255)
                                                                                     cTAGGATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttcatat           >  2:1747015/1‑100 (MQ=255)
                                                                                        ggATGAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATAtgtg        <  1:1075271/100‑1 (MQ=255)
                                                                                            gAGGTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCAtt    >  1:781674/1‑100 (MQ=255)
                                                                                              ggTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTa  >  1:1495276/1‑100 (MQ=255)
                                                                                              ggTATG‑CAAAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTa  >  2:958330/1‑100 (MQ=255)
                                                                                                      |                                                                                             
CGGAATTGATTTATACGAGTGAAAGTTTAAT‑AAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATG‑C‑AAAAA‑AACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTA  >  NC_000913/584156‑584347

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: ‑