Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,442,896 A→T F401Y (TTC→TAC)  secY ← preprotein translocase membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,442,8960AT90.5% 187.1 / 10.0 63F401Y (TTC→TAC) secYpreprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (2/4);  new base T (26/31);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.19e-01

GAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATCACACGGATGACACCATCACGCTTAACGATTTTGCAGTTACGGCATAATTTCTTGACGGAAGCACGAACTTTCATTTTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTC  >  NC_000913/3442600‑3443189
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                     
gAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATCACACGGATGACACCATCACGCTTAACGATTTTGCAGTTACGGCATAATTTCTTGACGGAAGCACGAACTTTCATTTTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTTGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCgtagta                                                                                                                                                                                                                                                                                                   >  2:387427/1‑301 (MQ=255)
      tatGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATCACACGGATGACACCATCACGCTTAACGATTTTGCAGTTACGGCATAATTTCTTGACGGAAGCACGAACTTTCATTTTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACgg                                                                                                                                                                                                                                                                                             <  2:122682/301‑1 (MQ=255)
      tatGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGAACTGCAAATCACACGGATGACACCATCAGGCTTAACGATTTTGCAGTTACGGCATAATTTCTTGACGGAAGCACGAACTTTCATGTTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGAACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACgg                                                                                                                                                                                                                                                                                             <  1:221647/301‑1 (MQ=255)
      aatGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATCACACGGATGACACCATCACGCTTAACGATTTTGCAGTTACGGCATAATTTCTTGACGGAAGCACGAACTTTCATTTTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACgg                                                                                                                                                                                                                                                                                             <  2:582143/300‑1 (MQ=255)
          cGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATCACACGGATGACACCATCACGCTTAACGATTTTGCAGTTAAGGCATAATTTCTTGACGGAAGCACGACCTTTCATTTTTACTCTCCGTAACTTCTCGGTCGACCAATCATCGTCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGTCATCATCAGAGTTTTCATTTGAACCATAAAGTCCTTAATCACGACAACAACGATAAGCAGTGAGGACACAACCTAGTAGAACGTTaat                                                                                                                                                                                                                                                                                         >  1:226635/1‑299 (MQ=255)
                                     tGCTTCGGCTCGGCACTGCAAATCACACGGATGACACCATCACGCTTAACGATTTTGCAGTTACGGCATAATTTCTTGACGGAAGCACGAACTTTCATTTTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACTGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCggg                                                                                                                                                                                                                                                              >  2:572382/1‑301 (MQ=255)
                                      gCTTCGGCTCGGCACTGCGAATCACACGGATGACACCATCACGCTTAACGATTTTGCAGTTACGGCATAATTTCTTGACGGAAGCACGAACTTTCATTTTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGa                                                                                                                                                                                                                                                             >  2:286311/1‑301 (MQ=255)
                                         tcggctcggcACTGCAAATCACACGGATGACACCATCACGCTTAACGATTTTGCAGTTACGGCATAATTTCTTGACGGAAGCACGAACTTTCATTTTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCa                                                                                                                                                                                                                                                          >  2:100785/1‑301 (MQ=255)
                                             ctcggcCCTGCAAATCACACGGATGACACCATCACGCTCTACGATTTTGCAGTTACGGCATAATTTCTTGACGGAAGAACGAACTTTCCTTTTTACTCTCCCTAACTTTGCGGGCGACCAATTATCCGCCATAGCCTTTCAGGTTCGCCTTCTTCTATTCGGACTCGTACTGACTGGAAGTCATCAGAATTTGCACTTTAACCATAAAGTCCGTAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCTCGCATGAACTCCGGGATCAGGc                                                                                                                                                                                                                                                       <  1:296051/300‑1 (MQ=255)
                                             ctcggcACTGCAAAGCACACGGATGACACCATCACGCTTAACGATTGTGCAGTTACGGCATAATTTCTTGACGGAAGCACGAACTTTCATTTTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTGTGCACTGGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGc                                                                                                                                                                                                                                                       <  2:596952/300‑1 (MQ=255)
                                                   ccTGCGAATAACACGGATGACACCATCACGCTTAACGATTTTGCAGTTACGGCTTAATTTCTTGACGGAAGCACGAACTTTCCTTTTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCATAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTAGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCAACGTAGTAGAAAGGTACTTTAATTGCATCAAGCATGAACTCCGGGATCAGGCAGATaaa                                                                                                                                                                                                                                                <  1:45689/300‑1 (MQ=255)
                                                          aTCACACGGATGACACCATCACGCTTAACGATTTTGCAGTTACGGCATAATTTCTTGACGGAAGCACGAACTTTCATTTTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAata                                                                                                                                                                                                                                         <  1:236706/301‑1 (MQ=255)
                                                                   aTGACACCATCACGCTTAACGATTTTGCAGTTACGGCATAATTTCTTGACGGAAGCACGAACTTTCATTTTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCa                                                                                                                                                                                                                                >  1:34072/1‑301 (MQ=255)
                                                                    tgaaaacagaaagcgaaaCGATAGTGCAGTGGTGGCATAATTTGTGTACGGAAGCACGAACTTTCAGTTTTACAAGCCTTAACTTCTCAGAAGACCGATAATCGACAGTAGAACTTAAGGTACGCATAATTCAAGGCAGAATCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCGAAAAAGGCAAGAATCACGACAACAACGATAAGCAGTGAGGTCACACCGTAGTAGAACGGTACGTGCATTGCATCACGCATGAACTGAGGGATAAGGCAGATAAAGGTAATATACAGCGCac                                                                                                                                                                                                                               <  1:363171/285‑1 (MQ=255)
                                                                      acacCATCACGCTTAACGATTTTGCAGTTACGGCATAATTTCTTGACGGAAGCACGAACTTTCATTTTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCacca                                                                                                                                                                                                                             >  2:414142/1‑301 (MQ=255)
                                                                      acacCATCACGCTTAACGATTTTGCAGTTACGGCATAATTTCTTGACGGAAGCACGAACTTTCATTTTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCacca                                                                                                                                                                                                                             >  2:630975/1‑301 (MQ=255)
                                                                      acacCATCACGCTTAACGATTTTGCAGTTACGGCATAATTTCTTGACGGAAGCACGAACTTTCATTTTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCacca                                                                                                                                                                                                                             >  2:136343/1‑301 (MQ=255)
                                                                         ccATCACGCTTAACGATTTTGCAGTTACGGCATAATTTCTTGACGGAAGCACGAACTTTCATTTTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCaccaac                                                                                                                                                                                                                           >  2:305552/1‑300 (MQ=255)
                                                                             cACGCTTAACGATTTTGCAGTTACGGCATAATTTCTTGACGGAAGCACGAACTTTCATTTTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAggg                                                                                                                                                                                                                      <  2:648245/301‑1 (MQ=255)
                                                                                                gTTACGGCATAATTTCTTGACGGAAGCACGAACTTTCATTTTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTAAGGCGGGTCATTACttt                                                                                                                                                                                                   >  1:297439/1‑301 (MQ=255)
                                                                                                      gCATAATTTCTTGACGGAAGCACGAACTTTCATTTTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGa                                                                                                                                                                                              <  2:285171/300‑1 (MQ=255)
                                                                                                      gCATAAGTTCTTGACGGAAGCACGAACTTTCATTGTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCATAAATGCAGACTCATACTGACTGGACATCATCAGGGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGGTAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACTCATGAACTCCGGGATCAGGAAGATAAAGGTAATATACAGCGCCCCAACCAGCGTAAGGCGGGTCATTACTTTATCGa                                                                                                                                                                                              <  1:387422/300‑1 (MQ=255)
                                                                                                               cTTGACGGAAGCACGAACTTTCATTTTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGc                                                                                                                                                                                    <  2:297443/301‑1 (MQ=255)
                                                                                                                 tGACGGAAGAACGAACTTTCATTTTTACTCTCCGTAACTTCTCGTGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCg                                                                                                                                                                                  <  1:100784/301‑1 (MQ=255)
                                                                                                                              aaCTTTCATTTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGa                                                                                                                                                                    <  2:310470/301‑1 (MQ=255)
                                                                                                                                          tACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTgg                                                                                                                                                         <  1:414135/301‑1 (MQ=255)
                                                                                                                                          tACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTgg                                                                                                                                                         <  2:384405/301‑1 (MQ=255)
                                                                                                                                              ctcCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTTAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTACTAGAACGCAACCTTCATTGCATCACGCCTGAACTCCGGGATCAGGCAGCTACCGGTAATATACAGCGCACCCACCCGCGCCCGCCGGGTCCTTACCCTCTCGATCTACTTCGCCGTTTGCTCCCCCGCACCAACTccgcctac                                                                                                                                                      >  1:452208/1‑294 (MQ=255)
                                                                                                                                                         ctcGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGAGGGTCATTACTTTATCGATATACTTCGCCGTTTGCGCCACCGGACGAATTCCTGGTACAAATGCACCGGa                                                                                                                                          >  1:409415/1‑301 (MQ=255)
                                                                                                                                                                   cAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAAAATTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGGACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCACACCAACCAGGGTCAGGCGGATCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCAGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAgg                                                                                                                                <  1:572371/301‑1 (MQ=255)
                                                                                                                                                                    aaTTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGt                                                                                                                               <  1:286307/301‑1 (MQ=255)
                                                                                                                                                                        ccgatctgTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATc                                                                                                                           <  2:271086/294‑1 (MQ=255)
                                                                                                                                                                        aTCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTAt                                                                                                                            >  2:358441/1‑300 (MQ=255)
                                                                                                                                                                           ggCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATctg                                                                                                                         >  2:431367/1‑300 (MQ=255)
                                                                                                                                                                               gTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGCCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCCCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATcagatcgg                                                                                                                    >  1:271082/1‑294 (MQ=255)
                                                                                                                                                                                   ccTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTc                                                                                                                  <  2:452215/299‑1 (MQ=255)
                                                                                                                                                                                       tCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCgg                                                                                                            <  2:34072/301‑1 (MQ=255)
                                                                                                                                                                                        cAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATACAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCggg                                                                                                           >  1:24076/1‑301 (MQ=255)
                                                                                                                                                                                           attCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGtt                                                                                                         <  1:41039/299‑1 (MQ=255)
                                                                                                                                                                                            ttCGCCTTCTTCAATGCAGCCTCATACTGACTGGACATCATCCGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTg                                                                                                        <  1:358436/300‑1 (MQ=255)
                                                                                                                                                                                             tCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGaa                                                                                                      >  1:541208/1‑301 (MQ=255)
                                                                                                                                                                                                 cttcttCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAAcc                                                                                                  <  2:24076/301‑1 (MQ=255)
                                                                                                                                                                                                   tcttcAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCaa                                                                                                <  2:574270/301‑1 (MQ=255)
                                                                                                                                                                                                     ttcAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACg                                                                                              <  2:541218/301‑1 (MQ=255)
                                                                                                                                                                                                       aaaTGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGcc                                                                                            <  1:136342/300‑1 (MQ=255)
                                                                                                                                                                                                              gACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACATATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTaga                                                                                     >  2:524898/1‑301 (MQ=255)
                                                                                                                                                                                                                 tCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGGCGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTagaaga                                                                                  <  2:226639/301‑1 (MQ=255)
                                                                                                                                                                                                                     ctcttccgatcttcatcaGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACa                                                                              <  2:414412/289‑1 (MQ=255)
                                                                                                                                                                                                                      ctgactgGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGAGGGTCATTACTTTACCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCATCGGACTTCTACAGGCTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGCAGAAACag                                                                             >  1:139090/1‑301 (MQ=255)
                                                                                                                                                                                                                         actgGAGATCATAAGAGTTTGCACTTGAGCCATAAAGACCAGCATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCTCGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCATGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACagaag                                                                          <  1:279190/301‑1 (MQ=255)
                                                                                                                                                                                                                           tgGACATCAGCAGAGTTTGCACTTGAGCCATAAAATCCATAATCACGACCACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATAAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGAGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACagaagaa                                                                        <  2:139091/301‑1 (MQ=255)
                                                                                                                                                                                                                                 tcatcaGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGGAGaacaaacaagatcggaagag                                                                  >  1:414405/1‑289 (MQ=255)
                                                                                                                                                                                                                                              cACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCATGGTACAAATGCACCGGACTTCTCTAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTcgaagatcggaagagcgtcgggtagggaaagagtg                                                     >  1:128298/1‑272 (MQ=255)
                                                                                                                                                                                                                                               aCTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATaga                                                    >  1:548199/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                    aTCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAATGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGAACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGCCTTCTTCAGGTTATCTGCTGTTTCACTCGGGTCGAAAACCAACGCCGTGTAGAAGCAACTGTAGAAGATGATTGAAGACTCACAGAGTAACACATAAAGCGGTTGc                                >  1:572288/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                    aTCACGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAAATCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAAAAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGc                                >  1:404933/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                       aCGACAACAACGATAAGCAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAg                            >  2:611331/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                     aaGCAGTGAGGTCCCAACGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGGTCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCAATACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACATATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGc              <  1:524888/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        cAGTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAAATCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCg             <  2:404940/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          gTGAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAAt         >  2:128240/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                            gAGGTCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCTGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGAAAATATAGCGAAAttt       >  1:402466/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                tCCCACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAAttgttg    <  1:128239/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                 cccACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGTGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTc  <  2:409422/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                 cccACCGTAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATtgttgt   >  2:115905/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                     
GAAAAATATGCGAAAAAATCAGCCTTGGCGCTGTTTATGCTTCGGCTCGGCACTGCAAATCACACGGATGACACCATCACGCTTAACGATTTTGCAGTTACGGCATAATTTCTTGACGGAAGCACGAACTTTCATTTTTACTCTCCGTAACTTCTCGGGCGACCAATTATCGGCCGTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATATACAGCGCACCAACCAGGGTCAGGCGGGTCATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTC  >  NC_000913/3442600‑3443189

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