Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,093,973 C→G 5.5% A169G (GCT→GGT)  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,9730CG 5.5% 391.4 / 11.0 110A169G (GCT→GGT) hisBhistidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (50/54);  new base G (4/2);  total (54/56)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

Only 100 of 111 total aligned reads displayed.

ATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGAAAT  >  NC_000913/2093825‑2094120
                                                                                                                                                    |                                                                                                                                                   
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)
                    tgtgATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCa                                                                                                                                 >  1:834863/1‑149 (MQ=255)
                      tgATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCaaa                                                                                                                               <  4:526180/149‑1 (MQ=255)
                      tgATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCaaa                                                                                                                               <  4:549222/149‑1 (MQ=255)
                      tgATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCaaa                                                                                                                               <  1:55876/149‑1 (MQ=255)
                      tgATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCaaa                                                                                                                               <  1:34432/149‑1 (MQ=255)
                        aTTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAaga                                                                                                                             <  4:846786/149‑1 (MQ=255)
                        aTTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAaga                                                                                                                             <  1:874004/149‑1 (MQ=255)
                              gATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCa                                                                                                                       >  3:783310/1‑149 (MQ=255)
                                  gcgcgACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGAtt                                                                                                                   <  1:821304/149‑1 (MQ=255)
                                  gcgcgACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGAtt                                                                                                                   <  2:623211/149‑1 (MQ=255)
                                          gacATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCa                                                                                                           >  3:577798/1‑149 (MQ=255)
                                          gacATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCa                                                                                                           <  2:834863/149‑1 (MQ=255)
                                            cATTGAACTGGAGGAAAACATAGCTATTACTGGTTTACGCTACGACCGCGAAACCCTGATCTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGGAGATTGACGTTCAgg                                                                                                         <  4:886476/149‑1 (MQ=255)
                                             aTTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGgt                                                                                                        >  2:575949/1‑149 (MQ=255)
                                             aTTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGgt                                                                                                        >  1:681397/1‑149 (MQ=255)
                                               tCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAaga                                                                                                                             <  3:633388/126‑1 (MQ=255)
                                               tCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAaga                                                                                                                             >  4:633388/1‑126 (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)
                                                                 gggCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTg                                                                                    >  2:396823/1‑149 (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                                                                         >  4:837809/1‑149 (MQ=255)
                                                                            gTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGAtt                                                                         >  3:97007/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                                                  >  2:937414/1‑149 (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)
                                                                                                              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)
                                                                                                                         aGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGa                                                  >  4:332722/1‑127 (MQ=255)
                                                                                                                         aGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGa                                                  <  3:332722/127‑1 (MQ=255)
                                                                                                                            aaCTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTAcc                         >  1:560995/1‑149 (MQ=255)
                                                                                                                             aCTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGcttctt                                                     >  4:588273/1‑120 (MQ=255)
                                                                                                                             aCTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGcttctt                                                     <  3:588273/120‑1 (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                                                    <  3:540254/112‑1 (MQ=255)
                                                                                                                                     taCGTGAGCGTAACGGTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCttt                                                    >  4:540254/2‑113 (MQ=255)
                                                                                                                                      aCGTGAGCGTAACGGTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGAt                                     <  3:703729/127‑1 (MQ=255)
                                                                                                                                      aCGTGAGCGTAACGGTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGAt                                     >  4:703729/1‑127 (MQ=255)
                                                                                                                                      aCGTGAGCGTAACGGTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGtt               >  2:58541/1‑149 (MQ=255)
                                                                                                                                         tGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTcc            <  2:758590/149‑1 (MQ=255)
                                                                                                                                         tGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTcc            <  2:908142/149‑1 (MQ=255)
                                                                                                                                         tGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTcc            >  2:525395/1‑149 (MQ=255)
                                                                                                                                            ccGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCa         <  3:837809/149‑1 (MQ=255)
                                                                                                                                             cGTTACGCTCACGTAGTGCGTAATACCAAAGAGATGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCAt        >  2:43240/1‑149 (MQ=255)
                                                                                                                                             cGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCAt        <  1:575949/149‑1 (MQ=255)
                                                                                                                                                  cGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGaaa   >  4:268963/1‑149 (MQ=255)
                                                                                                                                                  cGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGaaa   <  2:560995/149‑1 (MQ=255)
                                                                                                                                                   gCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGAAAt  >  3:657832/1‑149 (MQ=255)
                                                                                                                                                   gCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGAAAt  <  1:525395/149‑1 (MQ=255)
                                                                                                                                                   gCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGAAAt  >  3:609368/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGAAAT  >  NC_000913/2093825‑2094120

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: