Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 503,650 A→G G349G (GGT→GGC fsr ← putative fosmidomycin efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913503,6500AG99.4% 566.7 / NA 153G349G (GGT→GGCfsrputative fosmidomycin efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (71/81);  total (72/81)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 153 total aligned reads displayed.

GAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGG  >  NC_000913/503503‑503797
                                                                                                                                                   |                                                                                                                                                   
gAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCg                                                                                                                                                  >  2:1185939/1‑151 (MQ=255)
 aaTATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGa                                                                                                                                                 <  2:591436/151‑1 (MQ=255)
   tatGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCg                                                                                                                                                  >  1:1467189/1‑148 (MQ=255)
          aaCATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTcc                                                                                                                                          <  2:1620498/149‑1 (MQ=255)
          aaCATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTg                                                                                                                                        <  2:1673108/151‑1 (MQ=255)
            cATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTgg                                                                                                                                       >  1:633546/1‑150 (MQ=255)
            cATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTgg                                                                                                                                       >  2:1175460/1‑150 (MQ=255)
             aTCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGTCGGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGTGTCCAGAAACCATGCCGATACGTCCTGGa                                                                                                                                      <  1:1628248/150‑1 (MQ=255)
                 ccAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAagcag                                                                                                                                 <  2:1071329/151‑1 (MQ=255)
                 ccAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAACCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAagc                                                                                                                                   >  1:321828/1‑149 (MQ=255)
                  cAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAagc                                                                                                                                   >  4:108615/1‑148 (MQ=255)
                   aaTAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCt                                                                                                                               <  1:467405/151‑1 (MQ=255)
                      aGTGGCAGGAAAGCACAGATTTTATAGCCTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCt                                                                                                                            <  4:532488/151‑1 (MQ=255)
                      aGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCt                                                                                                                            <  3:533971/151‑1 (MQ=255)
                      aGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAACCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTcc                                                                                                                             >  2:411224/1‑150 (MQ=255)
                         ggCAGTAAAGAACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCACAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGCCCTGGAAGCAGCTCCTGa                                                                                                                          <  1:1388008/150‑1 (MQ=255)
                         ggCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGa                                                                                                                          <  1:448390/150‑1 (MQ=255)
                          gCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGc                                                                                                                        <  4:437871/151‑1 (MQ=255)
                           cAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTg                                                                                                                           >  2:400922/1‑147 (MQ=255)
                           cAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCg                                                                                                                       >  2:1105542/1‑151 (MQ=255)
                              gAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGt                                                                                                                      >  4:311900/1‑149 (MQ=255)
                               aaaGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGa                                                                                                                   <  3:251283/151‑1 (MQ=255)
                                aaGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATAc                                                                                                                                              >  1:935036/1‑123 (MQ=255)
                                  gCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGAcc                                                                                                                 >  1:131646/1‑150 (MQ=255)
                                   cacaGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGAcc                                                                                                                 >  2:769893/1‑149 (MQ=255)
                                   cacaGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGAcc                                                                                                                 >  2:432206/1‑149 (MQ=255)
                                   cacaGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAg                                                                                                               <  1:281445/151‑1 (MQ=255)
                                   cacaGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAg                                                                                                               <  1:1487950/151‑1 (MQ=255)
                                            tataGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGc                                                                                                         >  1:293095/1‑148 (MQ=255)
                                              taGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCagag                                                                                                     <  1:876251/150‑1 (MQ=255)
                                              taGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCag                                                                                                       >  2:892985/1‑148 (MQ=255)
                                                 aCTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATg                                                                                                 >  2:956914/1‑151 (MQ=255)
                                                   tatcccGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATg                                                                                                 <  1:523792/144‑1 (MQ=255)
                                                      cTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGa                                                                                             <  1:47778/150‑1 (MQ=255)
                                                       tCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGa                                                                                             >  2:705299/1‑149 (MQ=255)
                                                        cGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGaa                                                                                            <  1:1503434/149‑1 (MQ=255)
                                                        cGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGc                                                                                          <  2:202337/151‑1 (MQ=255)
                                                        cGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGc                                                                                          <  2:1065507/151‑1 (MQ=255)
                                                          aTGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCg                                                                                         <  1:39169/150‑1 (MQ=255)
                                                           tGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCg                                                                                         >  1:768605/1‑149 (MQ=255)
                                                            gCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAgg                                                                                      <  4:225502/151‑1 (MQ=255)
                                                            gCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAg                                                                                       <  3:139429/150‑1 (MQ=255)
                                                             cTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAgg                                                                                      >  4:394476/1‑150 (MQ=255)
                                                             cTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAgg                                                                                      >  3:447343/1‑150 (MQ=255)
                                                              tGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAgg                                                                                      <  1:595838/149‑1 (MQ=255)
                                                               ggTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGAt                                                                                    <  3:513341/150‑1 (MQ=255)
                                                                gtgtGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATaa                                                                                  <  2:962046/151‑1 (MQ=255)
                                                                  gtgATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAAt                                                                                >  2:1448041/1‑151 (MQ=255)
                                                                    gATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATcc                                                                              <  1:1205706/151‑1 (MQ=255)
                                                                      tCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCa                                                                             >  1:1572968/1‑150 (MQ=255)
                                                                      tCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCa                                                                             >  2:532633/1‑150 (MQ=255)
                                                                       cGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCaat                                                                           >  1:856872/1‑151 (MQ=255)
                                                                       cGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCaa                                                                            <  4:41944/150‑1 (MQ=255)
                                                                          cGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCaataa                                                                         >  1:949091/1‑150 (MQ=255)
                                                                            aTAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATc                                                                       >  1:1105361/1‑150 (MQ=255)
                                                                            aTAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCa                                                                      <  1:526410/151‑1 (MQ=255)
                                                                             tAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCaataataa                                                                      >  1:658825/1‑148 (MQ=255)
                                                                             tAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCaataat                                                                        <  2:244797/148‑1 (MQ=255)
                                                                               aGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCg                                                                   >  3:70639/1‑151 (MQ=255)
                                                                                  cccAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGtt                                                                 >  1:1101826/1‑150 (MQ=255)
                                                                                    cAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTaaa                                                              <  2:433563/151‑1 (MQ=255)
                                                                                              gCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTcc                                                    >  1:443800/1‑151 (MQ=255)
                                                                                              gCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTcc                                                    <  3:188569/151‑1 (MQ=255)
                                                                                                tCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACATCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTcc                                                    >  2:1015578/1‑149 (MQ=255)
                                                                                                 cccAGACCTCCCATGCCAAAAGCAAATCCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCGCCGTTAAAACCCCCGTCCa                                                   <  1:1325294/149‑1 (MQ=255)
                                                                                                 cccAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAg                                                  >  1:1406805/1‑150 (MQ=255)
                                                                                                 cccAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGt                                                 <  3:311918/151‑1 (MQ=255)
                                                                                                     gACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAg                                             <  2:1091736/151‑1 (MQ=255)
                                                                                                      aCCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAgg                                            >  2:172422/1‑151 (MQ=255)
                                                                                                        cTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCt                                          >  1:1502519/1‑151 (MQ=255)
                                                                                                          cccATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCt                                          <  1:694102/149‑1 (MQ=255)
                                                                                                          cccATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTg                                         <  4:52877/150‑1 (MQ=255)
                                                                                                            cATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCg                                      <  1:307139/151‑1 (MQ=255)
                                                                                                             aTGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGt                                     <  2:1105496/151‑1 (MQ=255)
                                                                                                                ccAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAgg                                  >  2:341408/1‑151 (MQ=255)
                                                                                                                ccAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAgg                                  >  3:147583/1‑151 (MQ=255)
                                                                                                                 cAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAggg                                 <  2:583361/151‑1 (MQ=255)
                                                                                                                  aaaaGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGt                                <  2:1506047/151‑1 (MQ=255)
                                                                                                                   aaaGCAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTa                               <  1:1736154/151‑1 (MQ=255)
                                                                                                                       cAAAACCGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAAt                           >  2:1432241/1‑151 (MQ=255)
                                                                                                                             cGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTGAAATCAGCg                      >  1:1127821/1‑150 (MQ=255)
                                                                                                                             cGAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGt                     >  2:1518191/1‑151 (MQ=255)
                                                                                                                              gAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGt                     >  1:1258225/1‑150 (MQ=255)
                                                                                                                              gAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGt                     >  3:98606/1‑150 (MQ=255)
                                                                                                                              gAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGt                     >  4:342453/1‑150 (MQ=255)
                                                                                                                              gAAAAAGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTa                    <  3:51692/151‑1 (MQ=255)
                                                                                                                                aaaaGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTaa                   >  2:1266173/1‑150 (MQ=255)
                                                                                                                                 aaaGAGTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAAc                 >  2:928212/1‑151 (MQ=255)
                                                                                                                                      gTCCAGAAGCCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGGAGAGAATGCCGAAGCGAGGATAAATCCAAGAATCGCCGTTAAACCCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGgc              <  2:894883/149‑1 (MQ=255)
                                                                                                                                      gTCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGgcgc            <  4:409137/151‑1 (MQ=255)
                                                                                                                                       tCCAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGgcgc            >  2:487545/1‑150 (MQ=255)
                                                                                                                                        ccAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCaa          <  2:906421/151‑1 (MQ=255)
                                                                                                                                        ccAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCaa          <  4:33881/151‑1 (MQ=255)
                                                                                                                                        ccAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCaa          <  2:1238077/151‑1 (MQ=255)
                                                                                                                                        ccAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCa           <  3:585610/150‑1 (MQ=255)
                                                                                                                                         cAGAAACCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAAc         >  2:1533346/1‑151 (MQ=255)
                                                                                                                                            aaaCCATGCCGTTACGTCCTGGAAGCAGCTCCTGTGCGTAGGCCAGAATGTCAGAGAATGCCGAAGCGAGGATAAATCCAATAAGCACAGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTTAACTGCGCAACGc       >  1:399693/1‑150 (MQ=255)
                                                                                                                                             aaCCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCg     >  2:712421/1‑151 (MQ=255)
                                                                                                                                              aCCATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGa    <  1:1220946/151‑1 (MQ=255)
                                                                                                                                               ccATGCCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAg   <  4:579906/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGG  >  NC_000913/503503‑503797

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