Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,151,418 A→G intergenic (+14/‑197) fabG → / → acpP 3‑oxoacyl‑[acyl‑carrier‑protein] reductase/acyl carrier protein (ACP)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,151,4180AG98.4% 216.8 / ‑5.2 62intergenic (+14/‑197)fabG/acpP3‑oxoacyl‑[acyl‑carrier‑protein] reductase/acyl carrier protein (ACP)
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (27/34);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTT  >  NC_000913/1151272‑1151562
                                                                                                                                                  |                                                                                                                                                
ggCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACgaaa                                                                                                                                               <  1:909601/150‑1 (MQ=255)
 gcgcAGGTTCCTGCGGGTCGCCGCGGCGGGGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCAGACGGAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAAAGGCGGGAGGTACATGGTCTGACCGCGATTTGCACGAAAt                                                                                                                                              >  2:317940/1‑150 (MQ=255)
             gCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTgcgc                                                                                                                                 <  2:567633/151‑1 (MQ=255)
                ggTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTgcgcgc                                                                                                                               >  2:429523/1‑150 (MQ=255)
                 gTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAg                                                                                                                             <  1:100775/151‑1 (MQ=255)
                   cGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTc                                                                                                                           >  2:665125/1‑151 (MQ=255)
                     ccTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTg                                                                                                                         <  1:727677/151‑1 (MQ=255)
                         ggcggcGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCCCGAAATGCTCATGTTGCGCGCAGTCTGCg                                                                                                                       >  2:185651/1‑149 (MQ=255)
                         ggcggcGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTg                                                                                                                     >  2:897375/1‑151 (MQ=255)
                          gcggcgCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTgg                                                                                                                    <  2:569170/151‑1 (MQ=255)
                                cacaGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGa                                                                                                              <  2:738977/151‑1 (MQ=255)
                                 acaGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAg                                                                                                             <  1:634640/151‑1 (MQ=255)
                                 acaGGAAAGAGCAAACGCGGTTGCAATGCTGGCATCCGACGAAGCAGATTACGTCACGGGTGAAAATTTGCATGTGAACGGCGGGAGGGACATGGTCTTACCGCGATTTGCACAAAATGCTCAGGATGGGCGCAGGCTGCGGGGTTAt                                                                                                                >  2:367608/1‑148 (MQ=255)
                                    ggAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGtaa                                                                                                          <  1:539231/151‑1 (MQ=255)
                                     gAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGtaat                                                                                                         >  2:901859/1‑151 (MQ=255)
                                       aaTCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGATtaataa                                                                                                       <  1:777755/151‑1 (MQ=255)
                                         tCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAAtt                                                                                                     >  1:337779/1‑151 (MQ=255)
                                          cGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTa                                                                                                    <  1:386340/151‑1 (MQ=255)
                                           gCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAAtt                                                                                                     >  1:491306/1‑149 (MQ=255)
                                                   ggTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAAt                                                                                           >  1:835833/1‑151 (MQ=255)
                                                         aTTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTg                                                                                     <  1:131637/151‑1 (MQ=255)
                                                         aTTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTg                                                                                     <  2:713073/151‑1 (MQ=255)
                                                         aTTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTg                                                                                     <  1:283541/151‑1 (MQ=255)
                                                            ccTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGt                                                                                  <  2:562379/151‑1 (MQ=255)
                                                               ggCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTa                                                                                <  1:146065/150‑1 (MQ=255)
                                                               ggCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTAt                                                                               <  1:400899/151‑1 (MQ=255)
                                                                    ccGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTggg                                                                           <  2:39260/150‑1 (MQ=255)
                                                                    ccGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTggg                                                                           <  1:597936/150‑1 (MQ=255)
                                                                       acgaAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGt                                                                       >  2:1090/1‑151 (MQ=255)
                                                                          aaGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCTGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGGTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAg                                                                    <  2:11668/151‑1 (MQ=255)
                                                                            gcagcTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGc                                                                  <  1:429569/151‑1 (MQ=255)
                                                                               gcTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCt                                                                <  1:332581/150‑1 (MQ=255)
                                                                                  tACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCaa                                                             <  1:245490/150‑1 (MQ=255)
                                                                                   aCATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCaaa                                                            >  2:237906/1‑150 (MQ=255)
                                                                                     aTCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAAt                                                          >  2:859844/1‑150 (MQ=255)
                                                                                          gggTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTa                                                    <  1:379350/151‑1 (MQ=255)
                                                                                           ggTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATATGTGTAAAATGATTTGCGTTATTGGGGAGTAAGGCCTCAAAATAACg                                                      <  2:469514/148‑1 (MQ=255)
                                                                                              gAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAAt                                                >  1:646416/1‑151 (MQ=255)
                                                                                                 aCTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGt                                             >  1:101290/1‑151 (MQ=255)
                                                                                                   tttGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGGGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTgg                                           >  1:568964/1‑151 (MQ=255)
                                                                                                       cATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAg                                       <  1:901932/151‑1 (MQ=255)
                                                                                                             aaCGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGc                                 <  2:119929/151‑1 (MQ=255)
                                                                                                              aCGGCGGGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGcc                                >  1:39064/1‑151 (MQ=255)
                                                                                                               cggcggGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCg                               >  2:355567/1‑151 (MQ=255)
                                                                                                               cggcggGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCg                               <  2:257052/151‑1 (MQ=255)
                                                                                                                ggcggGATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCgg                              >  2:256423/1‑151 (MQ=255)
                                                                                                                   gggATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGt                                                     <  1:341560/125‑1 (MQ=255)
                                                                                                                   gggATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGt                                                     >  2:341534/1‑125 (MQ=255)
                                                                                                                   gggATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGAt                           >  2:401615/1‑151 (MQ=255)
                                                                                                                    ggATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGAtt                          <  2:852476/151‑1 (MQ=255)
                                                                                                                     gATGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCCGCCGGGAttt                         >  1:635026/1‑151 (MQ=255)
                                                                                                                      aTGTACATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTa                        <  1:59835/151‑1 (MQ=255)
                                                                                                                          aCATGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTg                    <  1:897448/151‑1 (MQ=255)
                                                                                                                            aTGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGc                   >  2:295830/1‑150 (MQ=255)
                                                                                                                             tGGTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCa                  >  2:713498/1‑150 (MQ=255)
                                                                                                                              ggTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCaa                 >  1:572422/1‑150 (MQ=255)
                                                                                                                              ggTCTGACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGCGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCaa                 >  2:698869/1‑150 (MQ=255)
                                                                                                                                   gACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAAttttt           <  1:377820/151‑1 (MQ=255)
                                                                                                                                   gACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAAttttt           >  1:583139/1‑151 (MQ=255)
                                                                                                                                   gACCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAAtt              <  2:149312/148‑1 (MQ=255)
                                                                                                                                    aCCGCGATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTc          <  1:20105/151‑1 (MQ=255)
                                                                                                                                       gcgATTTGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAAc       >  2:482310/1‑151 (MQ=255)
                                                                                                                                             tGCACGAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACAtttt  <  2:881272/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTT  >  NC_000913/1151272‑1151562

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: