Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 551,678 T→C 5.0% D306G (GAC→GGC)  purK ← N5‑carboxyaminoimidazole ribonucleotide synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913551,6780TC 5.0% 249.9 / 2.6 79D306G (GAC→GGC) purKN5‑carboxyaminoimidazole ribonucleotide synthase
Reads supporting (aligned to +/- strand):  ref base T (35/40);  new base C (2/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

CCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATCGGTAATCGCCCGCAGATGCAGCTCAAACT  >  NC_000913/551531‑551822
                                                                                                                                                   |                                                                                                                                                
ccGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGGCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTc                                                                                                                                                 <  5:37820/149‑1 (MQ=255)
 cGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCg                                                                                                                                                >  7:62777/1‑149 (MQ=255)
    aCTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGTCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTAc                                                                                                                                             <  6:212273/149‑1 (MQ=255)
      ttACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTAc                                                                                                                                             >  7:70332/1‑147 (MQ=255)
      ttACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTAc                                                                                                                                             >  4:114667/1‑147 (MQ=255)
          tctGCGCCTAAATCACGCCGCTGGCATATTCCTGCTGCAGAAGCTTGATTAAGGCTTCCTGAGTCGCAGTCAGACGATATGTGTAGCTTTCGGTCAAATTCAGATGCCCCACTTTACGCAACGGACTTACTTCTTTGTCGTACCa                                                                                                                                           <  8:395651/145‑1 (MQ=255)
          tctGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCa                                                                                                                                           >  7:395651/1‑145 (MQ=255)
          tctGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTg                                                                                                                                        >  1:437804/1‑148 (MQ=255)
            tGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGCCTTCTTTGTCGTACCAGTg                                                                                                                                        >  5:12829/1‑146 (MQ=255)
            tGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGc                                                                                                                                       <  6:175491/147‑1 (MQ=255)
              cgcCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCa                                                                                                                                      <  7:321801/146‑1 (MQ=255)
              cgcCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGAt                                                                                                                                   <  7:101284/149‑1 (MQ=255)
               gcCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCCCTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATg                                                                                                                                  >  3:99262/1‑149 (MQ=255)
               gcCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATg                                                                                                                                  >  8:319701/1‑149 (MQ=255)
                cccAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGc                                                                                                                                 >  2:348755/1‑149 (MQ=255)
                  cAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCAc                                                                                                                               <  7:319701/149‑1 (MQ=255)
                    aaTCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCAc                                                                                                                               >  6:258064/1‑147 (MQ=255)
                         cgccgcTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGc                                                                                                                        <  8:409767/149‑1 (MQ=255)
                          gccgcTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCa                                                                                                                       <  4:21209/149‑1 (MQ=255)
                          gccgcTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACACCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCa                                                                                                                       <  7:17259/149‑1 (MQ=255)
                                 gCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTc                                                                                                                  <  2:409594/147‑1 (MQ=255)
                                  cATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGc                                                                                                               <  1:348755/149‑1 (MQ=255)
                                  cATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGc                                                                                                               <  8:70332/149‑1 (MQ=255)
                                         cggcggCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCa                                                                                                         <  6:89724/148‑1 (MQ=255)
                                             ggcAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAAtt                                                                                                    <  1:365525/149‑1 (MQ=255)
                                                  cagcGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACAt                                                                                               >  8:163656/1‑149 (MQ=255)
                                                     cGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGc                                                                                            <  5:258064/149‑1 (MQ=255)
                                                        gATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGc                                                                                                                                 <  3:75361/109‑1 (MQ=255)
                                                        gATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGc                                                                                                                                 >  4:75361/1‑109 (MQ=255)
                                                        gATTAAGGCTTACAGCGTCTAAGTCAGTAGCGATGTGTCGCTGTATGTCAAATTCAGATTCCCCACTTTTCGCCGAGGTCGGAGTTCTTTGTCGTAGCAGATAAGATGCACGAGCGGCAGTTTCAGCCATTCATAATTCACATCGCTAc                                                                                         <  8:111581/149‑1 (MQ=255)
                                                         aTTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTAc                                                                                         >  1:248636/1‑148 (MQ=255)
                                                             aGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGt                                                                                                                                         <  6:377530/96‑1 (MQ=255)
                                                             aGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGt                                                                                                                                         >  5:377530/1‑96 (MQ=255)
                                                             aGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATc                                                                                    <  6:12829/149‑1 (MQ=255)
                                                              ggCTTCCGGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAAt                                                                                     >  8:447081/1‑147 (MQ=255)
                                                              ggCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAAt                                                                                     >  4:26944/1‑147 (MQ=255)
                                                               gCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAAt                                                                                     >  4:208656/1‑146 (MQ=255)
                                                                 ttCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGAt                                                                                >  1:457917/1‑149 (MQ=255)
                                                                    cAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGa                                                                             <  3:26944/149‑1 (MQ=255)
                                                                     aGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGat                                                                            <  7:147215/149‑1 (MQ=255)
                                                                     aGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGat                                                                            >  4:278928/1‑149 (MQ=255)
                                                                      gCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGatc                                                                           <  7:163656/149‑1 (MQ=255)
                                                                         tCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGCTGCACCAGCGGCAGTTTCCGCCAGTCATAATTCACATCGCTACCAATCAGATTGatcatc                                                                        >  2:295844/1‑149 (MQ=255)
                                                                         tCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGatcatc                                                                        <  3:278928/149‑1 (MQ=255)
                                                                         tCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGatcatc                                                                        <  7:458441/149‑1 (MQ=255)
                                                                           gCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCAc                                                                      <  3:117957/149‑1 (MQ=255)
                                                                           gCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGACATAATTCACATCGATACCAATCAGATTGATCATCAc                                                                      >  1:402876/1‑149 (MQ=255)
                                                                                 aGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAAtt                                                                                                    >  4:244969/1‑113 (MQ=255)
                                                                                 aGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAAtt                                                                                                    <  3:244969/113‑1 (MQ=255)
                                                                                        aTGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGattatt                                                          >  5:391358/1‑148 (MQ=255)
                                                                                           tgtCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTcac                                                       <  8:62777/148‑1 (MQ=255)
                                                                                              cGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACt                                                   >  4:49216/1‑149 (MQ=255)
                                                                                                     ggTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCg                                            <  6:391358/149‑1 (MQ=255)
                                                                                                     ggTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGGCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGc                                             <  4:99262/148‑1 (MQ=255)
                                                                                                      gTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCgg                                           <  2:457917/149‑1 (MQ=255)
                                                                                                         aaaTTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCAACGACGGATTATTCACCACTCGTTGCGGTaa                                        >  4:220237/1‑149 (MQ=255)
                                                                                                          aaTTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAAc                                       <  2:437804/149‑1 (MQ=255)
                                                                                                            ttCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACgg                                     <  5:61032/149‑1 (MQ=255)
                                                                                                                     cccACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGatt                                                             <  6:327408/116‑1 (MQ=255)
                                                                                                                     cccACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGt                                          <  5:108063/135‑1 (MQ=255)
                                                                                                                     cccACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGt                                          >  6:108063/1‑135 (MQ=255)
                                                                                                                     cccACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACAGCGCTACCAATCAGATTGATCATCACCGACGGatt                                                             >  5:327408/1‑116 (MQ=255)
                                                                                                                             aCGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATCGGTAATCGcc                     <  8:390302/148‑1 (MQ=255)
                                                                                                                                ccccGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATCGGTAATCGCCCGCa                 <  2:248636/149‑1 (MQ=255)
                                                                                                                                 cccGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGc                                    <  3:27824/129‑1 (MQ=255)
                                                                                                                                 cccGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGc                                    >  4:27824/1‑129 (MQ=255)
                                                                                                                                 cccGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACAACTGGTTGCGGTAACGGCAGATCGGTAATCGCCCg                   >  2:455034/1‑146 (MQ=255)
                                                                                                                                   cggacggaCTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATATTTCACCACTGGTTGCGGTAACGGCAGATCGGTAATCGCCCGCAGAt              >  7:421008/1‑149 (MQ=255)
                                                                                                                                     gacggacTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATCGGTAATCGCCCGCAGATGc            <  7:447081/149‑1 (MQ=255)
                                                                                                                                       cggacTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATCGGTAATCGCCCGCAGATGCAg          >  3:19515/1‑149 (MQ=255)
                                                                                                                                          acTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATCGGTAATCGCCCGCAGATGCAGCt        >  3:173989/1‑148 (MQ=255)
                                                                                                                                          acTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATCGGTAATCGCCCGCAGATGCAGCTc       <  1:455034/149‑1 (MQ=255)
                                                                                                                                            ttcttTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGCCGGATTATTCACCACTGGTTGCGGTAACGGCAGATCGGTAATCGCCCGCAGATGCAGCTCa      >  5:26079/1‑148 (MQ=255)
                                                                                                                                            ttcttTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATCGGTAATCGCCCGCAGATGCAGCTCa      >  1:45010/1‑148 (MQ=255)
                                                                                                                                               tttGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATCGGTAATCGCCCGCAGATGCAGCTCAAACt  >  6:344801/1‑149 (MQ=255)
                                                                                                                                               tttGCCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATCGGTa                           >  1:429759/1‑124 (MQ=255)
                                                                                                                                               tttGCCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATCGGTa                           >  1:436257/1‑124 (MQ=255)
                                                                                                                                               tttGCCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATCGGTa                           <  2:436257/124‑1 (MQ=255)
                                                                                                                                               tttGCCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATCGGTa                           <  2:429759/124‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCAGCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTCAGCCAGTCATAATTCACATCGCTACCAATCAGATTGATCATCACCGACGGATTATTCACCACTGGTTGCGGTAACGGCAGATCGGTAATCGCCCGCAGATGCAGCTCAAACT  >  NC_000913/551531‑551822

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑