Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,165,499 C→T 6.1% A159T (GCC→ACC)  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,4990CT 6.1% 206.4 / 3.0 65A159T (GCC→ACC) parCDNA topoisomerase IV, subunit A
Reads supporting (aligned to +/- strand):  ref base C (28/33);  new base T (2/2);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

AGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGATTTCGGATCGTC  >  NC_000913/3165354‑3165646
                                                                                                                                                 |                                                                                                                                                   
agcagcTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAg                                                                                                                                                   <  4:196575/148‑1 (MQ=255)
 gcagcTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAgg                                                                                                                                                  >  4:412791/1‑148 (MQ=255)
         tCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAg                                                                                                                                                   >  3:3657/1‑139 (MQ=255)
         tCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAg                                                                                                                                                   <  4:3657/139‑1 (MQ=255)
           gAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCAtttt                                                                                                                                         <  8:52603/147‑1 (MQ=255)
            aGCGTTGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCTTTAATCAAAATTTTTGGCAGACGGGCAGGTAGCATTTTCg                                                                                                                                       <  8:320782/148‑1 (MQ=255)
                     ttCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGc                                                                                                                               <  7:174436/147‑1 (MQ=255)
                       cggcggGTCGATTAATGCGATTGCCGCCTGAGACACTTCACGCAGGTTATGCGGTGGATTATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCaa                                                                                                                             <  8:219992/147‑1 (MQ=255)
                       cGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCaa                                                                                                                             <  3:204690/147‑1 (MQ=255)
                         gCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGt                                                                                                                          >  8:231716/1‑148 (MQ=255)
                              tCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCAttt                                                                                                                                          <  1:380043/127‑1 (MQ=255)
                              tCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCAttt                                                                                                                                          >  2:380043/1‑127 (MQ=255)
                              tCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGATTGTCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGt                                                                                                                     >  2:96832/1‑148 (MQ=255)
                                     aTGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGtt                                                                                                              >  2:192644/1‑148 (MQ=255)
                                        cGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTg                                                                                                            <  1:192644/147‑1 (MQ=255)
                                         gATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTgg                                                                                                           >  7:305156/1‑147 (MQ=255)
                                           ttGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGTCGGGTAGGTAGCAttt                                                                                                                                          <  4:262437/114‑1 (MQ=255)
                                           ttGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGTCGGGTAGGTAGCAttt                                                                                                                                          <  4:256959/114‑1 (MQ=255)
                                           ttGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGTCGGGTAGGTAGCAttt                                                                                                                                          >  3:262437/1‑114 (MQ=255)
                                           ttGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGTCGGGTAGGTAGCAttt                                                                                                                                          >  3:256959/1‑114 (MQ=255)
                                            tGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCAcc                                                                                                       >  4:208443/1‑148 (MQ=255)
                                              ccgccTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCAccc                                                                                                      <  8:122852/147‑1 (MQ=255)
                                               cgccTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAg                                                                                                    >  2:42406/1‑148 (MQ=255)
                                                      gCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCg                                                                                             <  8:305156/148‑1 (MQ=255)
                                                          cTTCACGCCGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGt                                                                                            <  6:11836/145‑1 (MQ=255)
                                                          cTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGt                                                                                            >  5:11836/1‑145 (MQ=255)
                                                          cTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTcc                                                                                          <  6:192907/147‑1 (MQ=255)
                                                          cTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTcc                                                                                          >  5:192907/1‑147 (MQ=255)
                                                           ttCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCaa                                                                                                                             <  3:200481/111‑1 (MQ=255)
                                                           ttCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCaa                                                                                                                             >  4:200481/1‑111 (MQ=255)
                                                            tCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCa                                                                                                                                                    <  3:201135/87‑1 (MQ=255)
                                                            tCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCa                                                                                                                                                    >  4:201135/1‑87 (MQ=255)
                                                            tCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCt                                                                                       >  6:279982/1‑148 (MQ=255)
                                                            tCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATCTTTGGCAGACGGCCAGGTAGCAAATTCGGCTCCTGCAAAGAGCCGTCGAACTATGGCACCCAGTCACCCGTCCCCt                                                                                       >  4:380613/1‑148 (MQ=255)
                                                                   ggTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAg                                                                                <  1:42406/148‑1 (MQ=255)
                                                                    gTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCa                                                                                                 <  2:6799/130‑1 (MQ=255)
                                                                    gTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCa                                                                                                 >  1:6799/1‑130 (MQ=255)
                                                                     ttATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCaa                                                                                                                             <  8:389682/101‑1 (MQ=255)
                                                                     ttATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCaa                                                                                                                             >  7:389682/1‑101 (MQ=255)
                                                                         gCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAgctcgc                                                                           <  3:208443/147‑1 (MQ=255)
                                                                         gCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGCCTCCTGCAAAGTGCCGTCGAAGTGTGGCACCCAGTAATCCGTCCCCTGCCCCAgctcgct                                                                          >  2:126641/1‑148 (MQ=255)
                                                                            gTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCCCCCAGTCAGCCGTCCCCTGCCCCAGCTCGCGCaa                                                                       >  8:158128/1‑148 (MQ=255)
                                                                            gTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCaa                                                                       >  4:256071/1‑148 (MQ=255)
                                                                               gAATATCGGTCGCCATGCCGACGGCAATACCGGTGGGGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATag                                                                    <  5:251252/148‑1 (MQ=255)
                                                                               gAATATAGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGAAGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCTTCCCCTGCCCCAGCTCGCTCAATag                                                                    <  1:126641/148‑1 (MQ=255)
                                                                                   atCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCAccc                                                                                                      <  8:228642/110‑1 (MQ=255)
                                                                                   atCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCAccc                                                                                                      >  7:228642/1‑110 (MQ=255)
                                                                                   atCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATagc                                                                   >  1:109477/1‑145 (MQ=255)
                                                                                   atCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATagc                                                                   <  2:109477/145‑1 (MQ=255)
                                                                                    tCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGGCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATagcagc                                                                <  1:96832/147‑1 (MQ=255)
                                                                                        tCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGa                                                           <  2:30846/148‑1 (MQ=255)
                                                                                           ccATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAata                                                        <  5:279982/148‑1 (MQ=255)
                                                                                             aTGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATAtt                                                      <  3:412791/148‑1 (MQ=255)
                                                                                               gCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGtt                                                                                                              >  2:252808/1‑90 (MQ=255)
                                                                                               gCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGtt                                                                                                              <  1:252808/90‑1 (MQ=255)
                                                                                                 cGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGa                                                  >  1:374197/1‑148 (MQ=255)
                                                                                                          tACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGa                                         >  1:262403/1‑148 (MQ=255)
                                                                                                           aCCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGt                                                                                                               <  7:22874/77‑1 (MQ=255)
                                                                                                           aCCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGt                                                                                                               >  8:22874/1‑77 (MQ=255)
                                                                                                                        ttAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCa                           <  3:380613/148‑1 (MQ=255)
                                                                                                                        ttAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCa                           <  2:262403/148‑1 (MQ=255)
                                                                                                                              aaaaTGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAAc                                              >  5:135034/1‑123 (MQ=255)
                                                                                                                              aaaaTGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAAc                                              <  6:135034/123‑1 (MQ=255)
                                                                                                                                        gCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGAtt            <  7:231716/147‑1 (MQ=255)
                                                                                                                                           gACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATACCAGCTCGCAATATTACGCCAACCGGGATTCGGAGTAACGCATTGACGCCAACCCTTTCg         >  4:259829/1‑147 (MQ=255)
                                                                                                                                                 cAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGATTTCGGATCGTc  <  2:374197/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
AGCAGCTGATCGAGCGTGGTTTTCGGCTGGTCGATTAATGCGATTGCCGCCTGAGCCACTTCACGCAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGATTTCGGATCGTC  >  NC_000913/3165354‑3165646

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: