Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,267,060 G→A S2S (AGC→AGT)  tdcA ← DNA‑binding transcriptional activator TdcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,267,0600GA99.8% 1778.9 / NA 416S2S (AGC→AGT) tdcADNA‑binding transcriptional activator TdcA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (231/184);  minor base C (0/1);  total (231/185)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 416 total aligned reads displayed.

TGAGTTAACCCTAATTCTTTTGCAGCCGAGCCGATAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTGCTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCACTCACAAATTAACTTTCAA  >  NC_000913/3266962‑3267157
                                                                                                  |                                                                                                 
 gAGTTAACCCTAATTCTTTTGCAGCCGAGCCGATAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTc                                                                                                >  2:6057955/1‑101 (MQ=255)
  aGTTAACCCTAATTCTTTTGCAGCCGAGCCGATAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCa                                                                                               >  1:6057957/1‑101 (MQ=255)
   gTTAACCCTAATTCTTTTGCAGCCGAGCCGATAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCAt                                                                                              >  1:6057964/1‑101 (MQ=255)
   gTTAACCCTAATTCTTTTGCAGCCGAGCCGATAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCAt                                                                                              <  1:6057262/101‑1 (MQ=255)
   gTTAACCCTAATTCTTTTGCAGCCGAGCCGATAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCAt                                                                                              >  2:6057963/1‑101 (MQ=255)
      aaCCCTAATTCTTTTGCAGCCGAGCCGATAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAAc                                                                                           >  1:6057970/1‑101 (MQ=255)
       aCCCTAATTCTTTTGCAGCCGAGCCGATAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACg                                                                                          >  2:6057975/1‑101 (MQ=255)
         ccTAATTCTTTTGCAGCCGAGCCGATAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGAc                                                                                        >  1:6057979/1‑101 (MQ=255)
            aaTTCTTTTGCAGCCGAGCCGATAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTa                                                                                     >  2:6057984/1‑101 (MQ=255)
             aTTCTTTTGCAGCCGAGCCGATAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTAc                                                                                    >  1:6057986/1‑101 (MQ=255)
                cTTTTGCAGCCGAGCCGATAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGtt                                                                                 >  2:6057994/1‑101 (MQ=255)
                 ttttGCAGCCGAGCCGATAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTa                                                                                >  1:6057997/1‑101 (MQ=255)
                 ttgtGCAGCCGAGCCGATAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTAc                                                                                    <  1:6057212/94‑1 (MQ=255)
                    tGCAGCCGAGCCGATAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAAtt                                                                             >  2:6058000/1‑101 (MQ=255)
                      cAGCCGAGCCGATAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTAc                                                                           >  1:6058005/1‑101 (MQ=255)
                        gccgagccgaTAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCt                                                                                      <  1:6057338/88‑1 (MQ=255)
                         ccgagccgaTAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTc                                                                        >  2:6058010/1‑101 (MQ=255)
                          ggagccgaTAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCa                                                                       <  2:6057219/100‑1 (MQ=255)
                           gagccgaTAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATa                                                                                             <  2:6057573/78‑1 (MQ=255)
                           gagccgaTAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGtt                                                                                 <  1:6057440/90‑1 (MQ=255)
                           gagccgaTAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGtt                                                                                 <  2:6057253/90‑1 (MQ=255)
                           gagccgaTAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTaa                                                                               <  2:6057290/92‑1 (MQ=255)
                           gagccgaTAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTAc                                                                           <  1:6057241/96‑1 (MQ=255)
                           gagccgaTAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTc                                                                        <  2:6057394/99‑1 (MQ=255)
                           gagccgaTAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCAt                                                                      >  1:6058017/1‑101 (MQ=255)
                           gagccgaTAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCAt                                                                      >  2:6058014/1‑101 (MQ=255)
                            agccgaTAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATaa                                                                                            <  1:6057473/78‑1 (MQ=255)
                            agccgaTAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGt                                                                                  <  2:6057422/88‑1 (MQ=255)
                            agccgaTAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTAc                                                                           <  2:6057235/95‑1 (MQ=255)
                            agccgaTAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCAtt                                                                     <  2:6057469/101‑1 (MQ=255)
                            agccgaTAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCAtt                                                                     <  1:6057421/101‑1 (MQ=255)
                              ccgaTAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGa                                                                   <  2:6057354/101‑1 (MQ=255)
                               cgaTAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGAc                                                                  <  1:6057515/101‑1 (MQ=255)
                               cgaTAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGAc                                                                  <  2:6057114/101‑1 (MQ=255)
                               cgaTAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGAc                                                                  >  1:6058022/1‑101 (MQ=255)
                                gaTAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGa                                                                   >  1:6058026/1‑99 (MQ=255)
                                  tAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGc                                                               >  2:6058031/1‑101 (MQ=255)
                                    gAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCAt                                                             >  1:6058032/1‑101 (MQ=255)
                                       ccACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGaa                                                          <  2:6057070/101‑1 (MQ=255)
                                       ccACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGaa                                                          >  2:6058039/1‑101 (MQ=255)
                                       ccACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGaa                                                          <  1:6057362/101‑1 (MQ=255)
                                        cACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAg                                                         >  1:6058040/1‑101 (MQ=255)
                                           ttCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTg                                                                    >  1:6058044/1‑87 (MQ=255)
                                           ttCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAgtgt                                                      <  2:6057336/101‑1 (MQ=255)
                                           ttCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAgtgt                                                      <  2:6056756/101‑1 (MQ=255)
                                            tCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTa                                                     <  2:6057495/101‑1 (MQ=255)
                                                aTGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCaa                                                 >  1:6058052/1‑101 (MQ=255)
                                                 tGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCaaa                                                <  2:6057627/101‑1 (MQ=255)
                                                  gACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCaaaa                                               >  2:6058057/1‑101 (MQ=255)
                                                   aCTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAaaat                                              >  1:6058059/1‑101 (MQ=255)
                                                      tCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAaaatgaa                                           >  2:6058063/1‑101 (MQ=255)
                                                      tCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAaaatgaa                                           <  1:6057294/101‑1 (MQ=255)
                                                       ccTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGa                                                                                         <  1:6057014/54‑1 (MQ=255)
                                                        cTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAaaatgaaat                                         <  2:6057313/101‑1 (MQ=255)
                                                         tGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAaaatgaaatg                                        >  2:6058069/1‑101 (MQ=255)
                                                          gAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAaatgaaatga                                       >  2:6058073/1‑101 (MQ=255)
                                                          gAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAaatgaaatga                                       >  1:6058075/1‑101 (MQ=255)
                                                          gAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAaatgaaatga                                       <  2:6057509/101‑1 (MQ=255)
                                                              gACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAg                                   >  2:6058081/1‑101 (MQ=255)
                                                              gACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAg                                   <  1:6057604/101‑1 (MQ=255)
                                                               aCTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAgg                                  >  1:6058083/1‑101 (MQ=255)
                                                               aCTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAgg                                  <  2:6057441/101‑1 (MQ=255)
                                                                cTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGa                                 <  2:6057596/101‑1 (MQ=255)
                                                                cTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGa                                 <  2:6057556/101‑1 (MQ=255)
                                                                 tACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGat                                >  1:6058087/1‑101 (MQ=255)
                                                                   ccAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGatat                              <  1:6057389/101‑1 (MQ=255)
                                                                   ccAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGatat                              >  2:6058091/1‑101 (MQ=255)
                                                                   ccAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGatat                              <  2:6057021/101‑1 (MQ=255)
                                                                    cAGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATAtg                             >  2:6058098/1‑101 (MQ=255)
                                                                     aGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATAtgt                            >  1:6058102/1‑101 (MQ=255)
                                                                     aGGTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATAtgt                            >  2:6058103/1‑101 (MQ=255)
                                                                      ggTGCTGCGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATAtgtg                           <  2:6057175/101‑1 (MQ=255)
                                                                        tgctgcGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCg                         >  2:6058110/1‑101 (MQ=255)
                                                                         gctgcGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGa                        >  1:6058114/1‑101 (MQ=255)
                                                                          ctgcGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGa                                 <  1:6057102/91‑1 (MQ=255)
                                                                          ctgcGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGAc                       >  2:6058116/1‑101 (MQ=255)
                                                                           tgcGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGAcc                      >  2:6058121/1‑101 (MQ=255)
                                                                            gcGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCa                     >  1:6058124/1‑101 (MQ=255)
                                                                            gcGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCa                     <  1:6056639/101‑1 (MQ=255)
                                                                             cGTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCAc                    >  2:6058128/1‑101 (MQ=255)
                                                                              gTTTTCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCACt                   >  1:6058136/1‑101 (MQ=255)
                                                                               ttttCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCACTc                  <  1:6057163/101‑1 (MQ=255)
                                                                               ttttCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCACTc                  >  2:6058137/1‑101 (MQ=255)
                                                                                tttCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCACTca                 <  2:6057370/101‑1 (MQ=255)
                                                                                 ttCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCACTcac                >  2:6058143/1‑101 (MQ=255)
                                                                                 ttCGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCACTcac                <  1:6057802/101‑1 (MQ=255)
                                                                                   cGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCACTCACaa              >  1:6058148/1‑101 (MQ=255)
                                                                                   cGGAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCACTCACaa              >  2:6058149/1‑101 (MQ=255)
                                                                                    ggAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCa                                                  <  2:6056940/64‑1 (MQ=255)
                                                                                    ggAAGAAGAATAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCa                     <  1:6057698/93‑1 (MQ=255)
                                                                                      aagaagaaTAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCACTCACAAAtt           >  1:6058155/1‑101 (MQ=255)
                                                                                        gaagaaTAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCACTCACAAATTaa         <  1:6057401/101‑1 (MQ=255)
                                                                                        gaagaaTAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCACTCACAAATTaa         >  1:6058160/1‑101 (MQ=255)
                                                                                         aagaaTAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCACTCACAAATTAAc        <  2:6057488/101‑1 (MQ=255)
                                                                                         aagaaTAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCACTCACAAATTAAc        >  2:6058161/1‑101 (MQ=255)
                                                                                          agaaTAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCACTCACAAATTAACt       <  2:6057840/101‑1 (MQ=255)
                                                                                           gaaTAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCACTCACAAATTAACtt      >  2:6058167/1‑101 (MQ=255)
                                                                                            aaTAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCACTCACAAATTAACttt     >  1:6058171/1‑101 (MQ=255)
                                                                                            aaTAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCACTCACAAATTAACttt     >  2:6058176/1‑101 (MQ=255)
                                                                                             aTAGTACTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCACTCACAAATTAACTTTc    >  1:6058180/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
TGAGTTAACCCTAATTCTTTTGCAGCCGAGCCGATAGAACCACTTCTAATGACTTCCTGAAAGACTACCAGGTGCTGCGTTTTCGGAAGAAGAATAGTGCTCATAACGACCTACGTTAATTACCTCATTGACGGCATGAAGTGTATCAAAATGAAATGAACAGGATATGTGCGACCACTCACAAATTAACTTTCAA  >  NC_000913/3266962‑3267157

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