Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,533,930 T→G 6.2% intergenic (‑79/‑94) yddH ← / → nhoA flavin reductase like‑protein/N‑hydroxyarylamine O‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,533,9300TG 6.2% 205.1 / 5.6 65intergenic (‑79/‑94)yddH/nhoAflavin reductase like‑protein/N‑hydroxyarylamine O‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (26/35);  new base G (2/2);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGAAACGCCTGATGACGCCCATTCTGAATCACTATTTTGCCCGTATTAACTGGTCGGGAGCT  >  NC_000913/1533786‑1534074
                                                                                                                                                |                                                                                                                                                
gACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCg                                                                                                                                               <  5:99281/148‑1 (MQ=255)
 aCCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCg                                                                                                                                               >  3:75365/1‑147 (MQ=255)
    gTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACg                                                                                                                                           <  6:142242/148‑1 (MQ=255)
         gCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATg                                                                                                                                      <  3:200651/148‑1 (MQ=255)
          ccATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGc                                                                                                                                     <  1:253880/148‑1 (MQ=255)
             tGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATg                                                                                                                                  <  1:198834/148‑1 (MQ=255)
                 ttAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCcg                                                                                                                              >  3:36768/1‑148 (MQ=255)
                    aCAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATCGCCTCGCGAAGAAGGGTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCcgcg                                                                                                                            <  3:268053/147‑1 (MQ=255)
                      aGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATc                                                                                                                         >  6:233534/1‑148 (MQ=255)
                       gACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCg                                                                                                                        >  5:393572/1‑148 (MQ=255)
                         cGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGca                                                                                                                      >  3:369203/1‑148 (MQ=255)
                            cTGGCATGGTGTAATTCTATGGGGATGACTCGGCTCATGGTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGcaca                                                                                                                    >  8:214910/1‑147 (MQ=255)
                                  tGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGc                                                                                                             >  8:87912/1‑148 (MQ=255)
                                    gtgtAATTCTATGGGGATGAATCGGCTCATGGTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCTATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCg                                                                                                           <  4:262229/148‑1 (MQ=255)
                                          ttCTATGGGGATGAATCGGCTCATGGTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGAt                                                                                                                          >  1:328054/1‑127 (MQ=255)
                                          ttCTATGGGGATGAATCGGCTCATGGTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGAt                                                                                                                          <  2:328054/127‑1 (MQ=255)
                                               tGGGGATGAATCGGCTCATGGTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTgcggcg                                                                                                   <  7:236912/145‑1 (MQ=255)
                                                  ggATGAATCGGCTCATGGTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTAtt                                                                                             <  4:369203/148‑1 (MQ=255)
                                                               cATGGTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGt                                                                                >  6:118399/1‑148 (MQ=255)
                                                               cATGGTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGt                                                                                <  7:214910/148‑1 (MQ=255)
                                                                 tagTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTg                                                                              <  3:135412/146‑1 (MQ=255)
                                                                 tGGTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGtt                                                                               >  8:267258/1‑147 (MQ=255)
                                                                 tGGTAGCGCCTCGCGAAGAAGGTGTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTg                                                                              <  2:403921/148‑1 (MQ=255)
                                                                  ggTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGtt                                                                               >  1:425392/1‑146 (MQ=255)
                                                                         ccTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTg                                                                                                                                         <  6:137067/81‑1 (MQ=255)
                                                                         ccTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTg                                                                                                                                         >  5:137067/1‑81 (MQ=255)
                                                                              cgAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTaa                                                                 <  7:87912/148‑1 (MQ=255)
                                                                                   aagGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGaa                                                            <  1:312108/148‑1 (MQ=255)
                                                                                       tttttttAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGGGTCTATTCCAGGTTATAAg                                                                                 >  7:118308/1‑123 (MQ=255)
                                                                                       tttttttAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAg                                                                                 <  8:118308/123‑1 (MQ=255)
                                                                                         tttttAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGAAACGcc                                                       <  5:118399/147‑1 (MQ=255)
                                                                                         tttttAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGAAACGcc                                                       <  2:89653/147‑1 (MQ=255)
                                                                                          ttttAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGAAACGCCt                                                      <  5:233534/147‑1 (MQ=255)
                                                                                              aaGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGAAACGCCtgatg                                                  <  4:75365/147‑1 (MQ=255)
                                                                                              aaGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGAAACGCCtgatg                                                  <  7:267258/147‑1 (MQ=255)
                                                                                              aaGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGAAACGCCtgatg                                                  <  2:425392/147‑1 (MQ=255)
                                                                                                    aGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGAAACGCCTGATGACGCCCa                                           >  7:270188/1‑148 (MQ=255)
                                                                                                      tCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCa                                                                     <  5:129650/120‑1 (MQ=255)
                                                                                                      tCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCa                                                                     >  6:129650/1‑120 (MQ=255)
                                                                                                       ccGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTgcgg                                                                                                     <  8:32590/87‑1 (MQ=255)
                                                                                                       ccGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTgcgg                                                                                                     >  7:32590/1‑87 (MQ=255)
                                                                                                         gTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTTTAAGTTGAGAAAACCACTAAGGGAACCGCCTGATGACGCCCATTCt                                       >  8:384009/1‑147 (MQ=255)
                                                                                                         gTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGAAACGCCTGATGACGCCCATTCt                                       >  8:435468/1‑147 (MQ=255)
                                                                                                            aCGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGAAACGCCTGATGACGCCCATTCTGaa                                    <  6:393572/147‑1 (MQ=255)
                                                                                                                  tAAGTAACGAAATTAACTGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCTGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGAAACGCCTGATGACGCCCATTCTGAATCACTAt                             <  3:27384/148‑1 (MQ=255)
                                                                                                                  tAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGAAACGCCTGATGACGCCCATTCTGAATCACTAt                             <  4:36768/148‑1 (MQ=255)
                                                                                                                   aaGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGaa                                                            <  3:106070/116‑1 (MQ=255)
                                                                                                                   aaGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGaa                                                            >  4:106070/1‑116 (MQ=255)
                                                                                                                   aaGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGAAACGCCTGATGACGCCCATTCTg                                      <  6:396303/138‑1 (MQ=255)
                                                                                                                   aaGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGAAACGCCTGATGACGCCCATTCTg                                      >  5:396303/1‑138 (MQ=255)
                                                                                                                         cGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGt                                                                                <  2:316769/90‑1 (MQ=255)
                                                                                                                         cGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGt                                                                                >  1:316769/1‑90 (MQ=255)
                                                                                                                         cGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGAAACGCCTGATGACGCCCATTCTGAATCACTATTTTGccc                      >  1:283201/1‑148 (MQ=255)
                                                                                                                             aTTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCACACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGAAACGCCTGATGACGCCCATTCTGAATCACTATTTTGCCCGTAt                  <  7:259831/148‑1 (MQ=255)
                                                                                                                              ttAACGGGATTGGCGATTGGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAg                                                                >  1:159620/1‑101 (MQ=255)
                                                                                                                              ttAACGGGATTGGCGATTGGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAg                                                                <  2:153561/101‑1 (MQ=255)
                                                                                                                              ttAACGGGATTGGCGATTGGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAg                                                                <  2:159620/101‑1 (MQ=255)
                                                                                                                              ttAACGGGATTGGCGATTGGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAg                                                                >  1:153561/1‑101 (MQ=255)
                                                                                                                                 aCGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGAAACGCCTGATGACGCCCATTCTGAATCACTATTTTGCCCGTATTaa               <  7:435468/147‑1 (MQ=255)
                                                                                                                                  cGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGAAACGCCTGATGACGCCCATTCTGAATCACTATTTTGCCCGTATTAACt             >  1:193543/1‑148 (MQ=255)
                                                                                                                                        tGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGAAACGCCTGATGACGCCCATTCTGAATCACTATTTTGCCCGTATTAACTGGTCgg       >  5:332389/1‑148 (MQ=255)
                                                                                                                                         ggCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGAAACGCCTGATGACGCCCATTCTGAATCACTATTTTGCCCGTATTAACTGGTCgg       >  6:100055/1‑147 (MQ=255)
                                                                                                                                          gCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGCGAAACGCCTGATGACGCCCATTCTGAATCACTATTTTGCCCGTATTAACTGGTCggg      <  4:14216/147‑1 (MQ=255)
                                                                                                                                            gATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGAAACGCCTGATGACGCCCATTCTGAATCACTATTTTGCCCGTATTAACTGGTCGGGAgc   <  8:270188/148‑1 (MQ=255)
                                                                                                                                             aTTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAACTTGAGAAAACCACGAAGGGAAACGCCTGATGACGCCCATTCTGAATCACTATTTTGCCCGTATTAACTGGTCGGGAgct  >  3:209024/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATCGCCTCGCGAAGAAGGTTTTTTTAAGCGTAGTCCGTAACGCAATAAGTAACGAAATTAACGGGATTGGCGATTTGCGAACGTGATGCATGTCCGCGATCGCACAAAATAGCCGGTGCGGCGTCTATTCCAGGTTATAAGTTGAGAAAACCACTAAGGGAAACGCCTGATGACGCCCATTCTGAATCACTATTTTGCCCGTATTAACTGGTCGGGAGCT  >  NC_000913/1533786‑1534074

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