Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,123,634 C→G 5.2% G296A (GGC→GCC) ‡ cpsB ← mannose‑1‑phosphate guanyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,123,6340CG 5.2% 306.9 / 6.1 91G296A (GGC→GCC) ‡cpsBmannose‑1‑phosphate guanyltransferase
Reads supporting (aligned to +/- strand):  ref base C (43/42);  new base G (2/3);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

TACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGCTT  >  NC_000913/2123486‑2123781
                                                                                                                                                    |                                                                                                                                                   
tACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGc                                                                                                                                                     <  4:5690/149‑1 (MQ=255)
tACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGc                                                                                                                                                     <  6:235433/149‑1 (MQ=255)
  ccGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGcc                                                                                                                                                    <  4:278428/148‑1 (MQ=255)
         aCGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTcggcg                                                                                                                                             <  1:243719/148‑1 (MQ=255)
         aCGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTcggcg                                                                                                                                             <  1:69447/148‑1 (MQ=255)
          cGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCTCCCTGGCAAAAGTTGACCTCGGCGgt                                                                                                                                           >  6:160011/1‑149 (MQ=255)
          cGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACACCGCCGAGGGCAACGTTTGCCAcggcgat                                                                                                                                           <  7:137518/149‑3 (MQ=255)
            gTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGGTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGgtgt                                                                                                                                         >  4:180542/1‑149 (MQ=255)
             tCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGtgtg                                                                                                                                        <  5:135888/149‑1 (MQ=255)
                  aaTCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGgcgc                                                                                                                                   <  5:58337/149‑1 (MQ=255)
                   aTCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCt                                                                                                                                  >  4:278708/1‑149 (MQ=255)
                   aTCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACACCGCCGAGGGCAACGTTTGCCACGGCGAtgtgattaat                                                                                                                                  <  1:338001/149‑7 (MQ=37)
                    tCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTcg                                                                                                                                                >  4:262793/1‑134 (MQ=255)
                    tCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTcg                                                                                                                                                <  3:262793/134‑1 (MQ=255)
                      aGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGa                                                                                                                                <  2:361865/148‑1 (MQ=255)
                       gCAGCGCATCTTTGGTCTGCACCACTACCAGATCTTTCCCCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGCCAAACGTGGCCCTCGGCGGGGTGGGCGCTGAtc                                                                                                                              <  2:205930/149‑1 (MQ=255)
                         aCCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGAtctc                                                                                                                            >  4:213894/1‑149 (MQ=255)
                          ccGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTcc                                                                                                                           >  5:623/1‑149 (MQ=255)
                             cATCTTTGGTCTGCACCACTACCAGATCTTTCACCACGACGGTGGAGACCAGGACAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCa                                                                                                                         >  4:263605/1‑148 (MQ=255)
                                 tttGGTCTGCACCACTACCAGATCTTTCCCCCCCACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTATTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCACTCGGCGGTGTggg                                                                                                                                      >  1:25757/1‑131 (MQ=255)
                                   tGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATaa                                                                                                                      <  8:179076/145‑1 (MQ=255)
                                    ggTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGaa                                                                                                                  <  1:192723/148‑1 (MQ=255)
                                     gTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGa                                                                                                                <  2:354005/149‑1 (MQ=255)
                                      tCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGAc                                                                                                               <  3:213894/149‑1 (MQ=255)
                                         gCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAg                                                                                                            >  3:315760/1‑149 (MQ=255)
                                         gCACCACTACCAGATATTTCACCCCGACGGTCGTGACCACGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATCAACACCAg                                                                                                            >  4:200160/1‑149 (MQ=255)
                                          caccacTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCa                                                                                                             <  6:623/147‑1 (MQ=255)
                                          caccacTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAg                                                                                                            <  1:105081/148‑1 (MQ=255)
                                           accacTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAg                                                                                                            >  2:331902/1‑147 (MQ=255)
                                           accacTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGa                                                                                                          >  5:292117/1‑149 (MQ=255)
                                           accacTACCAGATCTTTAACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCCGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAg                                                                                                            >  5:297835/1‑147 (MQ=255)
                                            ccacTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAg                                                                                                         >  7:279746/1‑149 (MQ=255)
                                             cacTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGc                                                                                                        >  8:17977/1‑149 (MQ=255)
                                             cacTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGc                                                                                                        <  1:125149/149‑1 (MQ=255)
                                             cacTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGc                                                                                                        <  7:184320/149‑1 (MQ=255)
                                              acTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGc                                                                                                        >  8:155765/1‑148 (MQ=255)
                                               cTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCa                                                                                                      <  4:196895/149‑1 (MQ=255)
                                                  ccAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCaa                                                                                                     <  1:275054/147‑1 (MQ=255)
                                                   cAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACAt                                                                                                  >  5:173205/1‑149 (MQ=255)
                                                   cAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACACCGCCGAGGGCAACGTTTGCCACGGCGatgt                                                                                                                                         >  8:312348/1‑110 (MQ=37)
                                                   cAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACACCGCCGAGGGCAACGTTTGCCACGGCGatgt                                                                                                                                         <  7:312348/110‑1 (MQ=37)
                                                    aGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATc                                                                                                 >  2:318146/1‑149 (MQ=255)
                                                     gATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCATCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTTGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCAATAATGAATACAAGCAGCCAACATCg                                                                                                >  5:286210/1‑149 (MQ=255)
                                                     gATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCg                                                                                                >  5:200997/1‑149 (MQ=255)
                                                           tCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTc                                                                                             <  1:331902/146‑1 (MQ=255)
                                                           tCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAg                                                                                          >  4:20776/1‑149 (MQ=255)
                                                              ccccGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGccc                                                                                       >  3:119108/1‑149 (MQ=255)
                                                                ccGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCg                                                                                      >  3:16/1‑148 (MQ=255)
                                                                      gtggtgACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACCTCGCTCCAGCCCGcatccat                                                                               <  6:111593/149‑1 (MQ=255)
                                                                       tggtgACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGcatccatc                                                                              <  6:292117/149‑1 (MQ=255)
                                                                       tggtgACCAGGCCAGATTCAGAATACACATAGCTGTTTTCAGTTTTTTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGcatccatc                                                                              <  6:297835/149‑1 (MQ=255)
                                                                         gtgACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCgg                                                                            >  5:97874/1‑149 (MQ=255)
                                                                            aCCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGca                                                                          <  6:173205/148‑1 (MQ=255)
                                                                                              tACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGCCCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCg                                                        >  5:385729/1‑148 (MQ=255)
                                                                                                    tAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTc                                                 >  2:410799/1‑149 (MQ=255)
                                                                                                    tAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTc                                                 <  1:410799/149‑1 (MQ=255)
                                                                                                      gCTGTTTTCTGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTcc                                                <  6:184942/148‑1 (MQ=255)
                                                                                                      gCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCa                                               <  3:20776/149‑1 (MQ=255)
                                                                                                      gCTGGTTTCAGTTTTTTGATTAATCACATCGCCGTGTCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGTCACCACAACAGCATCTGCCGTACGTTcc                                                <  5:160011/148‑1 (MQ=255)
                                                                                                       cTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCAt                                              >  7:273534/1‑149 (MQ=255)
                                                                                                       cTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCAt                                              <  8:279746/149‑1 (MQ=255)
                                                                                                        tGTTTTCAGTTTTGTGATTAATGACATCGCCGTGGCAAACGTTGCCCTCGGGGGTGTGGGCGGTGATCTCCCATTATGAAGACCAGGAGCGAACATCGCTCCAGCCCGCATCCATCGGCACCACAGCAGCATCTGCCGTACGTTACATg                                             >  8:137785/1‑149 (MQ=255)
                                                                                                         gTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGa                                            >  7:330536/1‑149 (MQ=255)
                                                                                                            ttCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGAcc                                          <  6:312176/148‑1 (MQ=255)
                                                                                                            ttCAGTTTTGGGATTAATCACATCGCCTTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGCTCTCCCATAATGAAGACCAGGAGCCCACATCGCTCCAGCCCGCATCCCTCGGCCCCACAACCGCAGCTGCCGTACGTTCCCTTACcg                                         >  4:79750/1‑149 (MQ=255)
                                                                                                               aGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGt                                      <  3:180542/149‑1 (MQ=255)
                                                                                                                gTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTc                                                                                             <  6:240050/93‑1 (MQ=255)
                                                                                                                gTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTc                                                                                             >  5:240050/1‑93 (MQ=255)
                                                                                                                gTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCAt                                                              >  2:272396/1‑124 (MQ=255)
                                                                                                                 ttttGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCACCATCGCTCCAGCCCGCATCCATCGGCACGCCAACAGCATCTGCCGTACGTTCCATGACCGCGTaa                                    >  6:305872/1‑149 (MQ=255)
                                                                                                                 ttttGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTaa                                    >  5:349771/1‑149 (MQ=255)
                                                                                                                   ttGTGATTAATCACATCGCCGTGGCAAAGGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAAt                                   <  1:318146/148‑1 (MQ=255)
                                                                                                                    tgtgATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATc                                  <  7:155765/148‑1 (MQ=255)
                                                                                                                      tgATTAATCACATCTCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACGACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCa                                <  4:16/148‑1 (MQ=255)
                                                                                                                      tgATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCAc                               >  8:386596/1‑149 (MQ=255)
                                                                                                                       gATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCCCAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCAcc                              >  5:9676/1‑149 (MQ=255)
                                                                                                                        aTTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCAcc                              <  6:349771/148‑1 (MQ=255)
                                                                                                                        aTTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCAcc                              >  5:371396/1‑148 (MQ=255)
                                                                                                                             tCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGccc                                                                                       >  2:119235/1‑86 (MQ=255)
                                                                                                                               acaTCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTc                      >  5:291700/1‑149 (MQ=255)
                                                                                                                               acaTCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTc                      <  3:278708/149‑1 (MQ=255)
                                                                                                                                 aTCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCg                    >  6:80131/1‑149 (MQ=255)
                                                                                                                                 aTCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCg                    <  7:386596/149‑1 (MQ=255)
                                                                                                                                 aTCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCg                    >  7:79676/1‑149 (MQ=255)
                                                                                                                                  tCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGCCTCTTCCgg                   >  6:390439/1‑149 (MQ=255)
                                                                                                                                   cGCCGAGGGCAACGTTTGCCACGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCggg                  >  5:320897/1‑149 (MQ=255)
                                                                                                                                      cGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGAGCAGGAGCAAACATCGCTCCAGCCCTCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAg               <  3:200160/149‑1 (MQ=255)
                                                                                                                                          gCAAACGTTGCCCTCGGCGTTGTGGGAGCTGATATCACATAATTAAGACCAGGAGCAAACATCGCTCCAGACCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCga           <  6:385729/149‑1 (MQ=255)
                                                                                                                                          gCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATTAAGACCAGGAGACAACATCGCTCCAGCCCGCATCCATCGGAAACTCATCAGCATCTGCCGTTCGTTGCATTACCGCGTTATCCACCGACTCTTCCGGGCAGGCg            <  6:286210/148‑1 (MQ=255)
                                                                                                                                            aaaCGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCAc                               <  4:202902/127‑1 (MQ=255)
                                                                                                                                            aaaCGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCAc                               >  3:202902/1‑127 (MQ=255)
                                                                                                                                               cGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAAc      <  6:291700/149‑1 (MQ=255)
                                                                                                                                                 ttGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATGCACCGACTCTTCCCGGCAGGCGAGAAACGc    <  6:200997/149‑1 (MQ=255)
                                                                                                                                                 ttGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGc    >  3:191089/1‑149 (MQ=255)
                                                                                                                                                    cccTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGctt  <  3:263605/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGCTT  >  NC_000913/2123486‑2123781

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