Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,947,332 A→T 15.3% A212A (GCT→GCA hdfR ← flhDC operon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,947,3320AT15.3% 455.7 / 71.1 182A212A (GCT→GCAhdfRflhDC operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (85/69);  new base T (17/11);  total (102/80)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

Only 100 of 211 total aligned reads displayed.

GCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGA  >  NC_000913/3947233‑3947431
                                                                                                   |                                                                                                   
gCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTca                                                                                                     >  1:3175836/1‑100 (MQ=255)
gCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTca                                                                                                     <  1:3643051/100‑1 (MQ=255)
gCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTca                                                                                                     >  2:3068519/1‑100 (MQ=255)
gCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTca                                                                                                     <  1:458347/100‑1 (MQ=255)
gCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTca                                                                                                     >  1:488144/1‑100 (MQ=255)
 cTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTcag                                                                                                    <  1:2906974/100‑1 (MQ=255)
 cTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAAACATTAAGCATCGCAATCTGTTGCTGTGCCAGTTcag                                                                                                    >  2:3641527/1‑100 (MQ=255)
      gACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGaa                                                                                               >  1:1075775/1‑100 (MQ=255)
      gACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGaa                                                                                               >  1:3128848/1‑100 (MQ=255)
        cAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAAc                                                                                              <  2:557372/99‑1 (MQ=255)
        cAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACt                                                                                             >  2:505846/1‑100 (MQ=255)
          aCGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGGCGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTgg                                                                                           >  1:2684798/1‑100 (MQ=255)
           cGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCTGCTGAACTggc                                                                                          <  1:1513053/100‑2 (MQ=255)
           cGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCTGCTGAACTggc                                                                                          <  1:1377128/100‑2 (MQ=255)
           cGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGt                                                                                          <  1:3585516/100‑1 (MQ=255)
            ggTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTg                                                                                            >  2:3739244/1‑97 (MQ=255)
             gTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCTGCTGAACTggcac                                                                                        <  1:1871685/100‑4 (MQ=255)
              tATGCAGGCCGCCTTTTTTACGCGCCCAGCTGATGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCTGCTGAACTggcaca                                                                                       <  2:2179085/100‑5 (MQ=255)
                tGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTACCAGTTCTGCTGAACTggcacagc                                                                                     <  1:574283/100‑7 (MQ=255)
                 gCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAg                                                                                    >  1:1792121/1‑100 (MQ=255)
                 gCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAg                                                                                    >  1:1023691/1‑100 (MQ=255)
                  cAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGa                                                                                   >  2:1940503/1‑100 (MQ=255)
                  cAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGa                                                                                   >  1:2457624/1‑100 (MQ=255)
                  cAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGATGTGCCAGTTCAGCAGAACTGGTTGTCAGa                                                                                   >  2:2403104/1‑100 (MQ=255)
                  cAGGCCGCCTTTTTTACGCGCCCACCTGACGGGTAGCCAGGTGCAACCATTAACCATCGCCCTCTGCTGCTGTGCCCGTTCAGCAGAACTGGTTGTCAGa                                                                                   >  2:628575/1‑100 (MQ=255)
                    ggCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCGGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGaa                                                                                  >  2:2362144/1‑99 (MQ=255)
                    ggCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAAt                                                                                 <  2:1681446/100‑1 (MQ=255)
                      ccgccTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATgg                                                                               <  1:1974303/100‑1 (MQ=255)
                       cgccTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATggg                                                                              >  2:292911/1‑100 (MQ=255)
                         gcctttttACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCa                                                                            <  1:3499953/97‑1 (MQ=255)
                         ccTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCa                                                                            >  1:2009230/1‑100 (MQ=255)
                                 aCGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCa                                                                    <  2:1737921/100‑1 (MQ=255)
                                  cgcgCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAg                                                                   >  1:911810/1‑100 (MQ=255)
                                  cgcgCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAg                                                                   >  1:733188/1‑100 (MQ=255)
                                  cgcgCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAg                                                                   <  1:3009798/100‑1 (MQ=255)
                                   gcgcCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGc                                                                  >  2:3261904/1‑100 (MQ=255)
                                    cgcCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCa                                                                 <  2:1703691/100‑1 (MQ=255)
                                        cAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGa                                                             <  2:3807457/100‑1 (MQ=255)
                                          gCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATc                                                           <  1:2493363/100‑1 (MQ=255)
                                            tGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCan                                                         <  2:1900398/100‑2 (MQ=255)
                                             gACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCaaa                                                        >  1:1349985/1‑100 (MQ=255)
                                               cGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAAcc                                                      <  1:67360/100‑1 (MQ=255)
                                                 ggTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTg                                                    >  1:2282790/1‑100 (MQ=255)
                                                   tAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGcc                                                  <  1:1830448/100‑1 (MQ=255)
                                                    aGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCt                                                 <  2:3539557/100‑1 (MQ=255)
                                                     gCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTc                                                >  1:3293445/1‑100 (MQ=255)
                                                       cAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCAt                                              >  2:1640070/1‑100 (MQ=255)
                                                         ggTGCAACCATTCCGCATCGTAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAACCCTGCCTCATGc                                            <  2:3605666/100‑1 (MQ=255)
                                                          gTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCt                                           >  1:1000926/1‑100 (MQ=255)
                                                             cAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGtt                                        >  2:3048465/1‑100 (MQ=255)
                                                             cAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGtt                                        <  1:1236670/100‑1 (MQ=255)
                                                              aaCCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTg                                       >  1:331008/1‑100 (MQ=255)
                                                                  aTTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGaaaa                                   >  2:279033/1‑100 (MQ=255)
                                                                   ttAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAAt                                  >  1:2345354/1‑100 (MQ=255)
                                                                    tAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATc                                 <  1:194666/100‑1 (MQ=255)
                                                                    tAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATc                                 <  2:1059522/100‑1 (MQ=255)
                                                                     aaGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCt                                <  2:3335066/100‑1 (MQ=255)
                                                                      aGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTg                               <  2:1345359/100‑1 (MQ=255)
                                                                       gCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTgg                              >  2:852382/1‑100 (MQ=255)
                                                                        cATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGc                             <  1:3146518/100‑1 (MQ=255)
                                                                        cATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGc                             <  2:2914181/100‑1 (MQ=255)
                                                                         aTCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATTCTGTTGAAAATGTgg                              >  2:269079/1‑98 (MQ=255)
                                                                         aTCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGcc                            <  2:208959/100‑1 (MQ=255)
                                                                         aTCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGcc                            <  2:1517344/100‑1 (MQ=255)
                                                                          tCGCAATCTGTTGCGGTACCAGTTCAGCAGAGCTTGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCCGCCTCACGCTGATGAAAATcgcgcgc                           >  2:2424021/1‑94 (MQ=255)
                                                                           cGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGcccc                          >  1:1439610/1‑100 (MQ=255)
                                                                            gCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGccccc                         <  1:2721295/100‑1 (MQ=255)
                                                                             cAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAGTGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTggagatcg                        >  1:3559770/1‑94 (MQ=255)
                                                                               aTCTGTTGCTGTGCCAGTTCAGCAGATCTGGTTGTCAGAATGGGCACTTCGTCCGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCAc                       <  2:1240746/99‑1 (MQ=255)
                                                                               aTCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACt                      >  1:1202081/1‑100 (MQ=255)
                                                                               aTCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACt                      <  2:331006/100‑1 (MQ=255)
                                                                                tCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTc                     >  2:320026/1‑100 (MQ=255)
                                                                                   gTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAg                  <  2:1096108/100‑1 (MQ=255)
                                                                                   gTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAg                  <  1:309235/100‑1 (MQ=255)
                                                                                    ttGCTGTGCCAGTTCAGCAGCACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGt                 >  2:1816534/1‑100 (MQ=255)
                                                                                       gacaaccAGTTCTGCTGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGGCTCATACTGTTGAAGATCTGGCCCCCACTCAAGTCGc              >  2:2575822/6‑100 (MQ=255)
                                                                                       gacaaccAGTTCTGCTGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAACATCTGGCCCCCACTCAAGTCGc              >  1:3060291/6‑100 (MQ=255)
                                                                                       cTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGc              >  1:3392265/1‑100 (MQ=255)
                                                                                        tgtgCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCa             <  2:1320018/100‑1 (MQ=255)
                                                                                        tctgCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCa             >  1:2580094/3‑100 (MQ=255)
                                                                                         gtgCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAg            >  2:397466/1‑100 (MQ=255)
                                                                                         gtgCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAg            <  2:2213114/100‑1 (MQ=255)
                                                                                         caaccAGTTCTGCTGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAg            >  1:1294425/4‑100 (MQ=255)
                                                                                         caaccAGTTCTGCTGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAg            >  1:1668872/4‑100 (MQ=255)
                                                                                         caaccAGTTCTGCTGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAg            >  1:2825965/4‑100 (MQ=255)
                                                                                         caaccAGTTCTGCTGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAg            >  1:536159/4‑100 (MQ=255)
                                                                                         caaccAGTTCTGCTGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAg            >  1:2150565/4‑100 (MQ=255)
                                                                                          tgCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCTCTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGa           <  2:2010337/100‑1 (MQ=255)
                                                                                          tgCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGa           >  2:2900288/1‑100 (MQ=255)
                                                                                          aaccAGTTCTGCTGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGa           >  1:3244521/3‑100 (MQ=255)
                                                                                           accAGTTCTGCTGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGa           >  1:2649792/2‑99 (MQ=255)
                                                                                             cAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATaa        >  2:244574/1‑100 (MQ=255)
                                                                                               gTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAAtt      >  1:397215/1‑100 (MQ=255)
                                                                                                  cagcagAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAg   >  2:3491634/1‑100 (MQ=255)
                                                                                                  cagcagAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAg   >  1:252741/1‑100 (MQ=255)
                                                                                                  cagcagAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAg   <  2:2544985/100‑1 (MQ=255)
                                                                                                  cagcagAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAg   <  1:3484439/100‑1 (MQ=255)
                                                                                                   agcagAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGa  >  2:1409896/1‑100 (MQ=255)
                                                                                                   agcagAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGa  >  1:1406691/1‑100 (MQ=255)
                                                                                                   agcagAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGa  >  2:2495836/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGA  >  NC_000913/3947233‑3947431

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