Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,123,625 C→G 5.4% C299S (TGC→TCC) ‡ cpsB ← mannose‑1‑phosphate guanyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,123,6250CG 5.4% 318.4 / 8.8 93C299S (TGC→TCC) ‡cpsBmannose‑1‑phosphate guanyltransferase
Reads supporting (aligned to +/- strand):  ref base C (45/43);  new base G (2/3);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

ACATCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGA  >  NC_000913/2123478‑2123772
                                                                                                                                                   |                                                                                                                                                   
acaTCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGc                                                                                                                                                     <  2:152238/148‑1 (MQ=255)
acaTCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGc                                                                                                                                                     >  2:125149/1‑148 (MQ=255)
 caTCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCCCCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCa                                                                                                                                                    >  3:5690/1‑148 (MQ=255)
  aTCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCaaa                                                                                                                                                  >  3:278428/1‑149 (MQ=255)
     cTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGt                                                                                                                                               >  2:105081/1‑149 (MQ=255)
        tACCGCGTTATGGTCGGCAATGAGCACCGCATCTTTTGTCTTCACCACTACCATATCTTTCTCCCCGACGTTTTTGACCAGGCCAGATTCAGCTTACACATATCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGtt                                                                                                                                              <  1:104988/147‑1 (MQ=255)
        tACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGc                                                                                                                                            <  6:235433/149‑1 (MQ=255)
        tACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGc                                                                                                                                            <  4:5690/149‑1 (MQ=255)
          ccGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGcc                                                                                                                                           <  4:278428/148‑1 (MQ=255)
                 aCGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTcggcg                                                                                                                                    <  1:69447/148‑1 (MQ=255)
                 aCGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTcggcg                                                                                                                                    <  1:243719/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                                                                                               <  1:125149/149‑1 (MQ=255)
                                                     cacTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGc                                                                                               <  7:184320/149‑1 (MQ=255)
                                                     cacTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGc                                                                                               >  8:17977/1‑149 (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                                        <  1:410799/149‑1 (MQ=255)
                                                                                                            tAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTc                                        >  2:410799/1‑149 (MQ=255)
                                                                                                              gCTGTTTTCTGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTcc                                       <  6:184942/148‑1 (MQ=255)
                                                                                                              gCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCa                                      <  3:20776/149‑1 (MQ=255)
                                                                                                              gCTGGTTTCAGTTTTTTGATTAATCACATCGCCGTGTCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGTCACCACAACAGCATCTGCCGTACGTTcc                                       <  5:160011/148‑1 (MQ=255)
                                                                                                               cTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCAt                                     <  8:279746/149‑1 (MQ=255)
                                                                                                               cTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCAt                                     >  7:273534/1‑149 (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             <  3:278708/149‑1 (MQ=255)
                                                                                                                                       acaTCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTc             >  5:291700/1‑149 (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)
                                                                                                                                                   |                                                                                                                                                   
ACATCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGA  >  NC_000913/2123478‑2123772

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: