Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,322,700 A→G 6.9% intergenic (‑195/+33) secG ← / ← glmM preprotein translocase membrane subunit/phosphoglucosamine mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,322,7000AG 6.9% 292.6 / 8.0 87intergenic (‑195/+33)secG/glmMpreprotein translocase membrane subunit/phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base A (40/41);  new base G (2/4);  total (42/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTAC  >  NC_000913/3322563‑3322839
                                                                                                                                         |                                                                                                                                           
aGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCtt                                                                                                                                  <  2:746248/149‑1 (MQ=255)
 gTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttt                                                                                                                                 <  1:1198694/149‑1 (MQ=255)
  tttttAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCtttt                                                                                                                                <  4:1167472/149‑1 (MQ=255)
  tttttAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCtttt                                                                                                                                <  3:12582/149‑1 (MQ=255)
     ttAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCttttt                                                                                                                               >  4:169522/1‑147 (MQ=255)
     ttAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCttttt                                                                                                                               <  3:169522/147‑1 (MQ=255)
      tAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTa                                                                                                                              >  4:1321005/1‑147 (MQ=255)
      tAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTa                                                                                                                              <  3:1321005/147‑1 (MQ=255)
        aTCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttctca                                                                                                                          <  3:264472/149‑6 (MQ=255)
         tCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCtttt                                                                                                                                <  1:439146/142‑1 (MQ=255)
         tCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCtttt                                                                                                                                >  2:439146/1‑142 (MQ=255)
           tttCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGtt                                                                                                                       >  4:1148099/1‑149 (MQ=255)
              cccACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTa                                                                                                                              <  4:24545/139‑1 (MQ=255)
              cccACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTa                                                                                                                              >  3:24545/1‑139 (MQ=255)
               ccACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCt                                                                                                                   >  2:606942/1‑149 (MQ=255)
                cACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTa                                                                                                                              <  4:152452/137‑1 (MQ=255)
                cACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTa                                                                                                                              >  3:152452/1‑137 (MQ=255)
                 aCTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTaa                                                                                                                 <  3:1354333/149‑1 (MQ=255)
                           gtgtGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAAc                                                                                                       >  4:938610/1‑149 (MQ=255)
                              tgAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTa                                                                                                                              >  3:517754/1‑123 (MQ=255)
                              tgAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTa                                                                                                                              <  4:517754/123‑1 (MQ=255)
                                  tACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCt                                                                                                                   <  1:1177310/130‑1 (MQ=255)
                                  tACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCt                                                                                                                   >  2:1177310/1‑130 (MQ=255)
                                  tACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCtttt                                                                                                >  3:344539/1‑149 (MQ=255)
                                  tACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCtttt                                                                                                >  2:771352/1‑149 (MQ=255)
                                    cTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTAc                                                                                              <  1:606942/149‑1 (MQ=255)
                                     tAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgcc                                                                                                                                       >  4:254014/1‑107 (MQ=255)
                                     tAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgcc                                                                                                                                       <  3:254014/107‑1 (MQ=255)
                                     tAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTAc                                                                                              <  3:942788/148‑1 (MQ=255)
                                     tAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTAc                                                                                              >  4:942788/1‑148 (MQ=255)
                                      aaCCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTg                                                                                            >  1:5804/1‑149 (MQ=255)
                                        ccAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGc                                                                                           <  4:670254/148‑1 (MQ=255)
                                        ccAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGc                                                                                           >  3:670254/1‑148 (MQ=255)
                                           aaGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCg                                                                                       >  3:1029355/1‑149 (MQ=255)
                                            aGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAGCGCTTTAAACGGCTTTTACTGCATCgg                                                                                      >  3:1213091/1‑149 (MQ=255)
                                             gCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttt                                                                                                                                 <  3:696761/105‑1 (MQ=255)
                                             gCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttt                                                                                                                                 >  4:696761/1‑105 (MQ=255)
                                                   gACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCg                                                                               >  2:1434777/1‑149 (MQ=255)
                                                     ccTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgt                                                                             <  1:771352/149‑1 (MQ=255)
                                                     ccTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgt                                                                             <  2:5804/149‑1 (MQ=255)
                                                       ttCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgtgt                                                                           <  4:1029355/149‑1 (MQ=255)
                                                            aaGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCaaa                                                                      >  2:523176/1‑149 (MQ=255)
                                                               ggCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGAt                                                                                  <  3:1016782/134‑1 (MQ=255)
                                                               ggCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGAt                                                                                  >  4:1016782/1‑134 (MQ=255)
                                                                 cAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAg                                                                 <  4:344539/149‑1 (MQ=255)
                                                                  aaTTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGt                                                                <  1:523176/149‑1 (MQ=255)
                                                                       tATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCt                                                           <  3:938610/149‑1 (MQ=255)
                                                                             aTTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAg                                                                                                                         >  4:1268590/1‑81 (MQ=255)
                                                                             aTTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAg                                                                                                                         <  3:1268590/81‑1 (MQ=255)
                                                                             aTTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCAc                                                             <  4:768847/141‑1 (MQ=255)
                                                                             aTTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCAc                                                             >  3:768847/1‑141 (MQ=255)
                                                                              ttGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAGCGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTc                                                    <  4:1213091/149‑1 (MQ=255)
                                                                                 tATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTc                                                 >  1:1186054/1‑149 (MQ=255)
                                                                                    cAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTc                                              <  1:1434777/149‑1 (MQ=255)
                                                                                     aaCTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCg                                             >  3:5719/1‑149 (MQ=255)
                                                                                                aaaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAAc                                                                                                                <  4:1414909/71‑1 (MQ=255)
                                                                                                aaaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAAc                                                                                                                >  3:1414909/1‑71 (MQ=255)
                                                                                                 aaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTaa                                                                                                                 <  4:961287/69‑1 (MQ=255)
                                                                                                 aaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTaa                                                                                                                 >  3:961287/1‑69 (MQ=255)
                                                                                                     aGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCAcgc                             <  3:1148099/149‑1 (MQ=255)
                                                                                                              gATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTaa                      <  1:1113232/147‑1 (MQ=255)
                                                                                                              gATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGGCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTaa                      >  2:1113232/1‑147 (MQ=255)
                                                                                                                       gCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGt                  <  3:230948/142‑1 (MQ=255)
                                                                                                                       gCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGt                  >  4:230948/1‑142 (MQ=255)
                                                                                                                        cTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGcc          <  4:5719/149‑1 (MQ=255)
                                                                                                                        cTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGcc          >  1:920304/1‑149 (MQ=255)
                                                                                                                         tGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCg         <  2:1186054/149‑1 (MQ=255)
                                                                                                                         tGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCg         >  3:736289/1‑149 (MQ=255)
                                                                                                                           agaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCAcgc                             <  2:1360315/127‑1 (MQ=255)
                                                                                                                           agaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCAcgc                             >  1:1360315/1‑127 (MQ=255)
                                                                                                                           agaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGa       >  2:716349/1‑149 (MQ=255)
                                                                                                                            gaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTcca                                      <  2:1412783/117‑1 (MQ=255)
                                                                                                                            gaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTcca                                      >  1:1412783/1‑117 (MQ=255)
                                                                                                                             taaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGa       <  3:759438/146‑1 (MQ=255)
                                                                                                                             taaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGa       >  4:759438/2‑147 (MQ=255)
                                                                                                                                aaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCaaa                                                                      >  4:1362274/1‑81 (MQ=255)
                                                                                                                                aaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCaaa                                                                      <  3:1362274/81‑1 (MQ=255)
                                                                                                                                aaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCg                                             <  1:526893/106‑1 (MQ=255)
                                                                                                                                aaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCg                                             >  2:526893/1‑106 (MQ=255)
                                                                                                                                aaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTAc  <  1:716349/149‑1 (MQ=255)
                                                                                                                                aaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTAc  >  3:1336586/1‑149 (MQ=255)
                                                                                                                                 aaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAAc                     <  4:380364/129‑1 (MQ=255)
                                                                                                                                 aaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAAc                     >  3:380364/1‑129 (MQ=255)
                                                                                                                                      gACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTc                                                               >  4:326338/1‑82 (MQ=255)
                                                                                                                                      gACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTc                                                               <  3:326338/82‑1 (MQ=255)
                                                                                                                                        cAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTAc  <  3:942552/141‑1 (MQ=255)
                                                                                                                                        cAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTAc  >  4:942552/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                           
AGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTAC  >  NC_000913/3322563‑3322839

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