Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,334,705 C→T 8.5% intergenic (+24/‑204) ispB → / → sfsB octaprenyl diphosphate synthase/malPQ operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,334,7050CT 8.5% 227.7 / 8.2 73intergenic (+24/‑204)ispB/sfsBoctaprenyl diphosphate synthase/malPQ operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (31/35);  new base T (4/2);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

GGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTA  >  NC_000913/3334557‑3334824
                                                                                                                                                    |                                                                                                                        
ggACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgc                                                                                                                          <  3:667007/149‑1 (MQ=255)
     cgTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGTAGcgcgcg                                                                                                                    <  3:1740532/149‑1 (MQ=255)
      gTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCgg                                                                                                                   >  3:475044/1‑149 (MQ=255)
        cAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGa                                                                                                                 <  3:662228/149‑1 (MQ=255)
        cAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGa                                                                                                                 <  1:1538754/149‑1 (MQ=255)
          gCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATg                                                                                                               <  4:475044/149‑1 (MQ=255)
            gTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGAt                                                                                                                >  4:1193986/1‑146 (MQ=255)
            gTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGAt                                                                                                                <  3:1193986/146‑1 (MQ=255)
                 gAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAg                                                                                                        >  2:642981/1‑149 (MQ=255)
                    gaagaaGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGcgcgcg                                                                                                                    <  2:1718174/134‑1 (MQ=255)
                    gaagaaGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGcgcgcg                                                                                                                    >  1:1718174/1‑134 (MQ=255)
                     aagaagCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAg                                                                                                    >  4:1211495/1‑149 (MQ=255)
                     aagaagCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAg                                                                                                    <  1:642981/149‑1 (MQ=255)
                      agaagCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCCCATCCCGCGCGCTGCGCGGGATgagg                                                                                                            >  4:1411751/1‑137 (MQ=255)
                      agaagCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCCCATCCCGCGCGCTGCGCGGGATgagg                                                                                                            <  3:1411751/140‑4 (MQ=255)
                       gaagCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCCCATCCCGCGCGCTGCGCGGGATgaggg                                                                                                           <  1:1812179/140‑5 (MQ=255)
                       gaagCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCCCATCCCGCGCGCTGCGCGGGATgaggg                                                                                                           >  2:1812179/1‑136 (MQ=255)
                         agCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTcc                                                                                                >  4:1153907/1‑149 (MQ=255)
                          gCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCa                                                                                               >  3:1894295/1‑149 (MQ=255)
                                 aaGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGt                                                                                                             <  1:1182711/128‑1 (MQ=255)
                                 aaGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGt                                                                                                             >  2:1182711/1‑128 (MQ=255)
                                         gCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATaaaa                                                                                          <  3:1802096/139‑1 (MQ=255)
                                         gCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACAGCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATaaaa                                                                                          >  4:1802096/1‑139 (MQ=255)
                                            gCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAg                                                                              <  1:670836/148‑1 (MQ=255)
                                            gCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAg                                                                              >  2:670836/1‑148 (MQ=255)
                                               ttACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGctct                                                                          >  1:748158/1‑149 (MQ=255)
                                               ttACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCCCATCCCGCGCGCTGCGCGGGATgagg                                                                                                            >  4:894932/1‑112 (MQ=255)
                                               ttACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCCCATCCCGCGCGCTGCGCGGGATgagg                                                                                                            <  3:894932/115‑4 (MQ=255)
                                                                 acCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCta                                                                         <  4:585351/132‑1 (MQ=255)
                                                                 acCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCta                                                                         >  3:585351/1‑132 (MQ=255)
                                                                  ccccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAgg                                                       <  4:1894295/149‑1 (MQ=255)
                                                                  ccccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAgg                                                       <  4:1119022/149‑1 (MQ=255)
                                                                   cccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATaa                                                                      >  2:1132582/1‑133 (MQ=255)
                                                                   cccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATaa                                                                      <  1:1132582/133‑1 (MQ=255)
                                                                    ccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                     <  3:1211495/149‑1 (MQ=255)
                                                                     cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCGCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAt                                                    <  3:1153907/149‑1 (MQ=255)
                                                                      ttGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTAAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtg                                                   >  1:549141/1‑149 (MQ=255)
                                                                        ggCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                 <  2:748158/149‑1 (MQ=255)
                                                                         gCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                        >  1:489486/1‑141 (MQ=255)
                                                                         gCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                        <  2:489486/141‑1 (MQ=255)
                                                                           gagaAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCa                                                                                                         >  1:166414/1‑90 (MQ=255)
                                                                           gagaAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCa                                                                                                         <  2:166414/90‑1 (MQ=255)
                                                                                gCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAaca                                                                                       <  1:111324/103‑1 (MQ=255)
                                                                                gCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAaca                                                                                       >  2:111324/1‑103 (MQ=255)
                                                                                  aCTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAt                                                                                  >  3:78389/1‑106 (MQ=255)
                                                                                  aCTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAt                                                                                  <  4:78389/106‑1 (MQ=255)
                                                                                      aTCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                        <  2:599733/128‑1 (MQ=255)
                                                                                      aTCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                        >  1:599733/1‑128 (MQ=255)
                                                                                       tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCt                                                                 >  3:805188/1‑118 (MQ=255)
                                                                                       tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCt                                                                 <  4:805188/118‑1 (MQ=255)
                                                                                       tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTc                                  <  2:549141/149‑1 (MQ=255)
                                                                                        cGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTc                                                                                                           <  1:1746740/75‑1 (MQ=255)
                                                                                        cGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTc                                                                                                           >  2:1746740/1‑75 (MQ=255)
                                                                                           ccTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATaaaa                                                                                          >  2:1374657/1‑89 (MQ=255)
                                                                                           ccTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATaaaa                                                                                          <  1:1374657/89‑1 (MQ=255)
                                                                                            cTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                     >  1:879186/1‑125 (MQ=255)
                                                                                            cTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                     <  2:879186/125‑1 (MQ=255)
                                                                                                   acaTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCtat                                                                        <  1:1727966/99‑1 (MQ=255)
                                                                                                   acaTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCtat                                                                        >  2:1727966/1‑99 (MQ=255)
                                                                                                     aTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTc                                                                                                           >  2:513929/1‑62 (MQ=255)
                                                                                                     aTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTc                                                                                                           <  1:513929/62‑1 (MQ=255)
                                                                                                     aTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                        <  1:1284903/113‑1 (MQ=255)
                                                                                                     aTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                        >  2:1284903/1‑113 (MQ=255)
                                                                                                            ttCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCgg                                                                                                                   >  3:1322544/1‑47 (MQ=255)
                                                                                                            ttCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCgg                                                                                                                   <  4:1322544/47‑1 (MQ=255)
                                                                                                                aCGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTcc                                                                                                <  2:283569/62‑1 (MQ=255)
                                                                                                                aCGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTcc                                                                                                >  1:283569/1‑62 (MQ=255)
                                                                                                                                    ccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                   <  3:1809084/103‑1 (MQ=255)
                                                                                                                                    ccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                   >  4:1809084/1‑103 (MQ=255)
                                                                                                                                    ccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTc                                  <  1:1076817/104‑1 (MQ=255)
                                                                                                                                    ccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTc                                  >  2:1076817/1‑104 (MQ=255)
                                                                                                                                     cctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAg            <  2:1698331/125‑1 (MQ=255)
                                                                                                                                     cctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAg            >  1:1698331/1‑125 (MQ=255)
                                                                                                                                          tCCCGCGCGCTGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtg                                                   >  2:1508349/1‑81 (MQ=255)
                                                                                                                                          tCCCGCGCGCGGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATACCACCTGAAAAATAGGGAtg                                                   <  1:1508349/81‑1 (MQ=255)
                                                                                                                                                gcAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTa  <  2:1092336/124‑1 (MQ=255)
                                                                                                                                                gcAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTa  >  1:1092336/1‑124 (MQ=255)
                                                                                                                                                    |                                                                                                                        
GGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTA  >  NC_000913/3334557‑3334824

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: