Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,165,428 C→T 5.2% P182P (CCG→CCA parC ← DNA topoisomerase IV, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,165,4280CT 5.2% 245.0 / 2.3 76P182P (CCG→CCAparCDNA topoisomerase IV, subunit A
Reads supporting (aligned to +/- strand):  ref base C (27/45);  new base T (2/2);  total (29/47)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.31e-01

GATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGC  >  NC_000913/3165281‑3165569
                                                                                                                                                   |                                                                                                                                             
gATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCGGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCg                                                                                                                                              <  7:188509/149‑1 (MQ=255)
gATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGTTTAATGCGATTGCCGCCGGAGCCACTTCACGCAGGTTATGCg                                                                                                                                              >  7:275031/1‑149 (MQ=255)
    ttACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGGTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTg                                                                                                                                           <  4:318807/148‑1 (MQ=255)
     tACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCCAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTgg                                                                                                                                          >  3:231027/1‑148 (MQ=255)
        ggATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCAGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAata                                                                                                                                      >  2:265804/1‑149 (MQ=255)
         gagctcGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGTAatat                                                                                                                                     <  7:23461/146‑1 (MQ=255)
         gATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAatat                                                                                                                                     <  8:142237/149‑1 (MQ=255)
         gATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTGTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTAACGCAGGTTATGCGGTGGAatat                                                                                                                                     >  6:317351/1‑149 (MQ=255)
               ggCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTc                                                                                                                                <  7:425816/148‑1 (MQ=255)
               ggCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCg                                                                                                                               >  4:21795/1‑149 (MQ=255)
                gcgcgcGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCACCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCGGGTCGATTAATGCGATTGCCGCCTGAGCCACATCACGCAGGTGATGCGGTGGAATATCGGTCGc                                                                                                                              >  5:240579/1‑149 (MQ=255)
                   cgcgAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGGCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCAt                                                                                                                           <  2:262415/149‑1 (MQ=255)
                    gcgAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGGGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATg                                                                                                                          >  7:343599/1‑149 (MQ=255)
                    gcgAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGTGGTGGAATATCGGTCGCCATg                                                                                                                          >  8:186067/1‑149 (MQ=255)
                    gcgAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATg                                                                                                                          <  2:381239/149‑1 (MQ=255)
                      gAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATgac                                                                                                                        >  2:155010/1‑147 (MQ=255)
                      gAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGcc                                                                                                                        >  2:10562/1‑149 (MQ=255)
                           gATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCAGACgg                                                                                                                   >  2:320464/1‑149 (MQ=255)
                                    cGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCg                                                                                                           <  7:286303/148‑1 (MQ=255)
                                     gCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCggt                                                                                                         <  8:343599/149‑1 (MQ=255)
                                       ttCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGc                                                                                                                         <  1:431797/131‑1 (MQ=255)
                                       ttCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGc                                                                                                                         >  2:431797/1‑131 (MQ=255)
                                          aGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGgtggtg                                                                                                     <  8:398975/148‑1 (MQ=255)
                                          aGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGc                                                                                                    <  3:95704/149‑1 (MQ=255)
                                            tCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCg                                                                                                  <  3:383179/149‑1 (MQ=255)
                                            tCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTAGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCg                                                                                                  >  3:326421/1‑149 (MQ=255)
                                               gATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGtt                                                                                                >  6:166294/1‑148 (MQ=255)
                                                aTAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTaa                                                                                              <  3:178803/149‑1 (MQ=255)
                                                aTAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATGGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTaa                                                                                              <  1:332989/149‑1 (MQ=255)
                                                   aTCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGc                                                                                            <  2:55200/148‑1 (MQ=255)
                                                    tCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGTGGTGGAATATCGGTCGCCAt                                                                                                                           >  7:119284/1‑116 (MQ=255)
                                                    tCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGTGGTGGAATATCGGTCGCCAt                                                                                                                           <  8:119284/116‑1 (MQ=255)
                                                    tCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAatat                                                                                                                                     >  5:208433/1‑106 (MQ=255)
                                                    tCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAatat                                                                                                                                     <  6:208433/106‑1 (MQ=255)
                                                    tCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCa                                                                                           >  3:97843/1‑148 (MQ=255)
                                                    tCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCa                                                                                           >  6:400916/1‑148 (MQ=255)
                                                     cctgCCCCTGCACGATTTCGAGCATCTGCTCGAGCGTGGGTTTCGGCTGGTCGATTAATCCGATTGCCGCCTGCGCCACTTCACGCAGGTTATGCGGTGGACTATCGGTCGCCCTGCCGACGGCAATACCGGTGGTGCCGTTAAGCaaa                                                                                         <  6:240579/146‑1 (MQ=255)
                                                      cGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCaaa                                                                                         <  4:231027/148‑1 (MQ=255)
                                                           ccTGCACGATATCCCGCAGCTGCTCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGttt                                                                                   <  8:395676/149‑1 (MQ=255)
                                                           ccTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGttt                                                                                   <  5:94853/149‑1 (MQ=255)
                                                           ccTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGttt                                                                                   <  5:166294/149‑1 (MQ=255)
                                                           ccTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGttt                                                                                   <  3:21795/149‑1 (MQ=255)
                                                             tGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGttt                                                                                   <  8:245388/147‑1 (MQ=255)
                                                             tGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGttt                                                                                   >  7:245388/1‑147 (MQ=255)
                                                               cACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGTGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTgg                                                                                 <  7:186067/147‑1 (MQ=255)
                                                                 cGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGa                                                                             >  4:248583/1‑149 (MQ=255)
                                                                   atatCCAGCAGCTGAGCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACg                                                                           <  8:350142/149‑1 (MQ=255)
                                                                     atCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGACGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACggg                                                                         <  4:82906/149‑1 (MQ=255)
                                                                          gcagcTGATCGAGCGTGGTTTTCGGCTGGGCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTTGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGt                                                                    <  5:42273/149‑1 (MQ=255)
                                                                           cagcTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTa                                                                   <  5:400916/149‑1 (MQ=255)
                                                                           cagcTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTa                                                                   <  4:147640/149‑1 (MQ=255)
                                                                              cTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGTAATATCGGTCGCCATGCCGACGGCAATACCGgtggtg                                                                                                     <  8:74660/112‑1 (MQ=255)
                                                                              cTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGgtggtg                                                                                                     >  7:74660/1‑112 (MQ=255)
                                                                              cTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCa                                                                <  3:368433/149‑1 (MQ=255)
                                                                                 aTCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGGGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCAttt                                                             >  8:379825/1‑149 (MQ=255)
                                                                                 aTCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCAt                                                               <  1:22703/147‑1 (MQ=255)
                                                                                    gAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGGGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTc                                                           <  6:306959/148‑1 (MQ=255)
                                                                                          ggTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAg                                                                      >  7:137777/1‑131 (MQ=255)
                                                                                          ggTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAg                                                                      <  8:137777/131‑1 (MQ=255)
                                                                                            ttttCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTg                                                   >  2:324503/1‑148 (MQ=255)
                                                                                                          aTTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGACAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCg                                      >  6:393795/1‑147 (MQ=255)
                                                                                                              aTGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGACAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGttt                                <  5:393795/149‑1 (MQ=255)
                                                                                                                 cGATTGGCGCCTGAGCCACTTCACGCAGGTTATGCTGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGc                             <  1:265804/149‑1 (MQ=255)
                                                                                                                 cGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTgg                              <  4:326421/148‑1 (MQ=255)
                                                                                                                 cGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTgg                              <  1:320464/148‑1 (MQ=255)
                                                                                                                 cGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGc                             <  4:97843/149‑1 (MQ=255)
                                                                                                                 cGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGc                             <  1:324503/149‑1 (MQ=255)
                                                                                                                   aTTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCAc                           <  1:155010/149‑1 (MQ=255)
                                                                                                                          ccTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGt                      <  6:123212/147‑1 (MQ=255)
                                                                                                                          ccTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGt                      >  5:123212/1‑147 (MQ=255)
                                                                                                                                 cACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTcc             >  2:103530/1‑149 (MQ=255)
                                                                                                                                   cTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTcccc           >  4:332787/1‑149 (MQ=255)
                                                                                                                                       aCGCAGGTTATGCGGCGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGc        <  3:166797/148‑1 (MQ=255)
                                                                                                                                            ggTTATGCGGGGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGGAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAgc  <  5:317351/149‑1 (MQ=255)
                                                                                                                                             gTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGc                             >  2:276879/1‑121 (MQ=255)
                                                                                                                                             gTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGc                             <  1:276879/121‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
GATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGC  >  NC_000913/3165281‑3165569

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