Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,093,969 T→A 5.6% Y168N (TAC→AAC)  hisB → histidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,093,9690TA 5.6% 381.3 / 10.5 107Y168N (TAC→AAC) hisBhistidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (47/54);  new base A (4/2);  total (51/56)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

Only 100 of 107 total aligned reads displayed.

CGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCAT  >  NC_000913/2093823‑2094114
                                                                                                                                                  |                                                                                                                                                 
cGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTAc                                                                                                                                                 <  3:325446/149‑1 (MQ=255)
  aTGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCGTTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGc                                                                                                                                               <  3:671403/149‑1 (MQ=255)
     gATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCa                                                                                                                                            >  3:526180/1‑149 (MQ=255)
       tCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACg                                                                                                                                          >  3:778655/1‑149 (MQ=255)
         gcgcTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTa                                                                                                                                        <  4:199296/149‑1 (MQ=255)
             tAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGc                                                                                                                                    >  1:538399/1‑149 (MQ=255)
               aCAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGt                                                                                                                                  <  4:805766/149‑1 (MQ=255)
                cAGTTATGTGATTGGCGATCGCGCGACCTACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTa                                                                                                                                 <  4:911152/149‑1 (MQ=255)
                 aGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTaa                                                                                                                                >  3:467570/1‑149 (MQ=255)
                 aGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTaa                                                                                                                                <  4:432995/149‑1 (MQ=255)
                    tATGTGATTGGCGATCGCGCGACCGACATTCAACTTGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATAc                                                                                                                             <  1:84742/149‑1 (MQ=255)
                    tATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATAc                                                                                                                             >  3:673910/1‑149 (MQ=255)
                      tgtgATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCa                                                                                                                           >  1:834863/1‑149 (MQ=255)
                       gtgATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCaa                                                                                                                          >  2:233437/1‑149 (MQ=255)
                        tgATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCaaa                                                                                                                         <  4:526180/149‑1 (MQ=255)
                        tgATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCaaa                                                                                                                         <  1:34432/149‑1 (MQ=255)
                          aTTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAaga                                                                                                                       <  1:874004/149‑1 (MQ=255)
                          aTTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAaga                                                                                                                       <  4:846786/149‑1 (MQ=255)
                                gATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCa                                                                                                                 >  3:783310/1‑149 (MQ=255)
                                    gcgcgACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGAtt                                                                                                             <  2:623211/149‑1 (MQ=255)
                                    gcgcgACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGAtt                                                                                                             <  1:821304/149‑1 (MQ=255)
                                            gacATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCa                                                                                                     >  3:577798/1‑149 (MQ=255)
                                              cATTGAACTGGAGGAAAACATAGCTATTACTGGTTTACGCTACGACCGCGAAACCCTGATCTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGGAGATTGACGTTCAgg                                                                                                   <  4:886476/149‑1 (MQ=255)
                                               aTTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGgt                                                                                                  >  1:681397/1‑149 (MQ=255)
                                               aTTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGgt                                                                                                  >  2:575949/1‑149 (MQ=255)
                                                 tCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAaga                                                                                                                       >  4:633388/1‑126 (MQ=255)
                                                 tCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAaga                                                                                                                       <  3:633388/126‑1 (MQ=255)
                                                  cAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGtgtg                                                                                               <  2:153516/149‑1 (MQ=255)
                                                   aaCTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTgg                                                                                              >  3:764297/1‑149 (MQ=255)
                                                   aaCTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTgg                                                                                              >  4:195928/1‑149 (MQ=255)
                                                       ggCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGGGAAACCCTGAACTGGCAAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTgg                                                                                          <  4:723068/149‑1 (MQ=255)
                                                         cGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGAt                                                                                        <  1:102502/149‑1 (MQ=255)
                                                          ggAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATc                                                                                       <  2:776737/149‑1 (MQ=255)
                                                             aaaCATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTg                                                                                    <  2:745002/149‑1 (MQ=255)
                                                               aCATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGaa                                                                                  <  3:647314/149‑1 (MQ=255)
                                                               aCATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGaa                                                                                  <  3:195928/149‑1 (MQ=255)
                                                                cATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAg                                                                                 <  4:227820/149‑1 (MQ=255)
                                                                cATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAg                                                                                 <  1:838767/149‑1 (MQ=255)
                                                                  tGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGt                                                                               <  2:66863/149‑1 (MQ=255)
                                                                  tGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGt                                                                               <  4:741991/149‑1 (MQ=255)
                                                                   gggCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTg                                                                              >  2:477264/1‑149 (MQ=255)
                                                                    ggCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTgg                                                                             >  4:219780/1‑149 (MQ=255)
                                                                     gCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGgc                                                                            <  3:397927/149‑1 (MQ=255)
                                                                       aTTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGgcag                                                                          >  1:545455/1‑149 (MQ=255)
                                                                         tACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGgcagca                                                                        >  4:153703/1‑149 (MQ=255)
                                                                              gTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGAtt                                                                   >  3:874168/1‑149 (MQ=255)
                                                                              gTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGAtt                                                                   >  3:97007/1‑149 (MQ=255)
                                                                              gTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGAtt                                                                   >  4:837809/1‑149 (MQ=255)
                                                                                ttACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTaa                                                                 <  3:67417/149‑1 (MQ=255)
                                                                                ttACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTaa                                                                 <  1:529124/149‑1 (MQ=255)
                                                                                  aCGCTACGACCGCGAAACCCTGAACTGGCCAAAGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGGAAGATTAaca                                                               >  4:792182/1‑149 (MQ=255)
                                                                                     cTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCg                                                            >  4:119437/1‑149 (MQ=255)
                                                                                      tACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCgg                                                           >  4:546234/1‑149 (MQ=255)
                                                                                             gcgAAACCCTGAACTGACCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGc                                                    <  4:848370/149‑1 (MQ=255)
                                                                                                   cccTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCttt                                              >  1:163659/1‑149 (MQ=255)
                                                                                                    ccTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTg                                             >  2:174551/1‑149 (MQ=255)
                                                                                                     cTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGCAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGa                                            <  1:455696/149‑1 (MQ=255)
                                                                                                     cTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGa                                            >  3:482085/1‑149 (MQ=255)
                                                                                                     cTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGa                                            >  2:937414/1‑149 (MQ=255)
                                                                                                      tGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGAt                                           <  4:340575/149‑1 (MQ=255)
                                                                                                       gAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTgg                                                                                              <  1:886371/97‑1 (MQ=255)
                                                                                                       gAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTgg                                                                                              >  2:886371/1‑97 (MQ=255)
                                                                                                        aaCTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCa                                         >  3:685958/1‑149 (MQ=255)
                                                                                                            ggCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCCTATTTGATCATATg                                     <  3:68711/149‑1 (MQ=255)
                                                                                                             gCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGc                                    <  1:396823/149‑1 (MQ=255)
                                                                                                              ccAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCt                                   >  1:494691/1‑149 (MQ=255)
                                                                                                                aaTGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTgg                                 >  1:425404/1‑149 (MQ=255)
                                                                                                                aaTGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTgg                                 >  3:235633/1‑149 (MQ=255)
                                                                                                                 aTGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGa                                <  3:316036/149‑1 (MQ=255)
                                                                                                                   gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATc                              <  2:425404/149‑1 (MQ=255)
                                                                                                                      tGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGAGTGCCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGa                           <  1:294644/149‑1 (MQ=255)
                                                                                                                       ggCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGAt                          <  4:783310/149‑1 (MQ=255)
                                                                                                                        gCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATc                         <  1:174551/149‑1 (MQ=255)
                                                                                                                        gCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATc                         >  2:924537/1‑149 (MQ=255)
                                                                                                                         cGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCg                        <  2:613133/149‑1 (MQ=255)
                                                                                                                           aGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGa                                            <  3:332722/127‑1 (MQ=255)
                                                                                                                           aGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGa                                            >  4:332722/1‑127 (MQ=255)
                                                                                                                              aaCTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTAcc                   >  1:560995/1‑149 (MQ=255)
                                                                                                                               aCTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGcttctt                                               <  3:588273/120‑1 (MQ=255)
                                                                                                                               aCTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGcttctt                                               >  4:588273/1‑120 (MQ=255)
                                                                                                                                 tCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGTCTTCTTTGATc                                          <  3:812392/123‑1 (MQ=255)
                                                                                                                                 tCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATc                                          >  4:812392/1‑123 (MQ=255)
                                                                                                                                 tCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCAc                >  3:125660/1‑149 (MQ=255)
                                                                                                                                   aCCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCAcgg              >  2:163486/1‑149 (MQ=255)
                                                                                                                                    ccAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCAcggc             >  1:407670/1‑149 (MQ=255)
                                                                                                                                    ccAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCAcggc             >  1:477774/1‑149 (MQ=255)
                                                                                                                                     cAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCAcggcg            >  3:378672/1‑149 (MQ=255)
                                                                                                                                     cAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCAcggc             <  2:528735/148‑1 (MQ=255)
                                                                                                                                     cAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCAcggc             >  1:528735/1‑148 (MQ=255)
                                                                                                                                     actaCGTGAGCGTAACGGTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCAcggcg            >  2:683774/4‑149 (MQ=255)
                                                                                                                                       taCGTGAGCGTAACGGTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCttt                                              >  4:540254/2‑113 (MQ=255)
                                                                                                                                        aCGTGAGCGTAACGGTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGAt                               >  4:703729/1‑127 (MQ=255)
                                                                                                                                        aCGTGAGCGTAACGGTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGAt                               <  3:703729/127‑1 (MQ=255)
                                                                                                                                        aCGTGAGCGTAACGGTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGtt         >  2:58541/1‑149 (MQ=255)
                                                                                                                                           tGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTcc      <  2:758590/149‑1 (MQ=255)
                                                                                                                                           tGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTcc      >  2:525395/1‑149 (MQ=255)
                                                                                                                                           tGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTcc      <  2:908142/149‑1 (MQ=255)
                                                                                                                                               cGTTACGCTCACGTAGTGCGTAATACCAAAGAGATGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCAt  >  2:43240/1‑149 (MQ=255)
                                                                                                                                               cGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCAt  <  1:575949/149‑1 (MQ=255)
                                                                                                                                               cGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCAt  <  1:728453/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCAT  >  NC_000913/2093823‑2094114

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