Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,736,100 T→A 15.6% intergenic (+141/+21) sodB → / ← ydhP superoxide dismutase, Fe/putative MFS transporter, inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,736,1000TA15.6% 248.1 / 35.9 96intergenic (+141/+21)sodB/ydhPsuperoxide dismutase, Fe/putative MFS transporter, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (41/40);  new base A (7/8);  total (48/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 = 9.51e-01

TCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGAC  >  NC_000913/1736002‑1736198
                                                                                                  |                                                                                                  
tCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATc                                                                                                   >  1:1175142/1‑100 (MQ=255)
tCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTAAGCAACTGTAACGAAGAAGGTTATc                                                                                                   >  1:1175090/1‑100 (MQ=255)
    gCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTc                                                                                               >  2:1346234/1‑100 (MQ=255)
    gCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTc                                                                                               <  2:16910/100‑1 (MQ=255)
      cTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTg                                                                                             <  2:1901178/100‑1 (MQ=255)
        aCGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGTATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGCg                                                                                           <  1:1764042/100‑1 (MQ=255)
        aCGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGACGGTTATCCTTCCGc                                                                                            >  2:767819/1‑99 (MQ=255)
        aCGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGCg                                                                                           <  2:2030899/100‑1 (MQ=255)
        aCGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGCg                                                                                           <  1:1424874/100‑1 (MQ=255)
          gcgGATCATCGATGTAGGTCGGATAAGGCCCTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGTAACTGTAACGCAGAAGGATAACCTTCTGCGtt                                                                                         <  1:2265868/100‑1 (MQ=39)
          gcgGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGTAACTGTAACGCAGAAGGATAACCTTCTGCGtt                                                                                         <  1:2266124/100‑1 (MQ=39)
          gcgGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGCGtt                                                                                         <  2:74055/100‑1 (MQ=255)
            ggATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGtttt                                                                                       <  1:1129383/100‑1 (MQ=255)
            ggATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGtttt                                                                                       >  2:149929/1‑100 (MQ=255)
             gATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGCGttac                                                                                       <  1:2022843/99‑3 (MQ=255)
              atcatcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTg                                                                                     <  2:1532708/100‑1 (MQ=255)
              atcatcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTg                                                                                     <  2:909793/100‑1 (MQ=255)
              atcatcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTg                                                                                     <  1:179824/100‑1 (MQ=255)
                catcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGtttt                                                                                       >  2:2286253/1‑96 (MQ=255)
                 atcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCGGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGttt                                                                                  <  2:1834308/100‑1 (MQ=255)
                 atcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGttt                                                                                  <  1:1083379/100‑1 (MQ=255)
                 atcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGttt                                                                                  <  1:1173023/100‑1 (MQ=255)
                  tcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTa                                                                                 <  2:1513468/100‑1 (MQ=255)
                  tcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTa                                                                                 >  1:862699/1‑100 (MQ=255)
                  tcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTa                                                                                 >  2:417996/1‑100 (MQ=255)
                  tcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTa                                                                                 <  1:69300/100‑1 (MQ=255)
                   cGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTaa                                                                                >  1:1312638/1‑100 (MQ=255)
                   cGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTaa                                                                                >  1:2025377/1‑100 (MQ=255)
                      tGTCGGTCGGATAAGGCTATCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTa                                                                             <  2:587533/100‑1 (MQ=255)
                       gtatctCGGATCAGCCACTCGCCGCCCCCGGCAAGATATATCGCACGGTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTa                                                                             <  2:1954425/94‑1 (MQ=255)
                       gTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAg                                                                            >  2:30659/1‑100 (MQ=255)
                        tAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGCTTTTGTTTAATTCGc                                                                           <  2:702716/100‑1 (MQ=255)
                           gTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGt                                                                        >  1:791602/1‑100 (MQ=255)
                              ggATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAg                                                                     >  1:2338825/1‑100 (MQ=255)
                                    ggCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGc                                                               >  1:550959/1‑100 (MQ=255)
                                    ggCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGc                                                               >  1:1350834/1‑100 (MQ=255)
                                    ggCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGc                                                               >  2:1741068/1‑100 (MQ=255)
                                      cACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCaa                                                             >  1:1163967/1‑100 (MQ=255)
                                      cACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCaa                                                             >  2:2378046/1‑100 (MQ=255)
                                        cTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAAc                                                           >  1:1140924/1‑100 (MQ=255)
                                        cTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAAc                                                           <  1:1640650/100‑1 (MQ=255)
                                          cgccgcATCCGGCAAGATAAGTCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACt                                                          >  2:1459543/1‑99 (MQ=255)
                                          cgccgcATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTg                                                         >  1:1096493/1‑100 (MQ=255)
                                          cgccgcATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTg                                                         >  1:426301/1‑100 (MQ=255)
                                           gccgcATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGt                                                        >  2:195883/1‑100 (MQ=255)
                                             cgcATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGttt                                                      >  1:296591/1‑100 (MQ=255)
                                              gcATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTc                                                     <  2:1724380/100‑1 (MQ=255)
                                               cATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCa                                                    <  1:417998/100‑1 (MQ=255)
                                                  ccGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgt                                                 >  1:2061432/1‑100 (MQ=255)
                                                    ggCAAGATAAATTGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgttg                                               <  2:513774/100‑1 (MQ=255)
                                                    ggCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgttg                                               >  1:577451/1‑100 (MQ=255)
                                                      cAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgttgtt                                             <  2:1601296/100‑1 (MQ=255)
                                                      cAAGATAAATCGCACGTTGTCACCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgt                                                 <  1:1066242/96‑1 (MQ=255)
                                                       aaGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGttt                                            >  2:803214/1‑100 (MQ=255)
                                                          aTAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCt                                         >  1:1753986/1‑100 (MQ=255)
                                                             aaTCGCACGTTGTCCGCAACTGTAACGCAGAATGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGc                                      <  2:242627/100‑1 (MQ=255)
                                                              aTCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCt                                     <  2:296589/100‑1 (MQ=255)
                                                              aTCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCt                                     <  1:1436428/100‑1 (MQ=255)
                                                                 gCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGAc                                  >  1:1306454/1‑100 (MQ=255)
                                                                 gCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGAc                                  <  1:1133867/100‑1 (MQ=255)
                                                                 gCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGAc                                  >  1:325976/1‑100 (MQ=255)
                                                                 gCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGAc                                  <  1:343193/100‑1 (MQ=255)
                                                                 gCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTGAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGAc                                  <  1:2258612/100‑1 (MQ=255)
                                                                  cACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACa                                 >  1:1524745/1‑100 (MQ=255)
                                                                  cACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACa                                 <  2:1228302/100‑1 (MQ=255)
                                                                    cGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATa                               <  2:325974/100‑1 (MQ=255)
                                                                       tctccgCAACTTTAACGTAAAAGGTTATACTTCTGCGTTTTTGTTTAAGTTTCTTTTAGCAACGCAAACTGTTTCAGGTTGTTTTCGGGCTGACATAAac                            <  2:554809/95‑1 (MQ=255)
                                                                        gTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAaca                           <  1:1956210/100‑1 (MQ=255)
                                                                         tCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAacac                          <  2:192128/100‑1 (MQ=255)
                                                                          cAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACAcc                         >  1:1190365/1‑100 (MQ=255)
                                                                          cAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACAcc                         >  1:1238203/1‑100 (MQ=255)
                                                                          cAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACAcc                         >  1:1137105/1‑100 (MQ=255)
                                                                          attaaacaaaaaCGCAGAAGGATAACCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACAcc                         >  2:1979395/11‑100 (MQ=37)
                                                                           aGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCa                        >  2:1813658/1‑100 (MQ=255)
                                                                            gCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAg                       >  2:5666/1‑100 (MQ=255)
                                                                            gCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAg                       >  2:1850382/1‑100 (MQ=255)
                                                                             aaacaaaaaCGCAGAAGGATAACCTTCTGCGGTTTTGTTTAATTGGCTGACAGCAACGCAGACAGTTACAGGTTGGTTTCTGGCTGACATAAaca                           >  2:1033027/8‑95 (MQ=16)
                                                                                 aaaaaCGCAGAAGGATAACCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCaata                  >  2:1084790/4‑100 (MQ=255)
                                                                                  gTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCaataa                 >  1:1665914/1‑100 (MQ=255)
                                                                                  gTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCaataa                 >  1:1889534/1‑100 (MQ=255)
                                                                                   aaaCGCAGAAGGATAACCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCaataat                >  1:1868869/2‑100 (MQ=255)
                                                                                    aaCGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATg               <  1:710460/100‑1 (MQ=255)
                                                                                    aaCGCAGAAGGATAACCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATg               <  1:1979412/100‑1 (MQ=255)
                                                                                     aCGCAGAAGGATAACCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGc              >  1:623250/1‑100 (MQ=255)
                                                                                     aCGCAGAAGGATAACCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGc              >  1:624715/1‑100 (MQ=255)
                                                                                     aCGCAGAAGGATAACCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGc              >  1:1870177/1‑100 (MQ=255)
                                                                                      cGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGcc             <  2:1998737/100‑1 (MQ=255)
                                                                                      cGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGcc             >  2:560328/1‑100 (MQ=255)
                                                                                      cGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGcc             >  1:857284/1‑100 (MQ=255)
                                                                                        cAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAg           <  1:2146946/100‑1 (MQ=255)
                                                                                        cAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAg           <  1:195884/100‑1 (MQ=255)
                                                                                        cAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAg           <  1:1584870/100‑1 (MQ=255)
                                                                                              gTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCcgc     <  2:1501600/100‑1 (MQ=255)
                                                                                              gTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCcgc     <  2:1414333/100‑1 (MQ=255)
                                                                                               ttATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCcgcg    <  1:30660/100‑1 (MQ=255)
                                                                                                 aTCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGAc  <  1:1267473/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
TCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGAC  >  NC_000913/1736002‑1736198

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