Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 897,081 C→T 5.8% H316H (CAC→CAT potH → putrescine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913897,0810CT 5.8% 222.4 / 2.7 69H316H (CAC→CATpotHputrescine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (31/34);  new base T (2/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

ATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGG  >  NC_000913/896935‑897224
                                                                                                                                                  |                                                                                                                                               
atcatTATCGGGCGCGTTCTATGGCAAGAGTTTTTCAATAACCGAGACTGGCCGGTGTCCTCGGCGGTAGCGATCATCATTTTTCTGCTGATAATTGTTCCGATAATGTGGTTCCAAAAACACCAGCAAAAAAGCGTGGGAGAACACg                                                                                                                                                >  4:209542/1‑148 (MQ=255)
 tcatGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACg                                                                                                                                                >  8:45061/1‑147 (MQ=255)
   atGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGa                                                                                                                                              <  6:213062/147‑1 (MQ=255)
    tGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGAt                                                                                                                                             >  3:100893/1‑147 (MQ=255)
      aTCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGaa                                                                                                                                          >  7:346244/1‑148 (MQ=255)
             gcgTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTTCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTa                                                                                                                                   <  8:313502/148‑1 (MQ=255)
                       gCATGAGTTTTTCAATAACCGCGACTGGCCGTTGGACTCGGCGGTAGCGATCATCATGTTGCTGCTGGTAATTGTGCCGATAATGTTGTTCCACAAAGACCAGCAAAAAAGCGTTTGAGAACACGGATGAATAATTTACCGGTAgtt                                                                                                                          <  3:209542/147‑1 (MQ=255)
                        cAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGTATAATTTACCGGTAgttcg                                                                                                                        >  7:319129/1‑148 (MQ=255)
                          agagTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAgttcgtt                                                                                                                      <  7:45061/148‑1 (MQ=255)
                           gagTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAgttcgttc                                                                                                                     <  6:328124/148‑1 (MQ=255)
                               ttttCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGccc                                                                                                                 >  1:8889/1‑148 (MQ=255)
                                    aaTAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACgg                                                                                                                                               <  6:43670/113‑1 (MQ=255)
                                    aaTAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACgg                                                                                                                                               >  5:43670/1‑113 (MQ=255)
                                        aCCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTTCCGATAATGTGGGTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGAt                                                                                                         <  1:63749/147‑1 (MQ=255)
                                        aCCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGAt                                                                                                         <  4:263668/147‑1 (MQ=255)
                                             gACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGAt                                                                                                   <  4:100893/148‑1 (MQ=255)
                                                  gCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGGGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGAtttt                                                                                                <  6:401188/146‑1 (MQ=255)
                                                     ggTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGTATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTtgctgc                                                                                           <  8:319129/148‑1 (MQ=255)
                                                     ggTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCAGGCGGATTGTGATTTtgctgc                                                                                           >  3:203959/1‑148 (MQ=255)
                                                         gCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAAGAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTggg                                                                                       >  4:272133/1‑148 (MQ=255)
                                                         gCCGCGGCGGTAGCGATCATGATGTTGCTGCTGCTAATTGTGCCGTTAATGTGGATCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCAGATTGTGATTTTGCTGCTgg                                                                                        <  7:156380/147‑1 (MQ=255)
                                                                       gATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGt                                                                                                                              <  7:234883/95‑1 (MQ=255)
                                                                       gATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGt                                                                                                                              >  8:234883/1‑95 (MQ=255)
                                                                        atcttcATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGAttt                                                                                                 >  7:297923/1‑123 (MQ=255)
                                                                        atcttcATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGAttt                                                                                                 <  8:297923/123‑1 (MQ=255)
                                                                         tcatcaTGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCAACAAACACCAGCAAAAAAGCGTGGGAGAACACGCATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTAc                                                                       >  2:164985/1‑148 (MQ=255)
                                                                                       ctgctAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAatgctg                                                           <  7:366480/146‑1 (MQ=255)
                                                                                        tgctAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAatgctgatg                                                        >  5:370086/1‑148 (MQ=255)
                                                                                        tgctAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAatgctgatg                                                        >  3:288968/1‑148 (MQ=255)
                                                                                         gctAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTg                                                                                                               <  5:282260/92‑1 (MQ=255)
                                                                                         gctAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTg                                                                                                               >  6:282260/1‑92 (MQ=255)
                                                                                         gctAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAatgctgatgc                                                       >  4:398042/1‑148 (MQ=255)
                                                                                          ctAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTctct                                                                         >  2:65837/1‑129 (MQ=255)
                                                                                          ctAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTctct                                                                         <  1:65837/129‑1 (MQ=255)
                                                                                          ctAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAatgctgatgct                                                      >  2:438895/1‑148 (MQ=255)
                                                                                            aaTTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTgg                                                    >  7:80176/1‑148 (MQ=255)
                                                                                             aTTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCTGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTgg                                                    <  8:80176/147‑1 (MQ=255)
                                                                                             aTTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCAGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTgg                                                    <  4:203959/147‑1 (MQ=255)
                                                                                               tgtgCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGcc                                                                  <  1:396775/131‑1 (MQ=255)
                                                                                               tgtgCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGcc                                                                  >  2:396775/1‑131 (MQ=255)
                                                                                                 tgCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATc                                               <  2:8889/148‑1 (MQ=255)
                                                                                                  gCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATc                                               >  2:157045/1‑147 (MQ=255)
                                                                                                     gATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTAtt                                           >  6:135949/1‑148 (MQ=255)
                                                                                                     gATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTAtt                                           <  1:438895/148‑1 (MQ=255)
                                                                                                     gATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTAtt                                           >  7:220861/1‑148 (MQ=255)
                                                                                                     gATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTGCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTAtt                                           <  5:135949/148‑1 (MQ=255)
                                                                                                       tAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAatgctg                                                           <  6:69408/130‑1 (MQ=255)
                                                                                                       tAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAatgctg                                                           >  5:69408/1‑130 (MQ=255)
                                                                                                           gtgGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTGTACCTTTCTCTCCGCGCCAATGCTGATGCTGGTTATCTATTCGttt                                      >  2:315926/1‑147 (MQ=255)
                                                                                                            ttgTTCCACAAACACCAGCAAAAAATCGTGGGATACACCTTATGAATAATTTACCGGTAGTTCGTGCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATTCTGATGCTGGTTATCTATTCGTTTa                                     <  6:363473/145‑1 (MQ=255)
                                                                                                            tgGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACAGGTAGTTAGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTa                                     <  6:370086/147‑1 (MQ=255)
                                                                                                                    caAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCg                             >  3:394777/1‑147 (MQ=255)
                                                                                                                        cacCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCa                                                                 <  2:46268/107‑1 (MQ=255)
                                                                                                                        cacCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCa                                                                 >  1:46268/1‑107 (MQ=255)
                                                                                                                        cacCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCt                        >  5:184409/1‑148 (MQ=255)
                                                                                                                          ccAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCa                                                                 >  5:187767/1‑105 (MQ=255)
                                                                                                                          ccAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCa                                                                 <  6:187767/105‑1 (MQ=255)
                                                                                                                            aGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGt                     >  5:368209/1‑147 (MQ=255)
                                                                                                                            aGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTg                    <  1:164985/148‑1 (MQ=255)
                                                                                                                             gCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGa                   >  7:220736/1‑148 (MQ=255)
                                                                                                                                aaaaaGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACgg                <  8:220736/148‑1 (MQ=255)
                                                                                                                                 aaaaGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGgt               <  1:157045/148‑1 (MQ=255)
                                                                                                                                    aGCGTGGGAGAACATGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGtgtg            >  7:174193/1‑148 (MQ=255)
                                                                                                                                    aGCGTGGGAGAACATGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGtgtg            >  7:320389/1‑148 (MQ=255)
                                                                                                                                     gCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTgg           <  8:220861/148‑1 (MQ=255)
                                                                                                                                           gagaACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGc     <  6:368209/148‑1 (MQ=255)
                                                                                                                                            agaACATGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCt    <  8:174193/148‑1 (MQ=255)
                                                                                                                                            agaACATGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCt    <  8:320389/148‑1 (MQ=255)
                                                                                                                                              aaCACGGATGAATAATTTACCGGTAGTGCGGTCGCCCTGGCGGATTGTGAGTTTGCTGCGGGGCTTTACCTTTCTCTACGCGCCACTGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTgg  <  2:283784/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
ATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGG  >  NC_000913/896935‑897224

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