Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,902,494 G→T 5.5% G200V (GGA→GTA) ‡ ygcE → putative kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,902,4940GT 5.5% 374.3 / 9.3 109G200V (GGA→GTA) ‡ygcEputative kinase
Reads supporting (aligned to +/- strand):  ref base G (50/53);  new base T (4/2);  total (54/55)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01

Only 100 of 111 total aligned reads displayed.

ATTATAAGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTA  >  NC_000913/2902347‑2902641
                                                                                                                                                   |                                                                                                                                                   
aTTATAAGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAgg                                                                                                                                                     <  2:417254/148‑1 (MQ=255)
  tataAGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGAc                                                                                                                                                   >  2:249147/1‑148 (MQ=255)
   ataAGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTg                                                                                                                                                 >  3:909143/1‑149 (MQ=255)
    taAGAGCTGGGCATGGAGCGAAGACGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATTTGCAAATGCCTGGCACCGTCCTCGGTCATATCACACCACAAGCCGCTCTGGCGACACATTTCCCGGCAGGACTGc                                                                                                                                                >  3:599402/1‑149 (MQ=255)
     aaGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGcc                                                                                                                                               >  4:791073/1‑149 (MQ=255)
     aaGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGc                                                                                                                                                >  2:953724/1‑148 (MQ=255)
      agagCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGcc                                                                                                                                               >  1:638408/1‑148 (MQ=255)
         gCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgt                                                                                                                                            <  1:882686/148‑1 (MQ=255)
            gggCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaaat                                                                                                                                        <  4:792407/149‑6 (MQ=255)
            gggCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgttgtt                                                                                                                                        <  4:214545/149‑1 (MQ=255)
            gggCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgttgt                                                                                                                                         <  2:203275/148‑1 (MQ=255)
             ggCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgttgtt                                                                                                                                        <  2:412179/148‑1 (MQ=255)
              gCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaaatgt                                                                                                                                      <  4:475811/149‑8 (MQ=255)
                            tgctgcGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGAt                                                                                                                         >  1:702233/1‑148 (MQ=255)
                              ctgcGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATaa                                                                                                                       <  2:644090/148‑1 (MQ=255)
                                gcGGTAATGGATAAGTTTATTATCCCCCGTCATATGCTGTTTGATGAGCAAATGGCTGGCAGCTTCCTCGCACATATCACACCACATGCCGCACTGGCGGCACATTTCCCGGCAGGACTGCCGGTTGTTTTTACCACTAGTGAAAAAcc                                                                                                                    <  4:186614/149‑1 (MQ=255)
                                gcGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAAc                                                                                                                     >  2:471272/1‑148 (MQ=255)
                                 cGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAAcc                                                                                                                    <  2:638408/148‑1 (MQ=255)
                                 cGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCg                                                                                                                   <  3:587876/149‑1 (MQ=255)
                                           tAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGct                                                                                                          <  2:855774/148‑1 (MQ=255)
                                               tttAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTggg                                                                                                     >  4:185732/1‑149 (MQ=255)
                                               tttAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTggg                                                                                                     >  4:558307/1‑149 (MQ=255)
                                               tttAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTggg                                                                                                     <  3:791073/149‑1 (MQ=255)
                                                ttAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTg                                                                                                                           <  1:835523/126‑1 (MQ=255)
                                                ttAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTg                                                                                                                           >  2:835523/1‑126 (MQ=255)
                                                ttAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTgggg                                                                                                    <  4:728360/149‑1 (MQ=255)
                                                    tatCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCg                                                                                                >  3:428966/1‑149 (MQ=255)
                                                      tCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgttgtt                                                                                                                                        <  1:124609/107‑1 (MQ=255)
                                                      tCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgttgtt                                                                                                                                        >  2:124609/1‑107 (MQ=255)
                                                       cccccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGa                                                                                              >  2:542090/1‑148 (MQ=255)
                                                       cccccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGa                                                                                              >  1:43418/1‑148 (MQ=255)
                                                       cccccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGa                                                                                              >  2:438378/1‑148 (MQ=255)
                                                          ccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTa                                                                                           <  2:681465/148‑1 (MQ=255)
                                                           cGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTAc                                                                                          >  2:194151/1‑148 (MQ=255)
                                                           cGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACt                                                                                         <  4:701630/149‑1 (MQ=255)
                                                           cGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGAGAAACCGGTAGAAGCTCTGGGGGCCGGATTAc                                                                                          <  2:168446/148‑1 (MQ=255)
                                                              cATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTgg                                                                                       >  2:744188/1‑148 (MQ=255)
                                                               atatGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGga                                                                                      <  1:827048/148‑1 (MQ=255)
                                                                tatGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGgat                                                                                     >  1:47852/1‑148 (MQ=255)
                                                                 atGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGgatg                                                                                    <  2:134989/148‑1 (MQ=255)
                                                                   gCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGgatgat                                                                                  >  2:595808/1‑148 (MQ=255)
                                                                    cTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGatgatg                                                                                 >  1:855479/1‑148 (MQ=255)
                                                                      gTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGaa                                                                               <  1:817505/148‑1 (MQ=255)
                                                                       tttGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGaatc                                                                             >  4:819541/1‑147 (MQ=255)
                                                                       tttGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAAc                                                                             >  4:517034/1‑149 (MQ=255)
                                                                       tttGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTACCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAAc                                                                             >  3:665274/1‑149 (MQ=255)
                                                                       tttGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAACCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAAc                                                                             >  4:908870/1‑149 (MQ=255)
                                                                           aTGTGCAAATGCCTGGCACCGTCCTGGGACATATCACACCACAAGCCGCACTGTCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAAcggc                                                                          >  1:501617/1‑148 (MQ=255)
                                                                           aTGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAAcggcg                                                                         <  4:476575/149‑1 (MQ=255)
                                                                           aTGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAAcggc                                                                          <  1:542090/148‑1 (MQ=255)
                                                                                  aaTGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAAttt                                                                  <  4:290298/149‑1 (MQ=255)
                                                                                  aaTGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAAtt                                                                   >  1:220730/1‑148 (MQ=255)
                                                                                      ccTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTtctc                                                               <  1:471272/148‑1 (MQ=255)
                                                                                       cTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTctct                                                              <  1:595808/148‑1 (MQ=255)
                                                                                       cTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTctct                                                              <  1:907918/148‑1 (MQ=255)
                                                                                        tGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCtt                                                             <  2:649798/148‑1 (MQ=255)
                                                                                         ggCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTg                                                            >  1:264615/1‑148 (MQ=255)
                                                                                           cACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGc                                                          >  2:496081/1‑148 (MQ=255)
                                                                                                tCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCTCGGCAGGACTGCCGGttgttt                                                                                                                                       >  4:587943/1‑66 (MQ=255)
                                                                                                tCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCTCGGCAGGACTGCCGGttgttt                                                                                                                                       <  3:587943/66‑1 (MQ=255)
                                                                                                tCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgttgtt                                                                                                                                        <  2:297334/65‑1 (MQ=255)
                                                                                                tCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgttgtt                                                                                                                                        >  1:297334/1‑65 (MQ=255)
                                                                                                   tCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACAt                                                  >  1:409902/1‑148 (MQ=255)
                                                                                                     ggACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGc                                               >  4:843095/1‑149 (MQ=255)
                                                                                                      gACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGc                                               <  2:332112/148‑1 (MQ=255)
                                                                                                        cATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAATCCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCAtt                                            <  4:768832/149‑1 (MQ=255)
                                                                                                          tatCACACCACAAGCCGCACTGGCGACACAGTTCCCGGCAGGACTGCCGGTTGTTTGTACCTCCAGTGATAAACCGGTAGAAGCTCTGGGGGGCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtga                                          <  3:843095/149‑1 (MQ=255)
                                                                                                           atCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgat                                         <  3:355044/149‑1 (MQ=255)
                                                                                                           atCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgat                                         <  4:387019/149‑1 (MQ=255)
                                                                                                           atCACACCACAAGCCGCACTGGCGACACAATTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgat                                         <  3:542513/149‑1 (MQ=255)
                                                                                                            tcaatccacAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgat                                         >  2:65709/6‑148 (MQ=255)
                                                                                                             cacaCCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGga                                                                                      >  1:732945/1‑102 (MQ=255)
                                                                                                             cacaCCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGga                                                                                      <  2:732945/102‑1 (MQ=255)
                                                                                                              acacCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatgat                                      >  3:587461/1‑149 (MQ=255)
                                                                                                              acacCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatgat                                      <  3:908870/149‑1 (MQ=255)
                                                                                                              acacCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatga                                       <  2:780092/148‑1 (MQ=255)
                                                                                                              acacCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatga                                       <  1:700292/148‑1 (MQ=255)
                                                                                                               caccacAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatgat                                      >  1:927631/1‑148 (MQ=255)
                                                                                                                 ccacAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGaa                                   <  4:52434/149‑1 (MQ=255)
                                                                                                                  cacaAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAAc                                  >  4:219659/1‑149 (MQ=255)
                                                                                                                  cacaAGCCGCACTGGCGACACATATCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAAc                                  >  4:61780/1‑149 (MQ=255)
                                                                                                                   acaAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAAc                                  >  1:465999/1‑148 (MQ=255)
                                                                                                                        ccGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCaa                             <  1:5296/148‑1 (MQ=255)
                                                                                                                                  gACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGa                  >  3:885309/1‑149 (MQ=255)
                                                                                                                                    cacaTTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGaa                 <  1:744188/148‑1 (MQ=255)
                                                                                                                                    cacaTTTCCCGGCAGGACTGCCGATTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGaaa                <  4:665274/149‑1 (MQ=255)
                                                                                                                                      aaacaaccGGCAGTCCTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGa              >  4:613584/7‑149 (MQ=255)
                                                                                                                                       aTTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGa              <  2:820631/148‑1 (MQ=255)
                                                                                                                                       aTTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGAt             <  3:517034/149‑1 (MQ=255)
                                                                                                                                        tttCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGAt             >  2:695315/1‑148 (MQ=255)
                                                                                                                                         ttCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGGACGGCAAAGCACTGCCGAAAGATc            >  1:662438/1‑148 (MQ=255)
                                                                                                                                          tCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATcc           <  2:702233/148‑1 (MQ=255)
                                                                                                                                          aaccGGCAGTCCTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATcc           >  2:147619/3‑148 (MQ=255)
                                                                                                                                          aaccGGCAGTCCTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCa          >  3:398738/3‑149 (MQ=255)
                                                                                                                                           cccGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCa          <  2:488697/148‑1 (MQ=255)
                                                                                                                                           cccGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCa          <  1:438378/148‑1 (MQ=255)
                                                                                                                                           accGGCAGTCCTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCa          >  1:659958/2‑148 (MQ=255)
                                                                                                                                            ccGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGATCTGGGGGCCGGATTAGTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAAGGGCAAAGCACTGCCGAAAGATCCAGt        >  4:18740/1‑149 (MQ=255)
                                                                                                                                               gCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAg      <  2:43418/148‑1 (MQ=255)
                                                                                                                                                  ggACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGCGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTa  >  3:114171/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATTATAAGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTA  >  NC_000913/2902347‑2902641

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