Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,255,991 A→T 6.2% K213* (AAA→TAA) ‡ yeiL → putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,255,9910AT 6.2% 333.6 / 8.8 95K213* (AAA→TAA) ‡yeiLputative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (47/42);  new base T (4/2);  total (51/44)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

GAGTATTTAGGCGTTTCTTATCGACATCTTTTATATGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATAAAGAACCCGGCAATAGTCGCACTCATCAGGTTGGAAA  >  NC_000913/2255844‑2256137
                                                                                                                                                   |                                                                                                                                                  
gagTATTTAGGCGTTTCTTATCGACATCTTTTATATGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATaa                                                                                                                                                   <  3:757010/149‑1 (MQ=255)
gagTATTTAGGCGTTTCTTATCGACATCTTTTATATGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATAGACCCGGAGAATaa                                                                                                                                                   >  2:639754/1‑149 (MQ=255)
 agTATTTAGGCGTTTCTTATCGACATCTTTTATATGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATaaa                                                                                                                                                  <  3:933999/149‑1 (MQ=255)
        aGGCGTTTCTTATCGACATCTTTTATATGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCg                                                                                                                                           <  3:635673/149‑1 (MQ=255)
          gCGTTTCTTATCGACATCTTTTATATGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCggg                                                                                                                                         >  2:960515/1‑149 (MQ=255)
          gCGTTTCTTATCGACATCTTTTATATGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCggg                                                                                                                                         >  2:146450/1‑149 (MQ=255)
             tttCTTATCGACATCTTTTATATGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGgatg                                                                                                                                      >  3:421443/1‑149 (MQ=255)
              ttCTTATCGACATCTTTTATATGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGgatga                                                                                                                                     <  1:281388/149‑1 (MQ=255)
                cTTATCGACATCTTTTATATGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATTTATTCTCCgggtccatc                                                                                                                                   <  3:786499/149‑7 (MQ=255)
                 ttATCGACATCTTTTATATGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGc                                                                                                                                  >  3:120840/1‑149 (MQ=255)
                         aTCTTTTATATGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTaaaaa                                                                                                                          >  4:280019/1‑149 (MQ=255)
                         aTCTTTTATATGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTaaaaa                                                                                                                          >  2:539339/1‑149 (MQ=255)
                         aTCTTTTATATGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTaaaaa                                                                                                                          >  2:133411/1‑149 (MQ=255)
                            ttttATATGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAatta                                                                                                                       >  2:985087/1‑149 (MQ=255)
                            ttttATATGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAatta                                                                                                                       <  2:619831/149‑1 (MQ=255)
                            ttttATATGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAatta                                                                                                                       >  2:540172/1‑149 (MQ=255)
                                atatGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTAtttt                                                                                                                   <  1:371683/149‑1 (MQ=255)
                                atatGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTAtttt                                                                                                                   <  1:540172/149‑1 (MQ=255)
                                atatGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTAtttt                                                                                                                   >  2:612051/1‑149 (MQ=255)
                                   tGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCaa                                                                                                                <  4:120840/149‑1 (MQ=255)
                                    gTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAAt                                                                                                               >  3:223221/1‑149 (MQ=255)
                                     ttCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATa                                                                                                              >  3:749526/1‑149 (MQ=255)
                                       ctcGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAgc                                                                                                            <  1:124865/149‑1 (MQ=255)
                                       ctcGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAgc                                                                                                            >  3:718862/1‑149 (MQ=255)
                                          gcgcAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGcgcg                                                                                                         >  4:1007098/1‑149 (MQ=255)
                                            gcAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGAt                                                                                                       <  4:564311/149‑1 (MQ=255)
                                                ttcattcaCGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGcc                                                                                                   <  3:795832/149‑1 (MQ=255)
                                                  cattcaCGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCgg                                                                                                 <  4:65041/149‑1 (MQ=255)
                                                    ttcaCGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGAt                                                                                               <  1:1073364/149‑1 (MQ=255)
                                                     tcaCGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATg                                                                                              >  1:592721/1‑149 (MQ=255)
                                                     tcaCGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATg                                                                                              <  4:5166/149‑1 (MQ=255)
                                                       acgacgGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCa                                                                                            >  3:183279/1‑149 (MQ=255)
                                                         gacgGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAAc                                                                                          >  2:437244/1‑149 (MQ=255)
                                                            ggTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCt                                                                                       <  4:673841/149‑1 (MQ=255)
                                                             gTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCtt                                                                                      >  1:532251/1‑149 (MQ=255)
                                                             gTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCtt                                                                                      >  1:260620/1‑149 (MQ=255)
                                                               ttattaATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTaa                                                                                    >  1:819333/1‑149 (MQ=255)
                                                                    aattaaAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGt                                                                               >  4:804569/1‑149 (MQ=255)
                                                                     attaaAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGtt                                                                              <  2:322992/149‑1 (MQ=255)
                                                                     attaaAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGtt                                                                              <  2:812915/149‑1 (MQ=255)
                                                                     attaaAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGtt                                                                              <  1:133411/149‑1 (MQ=255)
                                                                      ttaaAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGttt                                                                             <  4:103216/149‑1 (MQ=255)
                                                                      ttaaAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGttt                                                                             >  3:48855/1‑149 (MQ=255)
                                                                      ttaaAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGttt                                                                             <  2:592721/149‑1 (MQ=255)
                                                                       taaAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGtttt                                                                            >  4:650429/1‑149 (MQ=255)
                                                                         aaaGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTAt                                                                          >  3:41057/1‑149 (MQ=255)
                                                                         aaaGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTAt                                                                          >  4:899745/1‑149 (MQ=255)
                                                                         aaaGCAAGAAAGGGTATCTCATCAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTAt                                                                          >  1:621133/1‑149 (MQ=255)
                                                                          aaGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATc                                                                         <  3:20135/149‑1 (MQ=255)
                                                                           aGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGCTGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATcc                                                                        >  2:560263/1‑149 (MQ=255)
                                                                           aGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATcc                                                                        >  1:625056/1‑149 (MQ=255)
                                                                            gCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCg                                                                       >  1:446287/1‑149 (MQ=255)
                                                                              aaGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGt                                                                     >  3:1034292/1‑149 (MQ=255)
                                                                               aGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTc                                                                    <  2:833111/149‑1 (MQ=255)
                                                                                gAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCt                                                                   >  3:1051599/1‑149 (MQ=255)
                                                                                  aaGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTAc                                                                 <  4:156796/149‑1 (MQ=255)
                                                                                      gTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAAt                                                             <  1:326036/149‑1 (MQ=255)
                                                                                      gTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAAt                                                             >  2:1012329/1‑149 (MQ=255)
                                                                                       tATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATc                                                            <  1:612051/149‑1 (MQ=255)
                                                                                         tctcATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAATTTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGa                                                          <  3:63617/149‑1 (MQ=255)
                                                                                         tctcATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGa                                                          <  1:860273/149‑1 (MQ=255)
                                                                                         tctcATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGa                                                          <  3:577378/149‑1 (MQ=255)
                                                                                         tctcATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGa                                                          <  4:730999/149‑1 (MQ=255)
                                                                                         tctcATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGa                                                          <  2:1034624/149‑1 (MQ=255)
                                                                                          ctcATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAg                                                         <  3:804569/149‑1 (MQ=255)
                                                                                           tcATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGc                                                        >  2:499956/1‑149 (MQ=255)
                                                                                            cATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCa                                                       >  1:927546/1‑149 (MQ=255)
                                                                                             aTTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCAt                                                      >  4:896688/1‑149 (MQ=255)
                                                                                                 aaaaCAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACg                                                  >  1:288954/1‑149 (MQ=255)
                                                                                                    aCAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAattatt                                                                                                                     >  1:979997/1‑79 (MQ=255)
                                                                                                    aCAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAattatt                                                                                                                     <  2:979997/79‑1 (MQ=255)
                                                                                                       gAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGAc                                            <  2:625056/149‑1 (MQ=255)
                                                                                                       gAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGAc                                            >  3:159938/1‑149 (MQ=255)
                                                                                                              gTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATaaa                                     <  1:146450/149‑1 (MQ=255)
                                                                                                              gTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATaaa                                     <  4:16746/149‑1 (MQ=255)
                                                                                                               ttGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATAAAg                                    <  3:1007098/149‑1 (MQ=255)
                                                                                                                    aGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAAc                                                                                          >  2:521028/1‑90 (MQ=255)
                                                                                                                    aGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAAc                                                                                          <  1:521028/90‑1 (MQ=255)
                                                                                                                     ggACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATAAAGAACCCg                              >  1:144125/1‑149 (MQ=255)
                                                                                                                     ggACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATAAAGAACCCg                              <  1:539339/149‑1 (MQ=255)
                                                                                                                         tGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATAAAGAACCCGGCaa                          >  4:450141/1‑149 (MQ=255)
                                                                                                                           gcgcTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATAAAGAACCCGGCAATa                        <  2:314137/149‑1 (MQ=255)
                                                                                                                             gcTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATAAAGAACCCGGCAATAGt                      >  4:606506/1‑149 (MQ=255)
                                                                                                                              cTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATAAAGAACCCGGCAATAGTc                     <  1:873056/149‑1 (MQ=255)
                                                                                                                              cTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATAAAGAACCCGGCAATAGTc                     <  4:48855/149‑1 (MQ=255)
                                                                                                                                  agaTGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATAAAGAACCCGGCAATAGTCGCAc                 <  3:513649/149‑1 (MQ=255)
                                                                                                                                        tcccGGAGAATTTATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATAAAGAACCCGGCAATAGTCGCACtcatca           >  1:636771/2‑149 (MQ=255)
                                                                                                                                        aCCCGGATAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATAAAGAACCCGGCAATAGTCGCACtcatca           >  2:593095/1‑149 (MQ=255)
                                                                                                                                          ccGGAGAATTTATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGAt                                                                                               <  1:480954/63‑1 (MQ=38)
                                                                                                                                          ccGGAGAATTTATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGAt                                                                                               >  2:480954/1‑63 (MQ=38)
                                                                                                                                          ccGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATAAAGAACCCGGCAATAGTCGCACTCATCAgg         >  4:56560/1‑149 (MQ=255)
                                                                                                                                           cGGAGAATTTATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATAAAGAACCCGGCAATAGTCGCACTCATCAGGt        >  1:73059/1‑149 (MQ=255)
                                                                                                                                           cGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATAAAGAACCCGGCAATAGTCGCACTCATCAGGt        >  3:789293/1‑149 (MQ=255)
                                                                                                                                            ggAGAATTTATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATAAAGAACCCGGCAATAGTCGCACTCATCAGGtt       >  1:17811/1‑149 (MQ=255)
                                                                                                                                            ggAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATAAAGAACCCGGCAATAGTCGCACTCATCAGGtt       <  1:225780/149‑1 (MQ=255)
                                                                                                                                            ggAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATAAAGAACCCGGCAATAGTCGCACTCATCAGGtt       <  2:288954/149‑1 (MQ=255)
                                                                                                                                                aaTAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATAAAGAACCCGGCAATAGTCGCACTCATCAGGTTGGaa   >  1:1104688/1‑149 (MQ=255)
                                                                                                                                                 aTAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATAAAGAACCCGGCAATAGTCGCACTCATCAGGTTGGaaa  >  2:223919/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GAGTATTTAGGCGTTTCTTATCGACATCTTTTATATGTTCTCGCGCAGTTCATTCACGACGGTTTATTAATTAAAAGCAAGAAAGGGTATCTCATTAAAAACAGAAAGCAGTTGTCAGGACTGGCGCTGGAGATGGACCCGGAGAATAAATTCTCCGGGATGATGCAGTAAAAATTATTTTCAATAGCGCGATTGCCGGATGCAACGCTTAATACGTTTTATCCGGTCTACAAATCGAGCATTACGCCAGACCAATAAAGAACCCGGCAATAGTCGCACTCATCAGGTTGGAAA  >  NC_000913/2255844‑2256137

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