Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 399,801 A→G 5.3% Q70R (CAA→CGA)  yaiZ → DUF2754 family putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913399,8010AG 5.3% 406.1 / 9.0 131Q70R (CAA→CGA) yaiZDUF2754 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (56/68);  new base G (4/3);  total (60/71)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01

Only 100 of 138 total aligned reads displayed.

TATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTG  >  NC_000913/399654‑399949
                                                                                                                                                   |                                                                                                                                                    
tatTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCaa                                                                                                                                                     <  2:575712/149‑1 (MQ=255)
tatTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCaa                                                                                                                                                     <  1:983748/149‑1 (MQ=255)
tatTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACATATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCaa                                                                                                                                                     >  1:1417238/1‑149 (MQ=255)
 attcattcTTCATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAAt                                                                                                                                                    <  2:1463189/149‑1 (MQ=255)
    cattcTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAg                                                                                                                                                 >  1:355943/1‑149 (MQ=255)
     attcTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGa                                                                                                                                                >  1:203082/1‑149 (MQ=255)
     attcTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGa                                                                                                                                                >  2:198328/1‑149 (MQ=255)
     attcTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGa                                                                                                                                                >  1:346797/1‑149 (MQ=255)
     attcTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGa                                                                                                                                                <  1:1306440/149‑1 (MQ=255)
      ttcTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGaa                                                                                                                                               >  2:378297/1‑149 (MQ=255)
        cTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAata                                                                                                                                             >  1:680150/1‑149 (MQ=255)
         ttAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTTGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAatat                                                                                                                                            >  2:85314/1‑149 (MQ=255)
         ttAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAatat                                                                                                                                            >  2:1236006/1‑149 (MQ=255)
         ttAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAatat                                                                                                                                            >  2:1460777/1‑149 (MQ=255)
                   ggtggGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGt                                                                                                                                  <  1:1216746/149‑1 (MQ=255)
                    gtggGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTAGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGtt                                                                                                                                 >  1:108400/1‑149 (MQ=255)
                        ggTTACTTGGATTTGAAGCAAAAGGTTGGCAGACCTATGGCGTCGGTCTGGTGACGTGTGTGATTAGGTTCTGGCTGGCGGTGTTGATTATTCGTCGTCGCGATGAAGAAACTGTAAAAGCCCGATAATATTCTTATTGGGCGttttca                                                                                                                             <  2:898261/149‑5 (MQ=255)
                        ggTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTaaaa                                                                                                                             <  1:486581/149‑1 (MQ=255)
                         ggtACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCGATAATATTCTTATTGGGCGttttcag                                                                                                                            <  2:1327031/147‑6 (MQ=255)
                          ttACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATTCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAata                                                                                                                           <  2:330036/149‑1 (MQ=255)
                           tACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAatat                                                                                                                          <  2:1380685/149‑1 (MQ=255)
                             cTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGGTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTa                                                                                                                        <  1:1409681/149‑1 (MQ=255)
                              tGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTAc                                                                                                                       >  2:771037/1‑149 (MQ=255)
                              tGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTAc                                                                                                                       >  2:776329/1‑149 (MQ=255)
                              ggggATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGGTGATTAATCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTAc                                                                                                                       <  1:44582/148‑1 (MQ=255)
                                ggATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATTAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACAt                                                                                                                     >  2:128430/1‑149 (MQ=255)
                                   tttGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATtgt                                                                                                                  >  2:452811/1‑149 (MQ=255)
                                    ttGAAGCAAAAGGTTGGCAGACCTATGCCGTCTGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATtgtg                                                                                                                 <  1:1190272/149‑1 (MQ=255)
                                    ttGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATtgtg                                                                                                                 <  1:54935/149‑1 (MQ=255)
                                     tGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTgg                                                                                                                <  1:1460777/149‑1 (MQ=255)
                                     tGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTgg                                                                                                                <  1:1116817/149‑1 (MQ=255)
                                     tGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTgg                                                                                                                <  2:62168/149‑1 (MQ=255)
                                       aaGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGtt                                                                                                              >  2:1136711/1‑149 (MQ=255)
                                         gCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGtttt                                                                                                            >  2:511063/1‑149 (MQ=255)
                                         gCAAAAGGTTGGCAGACATATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCCCCATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGtttt                                                                                                            >  2:1048202/1‑149 (MQ=255)
                                          cAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTc                                                                                                           >  1:226317/1‑149 (MQ=255)
                                           aaaaGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCa                                                                                                          >  1:460617/1‑149 (MQ=255)
                                            aaaGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCaa                                                                                                         >  1:694723/1‑149 (MQ=255)
                                            aaaGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCaa                                                                                                         >  1:923868/1‑149 (MQ=255)
                                            aaaGGTTGGCAGACCTATGCCGTCGGTCTGGGGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCaa                                                                                                         <  1:515793/149‑1 (MQ=255)
                                             aaGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAAt                                                                                                        <  2:1504331/149‑1 (MQ=255)
                                               ggTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGc                                                                                                      >  1:1084541/1‑149 (MQ=255)
                                                     cAGACCTATGCCTGCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTAATCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATc                                                                                                <  1:809063/149‑1 (MQ=255)
                                                             tGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGt                                                                                        >  2:1561465/1‑149 (MQ=255)
                                                                cGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGGCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGAc                                                                                     <  2:1589882/149‑1 (MQ=255)
                                                                  tCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGc                                                                                   <  2:921784/149‑1 (MQ=255)
                                                                        tGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAgc                                                                             >  2:1157948/1‑149 (MQ=255)
                                                                         ggTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTGGATTATTCGGCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAgcg                                                                            <  2:625181/149‑1 (MQ=255)
                                                                          gTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAgcgc                                                                           <  1:784770/149‑1 (MQ=255)
                                                                          gTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAgcgc                                                                           >  1:268783/1‑149 (MQ=255)
                                                                            gACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCa                                                                         <  1:1136711/149‑1 (MQ=255)
                                                                               gTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGtt                                                                      <  1:378297/149‑1 (MQ=255)
                                                                                  ggTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTtcaa                                                                   <  1:789451/149‑1 (MQ=255)
                                                                                  ggTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTtcaa                                                                   <  1:432650/149‑1 (MQ=255)
                                                                                  ggTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTtcaa                                                                   <  2:680150/149‑1 (MQ=255)
                                                                                   tggATTAGTTTCTGGCTGGCTGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAATACGCCCAATAAGAATATTATCGGGCGTTAAAATATAACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTtcaat                                                                  <  1:239852/147‑1 (MQ=255)
                                                                                       ttAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCaaa                                                              <  1:1157948/149‑1 (MQ=255)
                                                                                        tAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTTGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAAc                                                             <  1:776329/149‑1 (MQ=255)
                                                                                        tAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAAc                                                             <  2:1181211/149‑1 (MQ=255)
                                                                                        tAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAAc                                                             >  1:478559/1‑149 (MQ=255)
                                                                                         aGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACg                                                            >  1:734475/1‑149 (MQ=255)
                                                                                            ttCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACgtgt                                                         <  2:268783/149‑1 (MQ=255)
                                                                                             tCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGtgtg                                                        >  1:381912/1‑149 (MQ=255)
                                                                                             tCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGtgtg                                                        >  1:387797/1‑149 (MQ=255)
                                                                                               tggctggcGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGAt                                                      <  2:1573524/149‑1 (MQ=255)
                                                                                                ggctggcGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATc                                                     <  2:785636/149‑1 (MQ=255)
                                                                                                     gcGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATc                                                <  2:986602/149‑1 (MQ=255)
                                                                                                     gcGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATc                                                >  2:1555977/1‑149 (MQ=255)
                                                                                                        gtgTTTATTATTCGTCGTAGAGATGAAGAAACTGAAATCGCACAATAAGAATATTTTCGGGCGTTAAAATATTTCATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGgt                                             <  1:198328/147‑1 (MQ=255)
                                                                                                        gggTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGgt                                             >  2:1423295/1‑149 (MQ=255)
                                                                                                        gggTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGgt                                             <  2:407118/149‑1 (MQ=255)
                                                                                                         ggTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATTCGTTATCCGCAGCGTGACGCTCCATCTCCAGTTCAATCAAACGTGTGATCAGATCGgtt                                            >  2:1130519/1‑148 (MQ=255)
                                                                                                           ttGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTa                                          >  2:1207725/1‑149 (MQ=255)
                                                                                                           ttGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTa                                          >  2:713649/1‑149 (MQ=255)
                                                                                                            tGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTaa                                         <  1:1087122/149‑1 (MQ=255)
                                                                                                            tGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTaa                                         <  2:942857/149‑1 (MQ=255)
                                                                                                              attattCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAAcc                                       <  2:505658/149‑1 (MQ=255)
                                                                                                                 attCTTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAAt                                                                                                        <  1:993866/81‑1 (MQ=255)
                                                                                                                 attCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAAt                                                                                                        >  2:758837/1‑81 (MQ=255)
                                                                                                                 attCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAAt                                                                                                        <  1:758837/81‑1 (MQ=255)
                                                                                                                      tcgtcgCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGc                               <  2:135442/149‑1 (MQ=255)
                                                                                                                           gcgATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCt                          <  1:128430/149‑1 (MQ=255)
                                                                                                                            cgATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCtt                         <  2:1097949/149‑1 (MQ=255)
                                                                                                                            cgATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCtt                         >  1:1588712/1‑149 (MQ=255)
                                                                                                                             gATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTtcaa                                                                   >  1:804786/1‑106 (MQ=255)
                                                                                                                             gATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTg                        >  1:219121/1‑149 (MQ=255)
                                                                                                                             gATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTg                        >  2:791819/1‑149 (MQ=255)
                                                                                                                                  agaaACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCcaca                   <  1:511063/149‑1 (MQ=255)
                                                                                                                                    aaaCTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCGATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGc                 >  2:948530/1‑149 (MQ=255)
                                                                                                                                      tattttaaCGCCCGATAATATTCTTATTGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGTTGTAACCCAGACCGCTGGCTTGCCACAGCtt               >  2:709033/7‑149 (MQ=38)
                                                                                                                                       cTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTc              <  2:694723/149‑1 (MQ=255)
                                                                                                                                       cTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTc              <  2:923868/149‑1 (MQ=255)
                                                                                                                                        agcgggcGCCCGATAATATTCTTATTGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGCGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCg             <  1:656290/143‑1 (MQ=255)
                                                                                                                                          aaaaCGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTgg                            >  2:799131/1‑132 (MQ=255)
                                                                                                                                          aaaaCGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGa           <  1:1595141/149‑1 (MQ=255)
                                                                                                                                          aaaaCGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGa           <  2:203082/149‑1 (MQ=255)
                                                                                                                                          aaaaCGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGa           <  2:346797/149‑1 (MQ=255)
                                                                                                                                             aCGCCCGATAATATTCTTATTGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGatac        >  1:973294/1‑149 (MQ=255)
                                                                                                                                              cGCCCGATAATATTCTTATTGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGataca       >  1:125067/1‑149 (MQ=255)
                                                                                                                                                  cAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACt   <  2:687298/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TATTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTG  >  NC_000913/399654‑399949

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