Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,416,948 A→G 33.6% intergenic (‑178/‑4) bsmA ← / → yjfP bioflm peroxide resistance protein/acyl CoA esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,416,9480AG33.6% 99.3 / 88.9 104intergenic (‑178/‑4)bsmA/yjfPbioflm peroxide resistance protein/acyl CoA esterase
Reads supporting (aligned to +/- strand):  ref base A (34/35);  new base G (19/16);  total (53/51)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

Only 100 of 106 total aligned reads displayed.

CGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTT  >  NC_000913/4416852‑4417047
                                                                                                |                                                                                                   
cGAAAGCGTAACCCACATCTCACCATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAATGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAg                                                                                                  >  2:1301499/1‑100 (MQ=255)
cGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAg                                                                                                  >  2:785290/1‑100 (MQ=255)
cGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAg                                                                                                  >  2:1258721/1‑100 (MQ=255)
cGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAg                                                                                                  >  2:469740/1‑100 (MQ=255)
cGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAg                                                                                                  >  2:2015199/1‑100 (MQ=255)
   aaGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAgatg                                                                                               >  2:1542304/1‑100 (MQ=255)
     gCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAgatgat                                                                                             <  2:21213/100‑1 (MQ=255)
     gCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAgatgat                                                                                             >  2:766704/1‑100 (MQ=255)
        tAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGa                                                                                          >  2:1015062/1‑100 (MQ=255)
          aCCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGaaa                                                                                        >  1:545192/1‑100 (MQ=255)
           cccACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGAAAAAAAGGAGAGACAGATGATTGAAAt                                                                                       <  2:1046257/100‑1 (MQ=255)
            ccACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAGATGATTGAAATa                                                                                      >  2:462719/1‑100 (MQ=255)
            ccACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAGATGATTGAAATa                                                                                      >  2:462751/1‑100 (MQ=255)
             cacaTCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAGATGATTGAAATan                                                                                     <  1:822525/100‑2 (MQ=255)
              acaTCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAGATGATTGAAATAGa                                                                                    <  1:1006896/100‑1 (MQ=255)
               caTCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGaa                                                                                   <  2:1319906/100‑1 (MQ=255)
                aTCTCAACATTATGAACTAATACCGCCAGCAGAAACTGGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGaa                                                                                   <  1:458635/99‑1 (MQ=255)
                 tttcAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATc                                                                                 >  1:1522472/3‑100 (MQ=255)
                 tctcAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATc                                                                                 <  1:1189578/100‑1 (MQ=255)
                  ctcAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCa                                                                                >  2:822373/1‑100 (MQ=255)
                  ctcAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCa                                                                                <  1:155769/100‑1 (MQ=255)
                   tcAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAGATGATTGAAATAGAATCAc                                                                               <  2:277122/100‑1 (MQ=255)
                    cAACATTATGAACTAATACCGCCAGCATAAACTTGTGCTAAGGATTTCTTCATGTAAAAAAAAAAAAAAAGGAGAGACAGATGATTGAAATAGaa                                                                                   <  2:150412/95‑1 (MQ=255)
                    cAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCAcg                                                                              >  2:2089344/1‑100 (MQ=255)
                    cAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAAGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCAcg                                                                              >  2:1226189/1‑100 (MQ=255)
                        aTTATGAACTAATACCGCCAGCAGAAACTTGTGTTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCAcgc                                                                             <  1:1538593/97‑1 (MQ=255)
                         ttATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAGATGATTGAAATAGAATCACGCGAGc                                                                         >  1:250166/1‑100 (MQ=255)
                           aTGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTg                                                                       >  2:568366/1‑100 (MQ=255)
                           aTGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTg                                                                       >  2:335284/1‑100 (MQ=255)
                             gAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGc                                                                     <  2:1394099/100‑1 (MQ=255)
                                cTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGa                                                                  >  2:58159/1‑100 (MQ=255)
                                cTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGa                                                                  <  2:1808983/100‑1 (MQ=255)
                                 tAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAGATGATTGAAATAGAATCACGCGAGCTGGCAGat                                                                 <  2:777224/100‑1 (MQ=255)
                                      ccGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTcc                                                            <  1:314520/100‑1 (MQ=255)
                                      ccGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTcc                                                            <  1:523176/100‑1 (MQ=255)
                                      ccGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTcc                                                            <  1:1608510/100‑1 (MQ=255)
                                      ccGCCAGCAGAAACTTGTGCTAAGGATTTCTTAATGTCATACCAGAAAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTcc                                                            <  2:1354925/100‑1 (MQ=255)
                                       cGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAGAGGATTGAAATAGAATCACGCGAGCTGGCAGATATTccc                                                           <  2:1195052/100‑1 (MQ=255)
                                         ccAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGt                                                         <  2:47486/100‑1 (MQ=255)
                                           agcagAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTACCGttc                                                       >  2:513580/1‑100 (MQ=255)
                                            gcagAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGttct                                                      <  2:1275888/100‑1 (MQ=255)
                                             cagAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGttctt                                                     >  1:379727/1‑100 (MQ=255)
                                               gAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCa                                                   >  2:635462/1‑100 (MQ=255)
                                                aaaaTTGTGCTAAGGATTTCTTCATGAAAAAAAAAAAAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGttcttc                                                    <  2:782328/98‑1 (MQ=255)
                                                  aCTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGc                                                <  2:1484813/100‑1 (MQ=255)
                                                  aCTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGc                                                <  1:2231504/100‑1 (MQ=255)
                                                  aCTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGc                                                <  1:330649/100‑1 (MQ=255)
                                                  aCTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGc                                                <  1:1848802/100‑1 (MQ=255)
                                                   cTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAGATGATTGAAATAGAATCACGTGAGCTGGCAGATATTCCCGTTCTTCATGc                                                <  2:1484428/99‑1 (MQ=255)
                                                   cTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGc                                                >  1:1168993/1‑99 (MQ=255)
                                                   cTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCt                                               <  2:851340/100‑1 (MQ=255)
                                                   cTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCt                                               <  2:790284/100‑1 (MQ=255)
                                                    ttGTGTTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCg                                                          <  1:78489/88‑1 (MQ=255)
                                                      gtgCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCCCGCGCGCTGGCAGATATTCCCGTTCTTCATGc                                                <  2:1565806/96‑1 (MQ=255)
                                                       tgCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTAt                                            >  1:1984926/1‑99 (MQ=255)
                                                        gCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATcc                                          >  1:1399156/1‑100 (MQ=255)
                                                        gCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATcc                                          >  2:7582/1‑100 (MQ=255)
                                                        gCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATcc                                          <  2:717233/100‑1 (MQ=255)
                                                        gCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATcc                                          >  2:103678/1‑100 (MQ=255)
                                                        gCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATcc                                          >  2:802566/1‑100 (MQ=255)
                                                        gCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATcc                                          <  2:545193/100‑1 (MQ=255)
                                                         cTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCt                                         >  2:109962/1‑100 (MQ=255)
                                                         cTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCt                                         >  1:1032916/1‑100 (MQ=255)
                                                         cTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCt                                         >  1:822429/1‑100 (MQ=255)
                                                           aaGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGt                                       >  1:1398254/1‑100 (MQ=255)
                                                            aGGATTTCTTCATGTCATACCATACAAAAACGAGAGACAGATGATTGAAATAGAATCACGCGAGCAGGCAGATATTCCCGTTCTTCATGCTTATcc                                          <  2:2082242/96‑1 (MQ=255)
                                                            aGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTc                                      <  1:109962/100‑1 (MQ=255)
                                                            aGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGt                                       >  1:2246481/1‑99 (MQ=255)
                                                             ggATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCg                                     >  2:1359968/1‑100 (MQ=255)
                                                             ggATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTTTCg                                     >  2:199776/1‑100 (MQ=255)
                                                              gATTTCTTCATGTCATACCAGACAAAAAGGAGAGGCAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCgg                                    >  1:486924/1‑100 (MQ=255)
                                                              gATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCgg                                    >  2:1837781/1‑100 (MQ=255)
                                                              gATTTCTTCATGTAAAAAAAAAAAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCgg                                    <  2:2183291/100‑1 (MQ=255)
                                                               atgtgttcATGTCATACCAGACAAAAGTGAGTGGTGGATGATTGAACTAGAATCACACGAGCTGGCAGAGATTCCCGTTCTTCATGGTTATCCTGTCggg                                   <  2:2160678/95‑1 (MQ=255)
                                                                tttCTTCATGTCATACCAGACAAAAAGGAGAGGCAGATGATTGAAATAGAATCACACGAGCTGTCAGATATTCCCGTTCTTCATGCTTATCCTGTCggg                                   >  2:30089/1‑99 (MQ=255)
                                                                tttCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGc                                  <  1:1091055/100‑1 (MQ=255)
                                                                 ttcttcATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCa                                 <  1:237959/100‑1 (MQ=255)
                                                                 ttcttcATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCa                                 <  1:1066970/100‑1 (MQ=255)
                                                                   cttcATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCaaa                               >  2:286660/1‑100 (MQ=255)
                                                                    ttcATGTCATACCAGACAAAAAGGAGAGGCAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCaaaa                              >  2:924257/1‑100 (MQ=255)
                                                                      cATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAg                            <  1:828879/100‑1 (MQ=255)
                                                                      cATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAg                            <  2:1093514/100‑1 (MQ=255)
                                                                       aTGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGa                           >  1:2029579/1‑100 (MQ=255)
                                                                         gTCATACCAGACAAAAAGGAGAGGCAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATa                         <  2:1400632/100‑1 (MQ=255)
                                                                         gTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATa                         >  2:895577/1‑100 (MQ=255)
                                                                           cATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATAcc                       >  1:319612/1‑100 (MQ=255)
                                                                            aTACCAGACAAAAAGGAGAGGCAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATa                         >  1:860459/1‑97 (MQ=255)
                                                                            aTACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATAccc                      >  2:877825/1‑100 (MQ=255)
                                                                            aTACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATAccc                      <  1:448055/100‑1 (MQ=255)
                                                                               ccAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGt                   <  1:672258/100‑1 (MQ=255)
                                                                                 aGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTa                 <  1:1909841/100‑1 (MQ=255)
                                                                                  gACAAAAAGGAGAGGCAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTAc                >  2:1021090/1‑100 (MQ=255)
                                                                                   acacataGGAAAGGCAGATGATTAAAATAGAATCACGCGAGCTGGCAGATATTCCCGTACTTCATGCTTATCATGTCGGGCAAAAAGATACACCGTTAc                >  2:931367/7‑99 (MQ=255)
                                                                                   aCAAAAAGGAGAGACAGCTGATTGAAATACAAACACGAGAGCTGGCAGATACCCCCGATCTTCATGCTTATTCTGTCCGGCAAAAAGATACCCCGTTAcc               >  2:63788/1‑100 (MQ=255)
                                                                                    cAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCg              >  1:621497/1‑100 (MQ=255)
                                                                                      aaaaGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTg            >  1:1489872/1‑100 (MQ=255)
                                                                                       aaaGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGc           <  1:1738862/100‑1 (MQ=255)
                                                                                        aaGGAGAGGCAGATGATTGAAATAGAATCACGCGAGCGGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCg          <  1:765833/100‑1 (MQ=255)
                                                                                              agaCAGATGATTGAAATAGAATCACGCGGGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAAttt    <  1:377843/100‑1 (MQ=255)
                                                                                                aCAGATGATTGAAAGAGAATCACGCGAGCTGGCAGATATCCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAAttttt  <  1:1821024/100‑1 (MQ=255)
                                                                                                |                                                                                                   
CGAAAGCGTAACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTT  >  NC_000913/4416852‑4417047

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