Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,162,333 A→G Y553H (TAT→CAT)  parC ← DNA topoisomerase IV subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,162,3330AG98.2% 339.3 / ‑3.2 109Y553H (TAT→CAT) parCDNA topoisomerase IV subunit A
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (49/58);  total (50/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 133 total aligned reads displayed.

CAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCA  >  NC_000913/3162243‑3162422
                                                                                           |                                                                                         
  aCGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGGCCACGAGCCGGCGGCAGGCTAATCGGGGTGATGGCATGGCtg                                                                                       >  1:1936748/1‑92 (MQ=255)
  aCGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTa                                                                                       <  2:3817978/93‑1 (MQ=255)
   cGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTAc                                                                                      <  2:1383398/93‑1 (MQ=255)
   cGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTAc                                                                                      <  1:3098493/93‑1 (MQ=255)
   cGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGGCAATGGCATGGCTAc                                                                                      >  1:1868986/1‑93 (MQ=255)
   cGGTCGCCCCAGGCGGCAACGGTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTAc                                                                                      <  2:1897547/93‑1 (MQ=255)
     gTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGGAATCGGGGCAATGGCATGGCTACGa                                                                                    >  2:3659418/1‑93 (MQ=255)
      gcGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGAc                                                                                   <  1:3704034/92‑1 (MQ=255)
       cGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGGAATCGGGGCAATGGCATGGCTACGAcc                                                                                  >  2:15583/1‑93 (MQ=255)
        gCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCg                                                                                 >  2:3370200/1‑93 (MQ=255)
          cccAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGt                                                                               <  1:3343453/93‑1 (MQ=255)
          cccAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGGAATCGGGTCAATGGCATGGCTACGACCggg                                                                               >  1:3282580/1‑92 (MQ=255)
           ccAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTg                                                                              <  2:1820236/93‑1 (MQ=255)
            cAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTgg                                                                             >  1:712228/1‑93 (MQ=255)
            cAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTgg                                                                             >  1:4326525/1‑93 (MQ=255)
            cAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCCGCCCACGCGCCGGCGGCCGCGTAATCGGGTCAATGGCATAGCTGCGACCGgggg                                                                             >  2:2855709/1‑93 (MQ=255)
             aGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGGCGACGGCCGCGCAATCGGGTCAATTGCAAGGCGACGACCCGCGGa                                                                            >  2:2476645/1‑93 (MQ=255)
                  gcAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGAACGGTGGGATcaa                                                                       >  1:4672211/1‑93 (MQ=255)
                  gcAACGTTAATTTGCCGGGGAGCGGCTCGACCGGGCGAAGAGCCCGGGGGACGCCGAATCGGTAAAAAGCATGGCTACGACCGGTGGAATtaa                                                                       >  2:4721838/1‑93 (MQ=255)
                   cAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCGCGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATcaac                                                                      <  2:2558480/93‑1 (MQ=255)
                    gaCGTTAATTTGCCTGTTGGCCGCTCGCCGTGGCCAGGCGGCGACGAGAACGTGATCGTGTCAAGGGCATTGCTATGACCGGTGGAATcaaca                                                                     <  2:4021973/92‑1 (MQ=255)
                    aaCGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGGTACGACCCGGGGGATcaaca                                                                     >  2:4288188/1‑93 (MQ=255)
                    aaCGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATcaaca                                                                     >  2:2307358/1‑93 (MQ=255)
                     aCGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCCGCGTAATCCGGTCAATGGCATGGCTACGACCCGTGGCATcaacaa                                                                    >  2:2586507/1‑93 (MQ=255)
                     aCGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCCACGGCAGCGCAAATCGGTCAATGGCACGGCTACGACCGGTGGAAGcaacaa                                                                    >  2:2412555/1‑93 (MQ=255)
                      cGTTAATTTGCCGGTGAGCGGCTCGCCCTTACCACGCGCCGACGGCAGCGCAAACGGGTCAATGGCATTGCTTCGACCGGGGGAATCaacaaa                                                                   >  1:3805415/1‑93 (MQ=255)
                      cGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCCGGGGAATCaacaaa                                                                   >  2:2214362/1‑93 (MQ=255)
                        ttAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGGCAATGGCATGGCTACCACCGGGGGAATCaaacaaca                                                                 >  1:859231/1‑93 (MQ=255)
                         tAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAacac                                                                <  2:1561829/93‑1 (MQ=255)
                          aaTTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACt                                                               <  2:1357842/93‑1 (MQ=255)
                            tttGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTAc                                                             <  2:2854203/93‑1 (MQ=255)
                              tGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCg                                                           <  1:3707439/93‑1 (MQ=255)
                                 cGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCggtt                                                        <  2:3271535/93‑1 (MQ=255)
                                  ggTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCggttg                                                       <  1:4210686/93‑1 (MQ=255)
                                  ggTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCggttg                                                       <  2:1819733/93‑1 (MQ=255)
                                  ggTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCggttg                                                       <  2:4234674/93‑1 (MQ=255)
                                     gAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCggttggtt                                                    <  1:181197/93‑1 (MQ=255)
                                     gAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCCTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCggttggtt                                                    <  2:1545764/93‑1 (MQ=255)
                                      aGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGgttggttg                                                   <  1:3610829/93‑1 (MQ=255)
                                      aGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGgttggttg                                                   <  1:2137944/93‑1 (MQ=255)
                                       gCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGGCAATGGCATGGCTACGACCGGGGGAATCAAAAAACACTACCCGTTGGTTGc                                                  >  1:3482436/1‑93 (MQ=255)
                                       gCGGCTCACCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGc                                                  >  2:1071317/1‑93 (MQ=255)
                                        cGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGct                                                 <  2:2463704/93‑1 (MQ=255)
                                          gCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGctct                                               <  2:4355385/93‑1 (MQ=255)
                                          gCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGctct                                               <  1:2650167/93‑1 (MQ=255)
                                          gCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGctct                                               <  2:1389769/93‑1 (MQ=255)
                                            tCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTa                                             >  2:1990450/1‑93 (MQ=255)
                                             cGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTAc                                            >  2:4407264/1‑93 (MQ=255)
                                              gCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTAcc                                           <  2:453800/93‑1 (MQ=255)
                                              gCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTAcc                                           >  2:2387392/1‑93 (MQ=255)
                                               cccTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCt                                          >  1:2212616/1‑93 (MQ=255)
                                               cccTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCt                                          >  2:1465905/1‑93 (MQ=255)
                                                ccTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCtt                                         >  1:227049/1‑93 (MQ=255)
                                                  tGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTc                                       >  1:1132629/1‑93 (MQ=255)
                                                  tGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTc                                       <  2:264878/93‑1 (MQ=255)
                                                   gACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTcc                                      >  1:1526606/1‑92 (MQ=255)
                                                      cACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCAccg                                   >  2:4435116/1‑93 (MQ=255)
                                                       aCGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCAccgc                                  >  2:3449296/1‑93 (MQ=255)
                                                       aCGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCAccgc                                  >  1:3246842/1‑93 (MQ=255)
                                                       aCGCGCCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCAccgc                                  >  2:2853370/1‑93 (MQ=255)
                                                        cgcgCCGACGGCAGCGAAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCAccgcc                                 <  1:4216202/93‑1 (MQ=255)
                                                         gcgcCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCAccgccg                                <  1:4021016/93‑1 (MQ=255)
                                                          cgcCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACcgccgc                               <  1:760909/93‑1 (MQ=255)
                                                          cgcCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACcgccgc                               <  1:4328068/93‑1 (MQ=255)
                                                           gcCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCt                              >  1:2365547/1‑93 (MQ=255)
                                                           gcCGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCt                              <  1:1438910/93‑1 (MQ=255)
                                                             ccGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCtt                             <  1:2899752/93‑1 (MQ=255)
                                                             ccGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCtt                             <  1:1228626/93‑1 (MQ=255)
                                                             ccGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCtt                             >  1:1316170/1‑93 (MQ=255)
                                                             ccGACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCtt                             >  1:1492924/1‑93 (MQ=255)
                                                               gACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTg                           >  2:1071185/1‑93 (MQ=255)
                                                               gACGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTg                           >  1:45885/1‑93 (MQ=255)
                                                                aCGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGa                          >  2:791434/1‑93 (MQ=255)
                                                                aCGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGa                          <  1:624499/93‑1 (MQ=255)
                                                                 cGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGaa                         <  2:4585744/93‑1 (MQ=255)
                                                                 cGGCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGCTGGTTGCTCTTACCTTTCACCGCCGCGTTGaa                         >  2:1247155/1‑93 (MQ=255)
                                                                   gCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGc                       >  1:2653194/1‑93 (MQ=255)
                                                                   gCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGc                       <  1:2815563/93‑1 (MQ=255)
                                                                   gCAGCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGc                       <  2:464638/93‑1 (MQ=255)
                                                                    cAGCGTAATCGGGTCAATGGCATGGCTACGACCGGGGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTCTAGCt                      >  2:2456050/1‑93 (MQ=255)
                                                                      gCGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTAt                    <  1:2199769/93‑1 (MQ=255)
                                                                       cGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATc                   >  2:2001797/1‑93 (MQ=255)
                                                                       cGTAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATc                   <  1:3702023/93‑1 (MQ=255)
                                                                         tAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCAc                 <  1:168/93‑1 (MQ=255)
                                                                         tAATCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCAc                 >  1:226935/1‑93 (MQ=255)
                                                                          aaTCGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCAcc                <  2:1094142/93‑1 (MQ=255)
                                                                             cGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGc             >  2:3858549/1‑93 (MQ=255)
                                                                             cGGGTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGc             <  2:443197/93‑1 (MQ=255)
                                                                              gggTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCt            <  2:1388597/93‑1 (MQ=255)
                                                                               ggTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCtt           >  1:4347966/1‑93 (MQ=255)
                                                                               ggTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCtt           <  2:3843929/93‑1 (MQ=255)
                                                                                gTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCttt          >  1:3854375/1‑93 (MQ=255)
                                                                                gTCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCttt          <  1:1174641/93‑1 (MQ=255)
                                                                                 tCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTta         >  2:2644778/1‑93 (MQ=255)
                                                                                 tCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTta         <  1:3864566/93‑1 (MQ=255)
                                                                                 tCAATGGCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTta         >  2:2202918/1‑93 (MQ=255)
                                                                                      ggCATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAAtt    >  2:197094/1‑93 (MQ=255)
                                                                                        cATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCa  <  1:3260250/93‑1 (MQ=255)
                                                                                        cATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCa  <  2:3874447/93‑1 (MQ=255)
                                                                                        cATGGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCa  <  2:599155/93‑1 (MQ=255)
                                                                                           |                                                                                         
CAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCA  >  NC_000913/3162243‑3162422

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