Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,695 T→G 15.0% intergenic (+47/+55) zntA → / ← tusA zinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,608,6950TG15.0% 136.4 / 5.1 76intergenic (+47/+55)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base T (40/19);  new base G (3/12);  total (44/32)
Fisher's exact test for biased strand distribution p-value = 1.15e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAG  >  NC_000913/3608549‑3608838
                                                                                                                                                   |                                                                                                                                                 
gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGGGCGCGATCGGTGCCCTCGCCCCTCTgg                                                                                                                                                >  1:441988/1‑149 (MQ=255)
  tGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgggg                                                                                                                                              <  1:546658/149‑1 (MQ=255)
    ttAGGGATGACCGAGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCATACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCGGGGgag                                                                                                                                            >  2:792884/1‑149 (MQ=255)
      aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGggag                                                                                                                                          >  1:280285/1‑149 (MQ=255)
      aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                          >  1:699596/1‑149 (MQ=255)
      aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCCCTGGGGagag                                                                                                                                          >  1:699518/1‑149 (MQ=255)
      aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCCGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                          <  1:660970/149‑1 (MQ=255)
      aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCCGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                          <  2:1070971/149‑1 (MQ=255)
       gggAGGGCCGGGGTGGGGCGGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCCAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                         <  2:122885/149‑1 (MQ=255)
           tGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACCGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                     >  1:575347/1‑149 (MQ=255)
           tGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                     <  2:639056/149‑1 (MQ=255)
             aCCGGGTTTGTGCTGGCGCTGCTGGCAGAGTCGGGGGCGGCGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGGTCAGGAGATAAGGCAAACCGCTCGCAACATTGGGCGCGATCGGTCCCCGCGCCCCTCTGGGGAGCGGGTTAgg                                                                                                                                   <  2:139813/149‑1 (MQ=255)
               cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCCCTGGGGAGGGGGGTAGGGt                                                                                                                                 >  1:1151173/1‑149 (MQ=255)
                     gtgGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCCGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgggg                                                                                                                           >  2:54815/1‑149 (MQ=255)
                       ggctggcCGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                        <  2:1016727/149‑1 (MQ=255)
                        gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTGGCGCAGGAGCTAAGGCAAACCGATCGCAACGTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTCGGGTGAGGGGaac                                                                                                                       >  2:1002952/1‑148 (MQ=39)
                         ctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAACCCGATCGCACCCTTGAGCGCGATCGCTCTCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                      <  1:1229562/149‑1 (MQ=39)
                         ctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                      <  2:1158882/149‑1 (MQ=255)
                          tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAg                                                                                                                     <  2:1082819/149‑1 (MQ=255)
                           ggcAGTGCTGGCGGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAgg                                                                                                                    <  1:770855/149‑1 (MQ=255)
                           ggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAgg                                                                                                                    <  1:1251733/149‑1 (MQ=255)
                            gcAGTGGTGGCAGATCCGGGGGCGACGGTGCTGGTGACGGCGAATGCGTTAATATTGTGGCGCAGGAGGTAAGGCGAACCGATCGCAACATGGAGCGGGATCGGTCCCCGCGCCCCGCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggc                                                                                                                   <  2:1228572/149‑1 (MQ=37)
                                tGCTGGCAGATACGGGGGCGACGGTGTTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAAGGCGGCa                                                                                                               >  1:1010480/1‑149 (MQ=255)
                                tGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGGGAGGGTTAGGGTGAGGGGAAAAGGCGGCa                                                                                                               >  1:231583/1‑149 (MQ=255)
                                tGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCa                                                                                                               >  1:1209326/1‑149 (MQ=255)
                                   tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGAGAGGGGGAAAGGCGGCATCg                                                                                                            >  1:812617/1‑149 (MQ=255)
                                   tGGCAGATACGGGGGCGACGGTGCTGGTACCAGCGAATGCGTTAAGATTGTTGCGCAGGAGATGAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCGGGGGAGAGGGTTAGGGTGAGGGGAGAAGGCGGCATCg                                                                                                            >  1:844533/1‑149 (MQ=39)
                                    ggCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATACGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                           >  1:865781/1‑149 (MQ=255)
                                     gCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATACGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCCGGGGAGAGGGTTAGGGCGAGGGGAAACGGGGGCATCGaa                                                                                                          >  2:1126898/1‑149 (MQ=39)
                                     gCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGaa                                                                                                          >  2:305686/1‑149 (MQ=255)
                                              ggggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGCTAGGGTGAGGGGGAAAGGCGGCATCGAAGCCAATCAg                                                                                                 >  1:462597/1‑149 (MQ=255)
                                                 gCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCGGccc                                                                                              >  1:599840/1‑149 (MQ=255)
                                                    aCGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGGAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAAGCAGccccca                                                                                           >  2:88156/1‑147 (MQ=255)
                                                     cGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAGGCCAATCAGCCCCTAt                                                                                          >  2:425973/1‑149 (MQ=255)
                                                          cTGGTGACAGAGAAGGCGTTAAGATTGTTGCGCAGGAGATAAGGAAAACCGATCGAAACATTGAGCGCGATCGGCCCCCTCGCCCCCCGGGGGGGGGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAcc                                                                                     <  1:918337/149‑1 (MQ=255)
                                                            ggTCACCGCGAATGCGGGAAGATTGGGGCGCGGGAGAGAAGGCAAACCGTGCGCAACATTGGGCGTGATGTTACCCCTCGCCCCTTGGGGGAGAGTGTTAGGGTGAGGGGAGAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                   <  2:355513/149‑1 (MQ=255)
                                                                     gAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAAAATTGAGCGCGATCGGTGACCTCGCACCTCTGGGGAGAGGGTGAGGGTGTGGGGAAAAGGCGGCATACAAGCCAATCAACCTCTATTAACCCCCTTTAcata                                                                          >  1:502063/1‑146 (MQ=255)
                                                                      aaTGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAAt                                                                         >  1:1178064/1‑149 (MQ=255)
                                                                      aaTGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTAACAACCGCCTTTCCGaaa                                                                         >  2:127267/1‑148 (MQ=255)
                                                                       aTGCGTTAAGATTGTTGCTCAGGAGATAAGGCAGACCGATCGCAACATTGAGCGCGCTCGATAGCCTCGCCCCTCTGGGTGGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATc                                                                        >  1:950424/1‑149 (MQ=255)
                                                                         gCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCaa                                                                      >  1:770648/1‑149 (MQ=255)
                                                                           gTTCAGATTGTTGCGCAGGAGATAAGGCAAACCGCTCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGGTAGGGTGAGGGGAAAAGGCGACATCAAAGCCCGTCAGCCCCCATAAACAGCCTTTAGGAATAAAAt                                                                    >  2:326438/1‑149 (MQ=37)
                                                                           gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAAt                                                                    >  1:230437/1‑149 (MQ=255)
                                                                           gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTGTACGAATCAAAt                                                                    >  2:870374/1‑149 (MQ=255)
                                                                           gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTCTCCGAATCaaaa                                                                    >  2:122019/1‑148 (MQ=255)
                                                                               aGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACg                                                                >  1:1242345/1‑149 (MQ=255)
                                                                                   tgttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATaa                                                            >  2:796048/1‑149 (MQ=255)
                                                                                     tggcgaAGGAGAAACGGAAACCGGTCGGCAACTTGGGGCGCGCTCGGCCCCCTCCCCCCTCGGGGGAGAGGGTTGGGGGGGGGGAGAAATGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                                                          <  1:864436/147‑1 (MQ=255)
                                                                                       gcgcGGGGGATAGGGCAACCGGTGCGCAACGTTGGGGGGGATCGGCCCCCTCGCCCCTCGGGGGGGGGGGTTGGGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                        <  2:699596/149‑1 (MQ=255)
                                                                                       gcgcAGGAGAGACGGCACACCGCTCGCAATCTGGCGCGCGAGCGTGCCCCTCGCCTCTCGGGGGAGAGGGTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                        <  2:930065/149‑1 (MQ=255)
                                                                                       cggagggAGTTAGGGCAACCCAATCGAAAGATGGGCGGGGATCGGTCCCCCCGCCCCTCGGGGGGGGGGGTTAGGGCGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                        <  2:280285/144‑1 (MQ=255)
                                                                                        cgcAGGATATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCACTATCAACCGCCTTTACGAATCAAAAAACAATAAGGCAg                                                       >  1:921367/1‑149 (MQ=38)
                                                                                        cgcAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAg                                                       >  1:904648/1‑149 (MQ=255)
                                                                                                aTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCg                                               >  2:783118/1‑149 (MQ=255)
                                                                                                aTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCGGGGGAGAGGGCTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTGACGAATCAAATAACGATAAGGCAGTCCATCCg                                               >  2:993146/1‑149 (MQ=255)
                                                                                                 tagggCAACCCGATCGCAACATTGGGCGGGATCGGGCCCCCGGCCCCGCTGGGGAGGGGGTTAGGGGCGGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                                              <  2:152147/146‑1 (MQ=255)
                                                                                                  aaGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTAGCAACCGCCTTTACGAATCACATAACGATAAGGCAGTCCATCCGtc                                             >  2:239518/1‑149 (MQ=255)
                                                                                                   aGGCAAACCGATCGCAACATTGAGCGCGATCGGTCGCCTCGCCCCTCTGGGGAGAGGGTAAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCACCTATCAACCGCCGTGACGCATCAAATAACGAGAAGGCAGTCCCACCGtct                                            >  2:781519/1‑149 (MQ=38)
                                                                                                     gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGGGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGTTACGGCAGTCCATCCGTctct                                          >  1:295352/1‑149 (MQ=255)
                                                                                                     gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                          >  1:414750/1‑149 (MQ=255)
                                                                                                        agccGCTCGCAACATGGGGCCCGAGCGGCCCCCCCGCCCCAGTGGGGGGAGGGCTAGGTTGGGGGGAAAAGGCGGCATCGAAGCCAACCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                       <  2:462597/147‑1 (MQ=255)
                                                                                                        accgactCGAAACATGGGGCGGGACGGGCCCCGCCGCCCCATGGGGGGAAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                       <  1:307203/143‑1 (MQ=255)
                                                                                                        acccGATCCAAACATCGGCCCCGATCGGGCCCCCCCCCCTCCGGGGGGGGGGGTAGGGGGGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                       <  1:331156/147‑1 (MQ=255)
                                                                                                        aaCCGATCGCAACATTGAGCGCGATCGGTCCCCAAAAATCTCTTTTGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGAAGTCCATCCGTCTCtttc                                       >  1:633778/1‑148 (MQ=255)
                                                                                                         aCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGCGGGGTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGCCCATCCGTCTCCTTAg                                      >  2:248406/1‑149 (MQ=255)
                                                                                                                gaaaaCAGTGGGCGCGATGGCCCCCCCCCCCCTTTGGGGGAGGGGGTGGGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCa                                <  1:54815/146‑1 (MQ=255)
                                                                                                                  cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAGGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGt                              >  2:1145125/1‑149 (MQ=255)
                                                                                                                       tggcgcgcGACCGGCCCCCCCGCCCCTCTGGGGAGAGGGGTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCATTCCATCCGTCTCTTTAGCAACCAGTTCgtg                         <  1:681575/145‑1 (MQ=255)
                                                                                                                          ggcgcGATCGGTCCCCTCGCCCCGCGGGGGAGAGGGTGAGGGTGAGGGGAAAGGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTc                      <  2:1026877/148‑1 (MQ=255)
                                                                                                                          gccgcgTTGGGTCCCGTCCCCCCTCTGGGGAGGGGGTTGGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTc                      <  2:1072073/147‑1 (MQ=255)
                                                                                                                           gcggtaTCGGCCCCTCCGCCCCCTGGGGGAGGGGGTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTATGAATCATATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                     <  2:844533/144‑1 (MQ=255)
                                                                                                                              ggagcGGTCCCCCGGCCCCTAGGGGGAGGGGGTGAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAAGAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATa                  <  1:545568/145‑1 (MQ=255)
                                                                                                                               gATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                 >  1:108647/1‑149 (MQ=255)
                                                                                                                                          cccccccTCGGGGGGGAGGGTTGGGGTGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc      <  1:122019/146‑1 (MQ=255)
                                                                                                                                             ccccTGGGGGGGGGGGGTTGGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa   <  1:148687/149‑1 (MQ=255)
                                                                                                                                              cccTCTGGGGAGGGGGTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg  <  2:1128959/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAG  >  NC_000913/3608549‑3608838

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