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

GCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAG  >  NC_000913/3608549‑3608838
                                                                                                                                                   |                                                                                                                                                 
gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGGGCGCGATCGGTGCCCTCGCCCCTCT‑gg                                                                                                                                                >  1:441988/1‑149 (MQ=255)
  tGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑gggg                                                                                                                                              <  1:546658/149‑1 (MQ=255)
    ttAGGGATGACCGAGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCATACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCG‑GGGgag                                                                                                                                            >  2:792884/1‑149 (MQ=255)
      aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGggag                                                                                                                                          >  1:280285/1‑149 (MQ=255)
      aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGagag                                                                                                                                          >  1:699596/1‑149 (MQ=255)
      aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCCCT‑GGGGagag                                                                                                                                          >  1:699518/1‑149 (MQ=255)
      aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCCGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGagag                                                                                                                                          <  1:660970/149‑1 (MQ=255)
      aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCCGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGagag                                                                                                                                          <  2:1070971/149‑1 (MQ=255)
       gggAGGGCCGGGGTGGGGCGGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCCAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAgg                                                                                                                                         <  2:122885/149‑1 (MQ=255)
           tGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACCGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTa                                                                                                                                     >  1:575347/1‑149 (MQ=255)
           tGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTa                                                                                                                                     <  2:639056/149‑1 (MQ=255)
             aCCGGGTTTGTGCTGGCGCTGCTGGCAGAGTCGGGGGCGGCGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGGTCAGGAGATAAGGCAAACCGCTC‑GCAACATTGGGCGCGATCGGTCCCCGCGCCCCTCT‑GGGGAGCGGGTTAgg                                                                                                                                   <  2:139813/149‑1 (MQ=255)
               cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCCCT‑GGGGAGGGGGGTAGGGt                                                                                                                                 >  1:1151173/1‑149 (MQ=255)
                     gtgGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCCGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAgggg                                                                                                                           >  2:54815/1‑149 (MQ=255)
                       ggctggcCGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑aa                                                                                                                        <  2:1016727/149‑1 (MQ=255)
                        gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTGGCGCAGGAGCTAAGGCAAACCGATC‑GCAACGTTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGGGGGTTCGGGTGAGGGG‑aac                                                                                                                       >  2:1002952/1‑148 (MQ=39)
                         ctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAACCCGATC‑GCACCCTTGAGCGCGATCGCTCTCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑aaaa                                                                                                                      <  1:1229562/149‑1 (MQ=39)
                         ctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑aaaa                                                                                                                      <  2:1158882/149‑1 (MQ=255)
                          tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAg                                                                                                                     <  2:1082819/149‑1 (MQ=255)
                           ggcAGTGCTGGCGGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAgg                                                                                                                    <  1:770855/149‑1 (MQ=255)
                           ggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAgg                                                                                                                    <  1:1251733/149‑1 (MQ=255)
                            gcAGTGGTGGCAGATCCGGGGGCGACGGTGCTGGTGACGGCGAATGCGTTAATATTGTGGCGCAGGAGGTAAGGCGAACCGATC‑GCAACATGGAGCGGGATCGGTCCCCGCGCCCCGCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAggc                                                                                                                   <  2:1228572/149‑1 (MQ=37)
                                tGCTGGCAGATACGGGGGCGACGGTGTTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGGGGGTTAGGGTGAGGGG‑AAAAGGCGGCa                                                                                                               >  1:1010480/1‑149 (MQ=255)
                                tGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGGGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCa                                                                                                               >  1:231583/1‑149 (MQ=255)
                                tGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCa                                                                                                               >  1:1209326/1‑149 (MQ=255)
                                   tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGAGAGGGG‑GAAAGGCGGCATCg                                                                                                            >  1:812617/1‑149 (MQ=255)
                                   tGGCAGATACGGGGGCGACGGTGCTGGTACCAGCGAATGCGTTAAGATTGTTGCGCAGGAGATGAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCG‑GGGGAGAGGGTTAGGGTGAGGGG‑AGAAGGCGGCATCg                                                                                                            >  1:844533/1‑149 (MQ=39)
                                    ggCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATACGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGa                                                                                                           >  1:865781/1‑149 (MQ=255)
                                     gCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATACGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCC‑GGGGAGAGGGTTAGGGCGAGGGG‑AAACGGGGGCATCGaa                                                                                                          >  2:1126898/1‑149 (MQ=39)
                                     gCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGaa                                                                                                          >  2:305686/1‑149 (MQ=255)
                                              ggggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGCTAGGGTGAGGGG‑GAAAGGCGGCATCGAAGCCAATCAg                                                                                                 >  1:462597/1‑149 (MQ=255)
                                                 gCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCGGccc                                                                                              >  1:599840/1‑149 (MQ=255)
                                                    aCGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GGAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAAGCAGccccca                                                                                           >  2:88156/1‑147 (MQ=255)
                                                     cGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAGGCCAATCAGCCCCTAt                                                                                          >  2:425973/1‑149 (MQ=255)
                                                          cTGGTGACAGAGAAGGCGTTAAGATTGTTGCGCAGGAGATAAGGAAAACCGATC‑GAAACATTGAGCGCGATCGGCCCCCTCGCCCCCCG‑GGGGGGGGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAcc                                                                                     <  1:918337/149‑1 (MQ=255)
                                                            ggTCACCGCGAATGCGGGAAGATTGGGGCGCGGGAGAGAAGGCAAACCGTGC‑GCAACATTGGGCGTGATGTTACCCCTCGCCCCTTG‑GGGGAGAGTGTTAGGGTGAGGGG‑AGAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                   <  2:355513/149‑1 (MQ=255)
                                                                     gAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAAAATTGAGCGCGATCGGTGACCTCGCACCTCT‑GGGGAGAGGGTGAGGGTGTGGGG‑AAAAGGCGGCATACAAGCCAATCAACCTCTATTAACCCCCTTTAcata                                                                          >  1:502063/1‑146 (MQ=255)
                                                                      aaTGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAAt                                                                         >  1:1178064/1‑149 (MQ=255)
                                                                      aaTGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTAACAACCGCCTTTCCGaaa                                                                         >  2:127267/1‑148 (MQ=255)
                                                                       aTGCGTTAAGATTGTTGCTCAGGAGATAAGGCAGACCGATC‑GCAACATTGAGCGCGCTCGATAGCCTCGCCCCTCT‑GGGTGGTGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATc                                                                        >  1:950424/1‑149 (MQ=255)
                                                                         gCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCaa                                                                      >  1:770648/1‑149 (MQ=255)
                                                                           gTTCAGATTGTTGCGCAGGAGATAAGGCAAACCGCTC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGGTAGGGTGAGGGG‑AAAAGGCGACATCAAAGCCCGTCAGCCCCCATAAACAGCCTTTAGGAATAAAAt                                                                    >  2:326438/1‑149 (MQ=37)
                                                                           gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAAt                                                                    >  1:230437/1‑149 (MQ=255)
                                                                           gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTGTACGAATCAAAt                                                                    >  2:870374/1‑149 (MQ=255)
                                                                           gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTCTCCGAATCaaaa                                                                    >  2:122019/1‑148 (MQ=255)
                                                                               aGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACg                                                                >  1:1242345/1‑149 (MQ=255)
                                                                                   tgttgCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATaa                                                            >  2:796048/1‑149 (MQ=255)
                                                                                     tggcgaAGGAGA‑AACGGAAACCGGTCGGCAACTTGGGGCGCGCTCGGCCCCCTCCCCCCTCG‑GGGGAGAGGGTTGGGGGGGGGGA‑GAAATGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                                                          <  1:864436/147‑1 (MQ=255)
                                                                                       gcgcGGGGGATAGGGCAACCGGTGC‑GCAACGTTGGGGGGGATCGGCCCCCTCGCCCCTCG‑GGGGGGGGGGTTGGGGTGGGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                        <  2:699596/149‑1 (MQ=255)
                                                                                       gcgcAGGAGAGACGGCACACCGCTC‑GCAATCTGGCGCGCGAGCGTGCCCCTCGCCTCTCG‑GGGGAGAGGGTTAGGGGGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                        <  2:930065/149‑1 (MQ=255)
                                                                                       cggagggAGTTAGGGCAACCCAATC‑GAAAGATGGGCGGGGATCGGTCCCCCCGCCCCTCG‑GGGGGGGGGGTTAGGGCGGGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                        <  2:280285/144‑1 (MQ=255)
                                                                                        cgcAGGATATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGGGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCACTATCAACCGCCTTTACGAATCAAAAAACAATAAGGCAg                                                       >  1:921367/1‑149 (MQ=38)
                                                                                        cgcAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAg                                                       >  1:904648/1‑149 (MQ=255)
                                                                                                aTAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCg                                               >  2:783118/1‑149 (MQ=255)
                                                                                                aTAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCG‑GGGGAGAGGGCTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTGACGAATCAAATAACGATAAGGCAGTCCATCCg                                               >  2:993146/1‑149 (MQ=255)
                                                                                                 tagggCAACCCGATC‑GCAACATTGGGCGGGATCGGGCCCCCGGCCCCGCT‑GGGGAGGGGGTTAGGGGCGGGGG‑ACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                                              <  2:152147/146‑1 (MQ=255)
                                                                                                  aaGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTAGCAACCGCCTTTACGAATCACATAACGATAAGGCAGTCCATCCGtc                                             >  2:239518/1‑149 (MQ=255)
                                                                                                   aGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCGCCTCGCCCCTCT‑GGGGAGAGGGTAAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCACCTATCAACCGCCGTGACGCATCAAATAACGAGAAGGCAGTCCCACCGtct                                            >  2:781519/1‑149 (MQ=38)
                                                                                                     gCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGGGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGTTACGGCAGTCCATCCGTctct                                          >  1:295352/1‑149 (MQ=255)
                                                                                                     gCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                          >  1:414750/1‑149 (MQ=255)
                                                                                                        agccGCTC‑GCAACATGGGGCCCGAGCGGCCCCCCCGCCCCAGT‑GGGGGGAGGGCTAGGTTGGGGGG‑AAAAGGCGGCATCGAAGCCAACCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                       <  2:462597/147‑1 (MQ=255)
                                                                                                        accgactC‑GAAACATGGGGCGGGACGGGCCCCGCCGCCCCATG‑GGGGGAAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                       <  1:307203/143‑1 (MQ=255)
                                                                                                        acccGATC‑CAAACATCGGCCCCGATCGGGCCCC‑CCCCCCTCC‑GGGGGGGGGGGTAGGGGGGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                       <  1:331156/147‑1 (MQ=255)
                                                                                                        aaCCGATC‑GCAACATTGAGCGCGATCGGTCCCCAAAAATCTCT‑TTTGAGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGAAGTCCATCCGTCTCtttc                                       >  1:633778/1‑148 (MQ=255)
                                                                                                         aCCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGCGGGGTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGCCCATCCGTCTCCTTAg                                      >  2:248406/1‑149 (MQ=255)
                                                                                                                gaaaaCAGTGGGCGCGAT‑GGCCCCCCCCCCCCTTTGGGGGAGGGGGTGGGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCa                                <  1:54815/146‑1 (MQ=255)
                                                                                                                  cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAGGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGt                              >  2:1145125/1‑149 (MQ=255)
                                                                                                                       tggcgcgcGACCGGCCCCCCCGCCCCTCT‑GGGGAGAGGGGTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCATTCCATCCGTCTCTTTAGCAACCAGTTCgtg                         <  1:681575/145‑1 (MQ=255)
                                                                                                                          ggcgcGATCGGTCCCCTCGCCCCGCG‑GGGGAGAGGGTGAGGGTGAGGGG‑AAAGGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTc                      <  2:1026877/148‑1 (MQ=255)
                                                                                                                          gccgcgTTGGGTCCCGTCCCCCCTCT‑GGGGAGGGGGTTGGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTc                      <  2:1072073/147‑1 (MQ=255)
                                                                                                                           gcggtaTCGGCCCCTCCGCCCCCTG‑GGGGAGGGGGTTAGGGGGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTATGAATCATATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                     <  2:844533/144‑1 (MQ=255)
                                                                                                                              ggagcGGTCCCCCGGCCCCTAG‑GGGGAGGGGGTGAGGGGGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAAGAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATa                  <  1:545568/145‑1 (MQ=255)
                                                                                                                               gATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                 >  1:108647/1‑149 (MQ=255)
                                                                                                                                          cccccccTCG‑GGGGGGAGGGTTGGGGTGAGGGA‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc      <  1:122019/146‑1 (MQ=255)
                                                                                                                                             ccccTGG‑GGGGGGGGGGTTGGGGGGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa   <  1:148687/149‑1 (MQ=255)
                                                                                                                                              cccTCT‑GGGGAGGGGGTTAGGGGGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg  <  2:1128959/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATC‑GCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCT‑GGGGAGAGGGTTAGGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAG  >  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: ‑