Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 340,923 C→A 75.2% L253L (CTC→CTA yahI → carbamate kinase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913340,9230CA75.2% 178.1 / 79.0 97L253L (CTC→CTAyahIcarbamate kinase‑like protein
Reads supporting (aligned to +/- strand):  ref base C (12/12);  new base A (41/32);  total (53/44)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

AGATTACCAAAGCGTGGACGCGGTTATCGACAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTAGAACAAGGCGGCAAAGAAGTGATTATCACCACGCCGGAATGCCTGC  >  NC_000913/340779‑341068
                                                                                                                                                |                                                                                                                                                 
agaTTACCAAAGCGTGGACGCGGTTATCGACAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATc                                                                                                                                               <  4:363380/149‑1 (MQ=255)
 gaTTACCAAAGCGTGGACGCGGTTATCGACAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCg                                                                                                                                              >  4:228776/1‑149 (MQ=255)
    tACCAAAGCGTGGACGCGGTTATCGACAAAGATTTCTCTACCGCGCTGCTGACCCGTGAAATTCACGCCGACCTTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGATAGATCGGGt                                                                                                                                           <  2:234156/149‑1 (MQ=255)
    tACCAAAGCGTGGACGCGGTTATCGACAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCAGCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGt                                                                                                                                           >  2:802994/1‑149 (MQ=255)
    tACCAAAGCGTGGACGCGGTTATCGACAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGATCGGGt                                                                                                                                           <  1:924849/149‑1 (MQ=255)
    tACCAAAGCGTGGACGCGGTTATCGACAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGt                                                                                                                                           <  2:478177/149‑1 (MQ=255)
    tACCAAAGCGTGGACGCGGTTATCGACAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGt                                                                                                                                           <  4:906979/149‑1 (MQ=255)
      ccAAAGCGTGGACGCGGTTATCGACAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGTTCGGGTgg                                                                                                                                         >  2:766036/1‑149 (MQ=255)
      ccAAAGCGTGGACGCGGTTATCGACAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTgg                                                                                                                                         >  4:738054/1‑149 (MQ=255)
       cAAAGCGTGGACGCGGTTATCGACAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGa                                                                                                                                        >  4:715904/1‑149 (MQ=255)
         aaGCGTGGACGCGGTTATCGACAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGata                                                                                                                                      >  2:808620/1‑149 (MQ=255)
         aaGCGTGGACGCGGTTATCGACAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGata                                                                                                                                      <  3:355773/149‑1 (MQ=255)
          aGCGTGGACGCGGTTATCGACAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGatat                                                                                                                                     <  1:50247/149‑1 (MQ=255)
            cGTGGACGCGGTTATCGACAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTg                                                                                                                                   <  2:872294/149‑1 (MQ=255)
               ggACGCGGTTATCGACAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGATCGGGTGGATATTGcca                                                                                                                                <  3:100287/149‑1 (MQ=255)
               ggACGCGGTTATCGACAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGcca                                                                                                                                <  2:502996/149‑1 (MQ=255)
                         aTCGACAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGATCGGGTGGATATTGCCACCATGACCCg                                                                                                                      >  3:534162/1‑149 (MQ=255)
                          tCGACAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGAccc                                                                                                                       <  4:233403/147‑1 (MQ=255)
                          tCGACAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGAccc                                                                                                                       >  3:233403/1‑147 (MQ=255)
                           cGACAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTATGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCt                                                                                                                    >  3:650586/1‑149 (MQ=255)
                            gACAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCta                                                                                                                   >  3:871096/1‑149 (MQ=255)
                             aCAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTg                                                                                                                                          <  4:580880/125‑1 (MQ=255)
                             aCAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTg                                                                                                                                          >  3:580880/1‑125 (MQ=255)
                               aaaGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTatat                                                                                                                >  4:197477/1‑149 (MQ=255)
                               aaaGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTatat                                                                                                                >  3:829347/1‑149 (MQ=255)
                               aaaGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTatat                                                                                                                >  1:148547/1‑149 (MQ=255)
                                  gATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGATCGGGTGGATATTGCCACCATGACCCGCTATATGCa                                                                                                             <  3:52842/149‑1 (MQ=255)
                                   aTCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAg                                                                                                            >  3:61532/1‑149 (MQ=255)
                                   aTCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAg                                                                                                            <  2:81023/149‑1 (MQ=255)
                                     ctctctACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATc                                                                                                                                               <  3:661758/112‑1 (MQ=255)
                                     ctctctACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATc                                                                                                                                               >  4:661758/1‑112 (MQ=255)
                                     ctctctACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGa                                                                                                          >  4:295280/1‑149 (MQ=255)
                                      tctctACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGaa                                                                                                         >  1:182809/1‑149 (MQ=255)
                                      tctctACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGaa                                                                                                         >  2:390330/1‑149 (MQ=255)
                                       ctctACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAg                                                                                                        <  1:390330/149‑1 (MQ=255)
                                        tctACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAgg                                                                                                       <  1:575387/149‑1 (MQ=255)
                                         ctACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAggg                                                                                                      <  4:464749/149‑1 (MQ=255)
                                          tACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGc                                                                                                     <  1:95925/149‑1 (MQ=255)
                                            ccGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCAt                                                                                                   <  2:898630/149‑1 (MQ=255)
                                                  tgctgGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTcccg                                                                                             >  4:356371/1‑149 (MQ=255)
                                                  tgctgGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTcccg                                                                                             >  1:777408/1‑149 (MQ=255)
                                                     tgGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTcccgccc                                                                                          >  2:632050/1‑149 (MQ=255)
                                                     tgGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTcccgccc                                                                                          >  4:709092/1‑149 (MQ=255)
                                                      ggCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTcccgcccg                                                                                         >  4:603231/1‑149 (MQ=255)
                                                         ccGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGgca                                                                                      >  4:169323/1‑149 (MQ=255)
                                                          cGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGa                                                                                                                                                 <  3:412892/89‑1 (MQ=255)
                                                          cGTGAAATTCACGCCGACATGCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGa                                                                                                                                                 >  4:412892/1‑89 (MQ=255)
                                                              aaaTTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATg                                                                                 <  3:751283/149‑1 (MQ=255)
                                                                aTTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGtt                                                                               <  3:693741/149‑1 (MQ=255)
                                                                aTTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGtt                                                                               >  2:655752/1‑149 (MQ=255)
                                                                 ttCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTg                                                                              >  2:853591/1‑149 (MQ=255)
                                                                  tCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGc                                                                             >  2:146864/1‑149 (MQ=255)
                                                                  tCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTATGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGc                                                                             <  4:650586/149‑1 (MQ=255)
                                                                         gACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAat                                                                      >  4:320022/1‑149 (MQ=255)
                                                                         gACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAat                                                                      >  2:507596/1‑149 (MQ=255)
                                                                           cATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAatca                                                                    <  4:93158/149‑1 (MQ=255)
                                                                            aTTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAatcat                                                                   >  4:550192/1‑149 (MQ=255)
                                                                               cTTGGGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGc                                                                >  4:765882/1‑149 (MQ=255)
                                                                                ttGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGcc                                                               >  3:99553/1‑149 (MQ=255)
                                                                                ttGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGcc                                                               >  4:108818/1‑149 (MQ=255)
                                                                                 tgtgATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCa                                                              <  4:829347/149‑1 (MQ=255)
                                                                                     aTCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCt                                                          >  1:498427/1‑149 (MQ=255)
                                                                                      tCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTg                                                         >  2:904303/1‑149 (MQ=255)
                                                                                           acTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACAtt                                                    >  3:450500/1‑149 (MQ=255)
                                                                                             tGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACAtttt                                                  <  3:197477/149‑1 (MQ=255)
                                                                                              ggCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACAttttt                                                 <  3:108818/149‑1 (MQ=255)
                                                                                               gCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTa                                                <  3:709092/149‑1 (MQ=255)
                                                                                               gCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTa                                                <  1:655752/149‑1 (MQ=255)
                                                                                                    gAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGtt                                                                               <  3:640407/113‑1 (MQ=255)
                                                                                                    gAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGtt                                                                               >  4:640407/1‑113 (MQ=255)
                                                                                                    gAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCaaaa                                                                        <  2:295082/120‑1 (MQ=255)
                                                                                                    gAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCaaaa                                                                        >  1:295082/1‑120 (MQ=255)
                                                                                                    gAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTAGAACa                                           >  3:532432/1‑149 (MQ=255)
                                                                                                     aaaaaGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGt                                                                                >  2:386578/1‑111 (MQ=255)
                                                                                                     aaaaaGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGt                                                                                <  1:386578/111‑1 (MQ=255)
                                                                                                     aaaaaGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTAGAACaa                                          <  1:632050/149‑1 (MQ=255)
                                                                                                      aaaaGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTAGAACAAg                                         <  4:61532/149‑1 (MQ=255)
                                                                                                             tgtATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTAGAACAAGGCGGCaa                                  >  2:479609/1‑149 (MQ=255)
                                                                                                               tATTCACTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTAGAACAAGGCGGCAaag                                >  3:66210/1‑149 (MQ=255)
                                                                                                                 ttCCCTTTGGCAAACCGCAGCAGCAGGCGCTAGAGCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTAGAACAAGGCGGCAaagaa                              <  4:378867/149‑1 (MQ=255)
                                                                                                                    aCTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGtt                                                                               <  3:740825/97‑1 (MQ=255)
                                                                                                                    aCTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGtt                                                                               >  4:740825/1‑97 (MQ=255)
                                                                                                                    aCTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTAGAACAAGGCGGCAaagaa                              >  1:29001/1‑146 (MQ=255)
                                                                                                                    aCTTTGGCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTAGAACAAGGCGGCAaagaa                              <  2:29001/146‑1 (MQ=255)
                                                                                                                         ggCAAACCGCAGCAGCAGGCGCTCGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTAGAACAAGGCGGCAAAGAAGTGATTAt                      <  1:913218/149‑1 (MQ=255)
                                                                                                                          gCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTAGaa                                             <  1:569882/125‑1 (MQ=255)
                                                                                                                          gCAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTAGaa                                             >  2:569882/1‑125 (MQ=255)
                                                                                                                           cAAACCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCATCCTGACATTTTTAGAACAAGGCGGCAAAGAAGTGATTATca                    <  3:603231/149‑1 (MQ=255)
                                                                                                                            aaaCCGCAGCAGCAGGCGCTCGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTAGAACAAGGCGGCAAAGAAGTGATTATcac                   >  3:509789/1‑149 (MQ=255)
                                                                                                                            aaaCCGCAGCAGCAGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTAGAACAAGGCGGCAAAGAAGTGATTATcac                   >  1:405605/1‑149 (MQ=255)
                                                                                                                                 gcagcagcaGGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTAGAACAAGGCGGCAAAGAAGTGATTATCACCACGc              >  1:124514/1‑149 (MQ=255)
                                                                                                                                      agcagGCGCTCGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTAGAACAAGGCGGCAAAGAAGTGATTATCACCACGCCGGaa         <  4:338801/149‑1 (MQ=255)
                                                                                                                                         agGCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTAGAACAAGGCGGCAAAGAAGTGATTATCACCACGCCGGAAtgc      <  4:450500/149‑1 (MQ=255)
                                                                                                                                          ggCGCTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTAGAACAAGGCGGCAAAGAAGTGATTATCACCACGCCGGAAtgcc     <  3:356371/149‑1 (MQ=255)
                                                                                                                                             gcTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTAGAACAAGGCGGCAAAGAAGTGATTATCACCACGCCGGAAtgcctgc  >  1:623473/1‑149 (MQ=255)
                                                                                                                                             gcTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTAGAACAAGGCGGCAAAGAAGTGATTATCACCACGCCGGAAtgcctgc  >  1:320544/1‑149 (MQ=255)
                                                                                                                                             gcTAGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTAGAACAAGGCGGCAAAGAAGTGATTATCACCACGCCGGAAtgcctgc  >  2:541600/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
AGATTACCAAAGCGTGGACGCGGTTATCGACAAAGATCTCTCTACCGCGCTGCTGGCCCGTGAAATTCACGCCGACATTCTTGTGATCACCACTGGCGTCGAAAAAGTGTGTATTCACTTTGGCAAACCGCAGCAGCAGGCGCTCGATCGGGTGGATATTGCCACCATGACCCGCTATATGCAGGAAGGGCATTTCCCGCCCGGCAGCATGTTGCCAAAAATCATCGCCAGCCTGACATTTTTAGAACAAGGCGGCAAAGAAGTGATTATCACCACGCCGGAATGCCTGC  >  NC_000913/340779‑341068

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: