Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,752,107:1 +G 15.4% coding (116/975 nt) yiaU → putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,752,1071.G15.4% 585.5 / 116.4 175coding (116/975 nt)yiaUputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (68/80);  new base G (9/18);  total (77/98)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

TTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTA  >  NC_000913/3752008‑3752206
                                                                                                    |                                                                                                   
ttAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCaa                                                                                                      <  2:3583811/100‑1 (MQ=255)
   aGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAata                                                                                                  <  1:3948357/100‑1 (MQ=255)
    gTACCGGGAGTTAAAAATTATCTCGGTGATCGCTGCCAGTGCAAATAACAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCCACCTCCGCAActat                                                                                                 >  1:1084796/1‑100 (MQ=255)
      acagggAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAata                                                                                                  <  1:1424998/94‑1 (MQ=255)
      aCCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCt                                                                                               <  2:887850/100‑1 (MQ=255)
        cGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTg                                                                                             >  2:4054298/1‑100 (MQ=255)
           gAGTTAAAAATTATCTCGGTAATCTCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGc                                                                                            <  1:481080/98‑1 (MQ=255)
           gAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTg                                                                                          <  1:3372593/100‑1 (MQ=255)
            aGTTAAAAATTATCTCGGTAATCGCTGCGTGTGAAAATAACAGCCATGCCGCGACGGTACTTGGCATCGCAAAGGCCAACGTCAGCAAATATCTTGc                                                                                            <  1:1853843/97‑1 (MQ=255)
             gTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAGATATTTGCTGac                                                                                        <  2:4296305/100‑2 (MQ=255)
             gTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAGATATTTGCTGac                                                                                        <  2:4316244/100‑2 (MQ=255)
             gTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGACCAACGTCAGCAAATATCTTGCTGAt                                                                                        >  1:2184105/1‑100 (MQ=255)
               tAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAGATATTTGCTGacgt                                                                                      <  2:4074396/100‑4 (MQ=255)
               tAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGAttt                                                                                      >  2:3660999/1‑100 (MQ=255)
                aaaaaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAGATATTTGCTGacgtt                                                                                     <  1:1281205/100‑5 (MQ=255)
                aaaaaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGAtttt                                                                                     <  2:2690403/100‑1 (MQ=255)
                 aaatTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAGATATTTGCTGacgttg                                                                                    <  1:4230385/100‑6 (MQ=255)
                 aaaaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTg                                                                                    <  1:2617046/100‑1 (MQ=255)
                 aaaaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTg                                                                                    <  2:2263658/100‑1 (MQ=255)
                 aaaaTTATCTCGGTAATCGCTGCCAGTGAAAATATAAGCAATGCCGCGACTGTACTTGGCATCGCACAGTCCAACGTCAGCAAGATATTTGCTGacgttg                                                                                    <  1:3563355/100‑6 (MQ=255)
                  aaaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAGATATTTGCTGacgttgg                                                                                   <  2:2334328/100‑7 (MQ=255)
                  aaaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAGATATTTGCTGacgttgg                                                                                   <  2:3393473/100‑7 (MQ=255)
                   aaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAGATATTTGCTGacgttggc                                                                                  <  1:3576582/100‑8 (MQ=38)
                   aaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGaa                                                                                  <  2:1058507/100‑1 (MQ=255)
                    aTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTATTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAAt                                                                                 <  1:432462/100‑1 (MQ=255)
                    aTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAGATATTTGCTGacgttggcc                                                                                 <  2:1857516/100‑9 (MQ=38)
                    aTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAAt                                                                                 <  2:3754602/100‑1 (MQ=255)
                      tATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACATCAGCAAGATATTTGCTGacgttggcctg                                                                               <  2:2126332/100‑11 (MQ=37)
                       aTCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCaa                                                                              <  1:170016/100‑1 (MQ=255)
                          tcGGTAATCGCTGCCAGTTAAAATATCAGCCATTCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAg                                                                           <  1:3445501/100‑1 (MQ=255)
                           cGGTAATCGCTGCCAGTTAAAATAACAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAg                                                                           >  1:3125963/1‑99 (MQ=255)
                           cGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGt                                                                          >  1:2900706/1‑100 (MQ=255)
                             ttAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTgg                                                                        <  1:2289759/99‑1 (MQ=255)
                               aaTCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGt                                                                      <  1:3180548/100‑1 (MQ=255)
                               aaTCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGt                                                                      <  1:3155145/100‑1 (MQ=255)
                               aaTCACAGCCAGTGAATATATCAGTCACGCCGCAACTGTACTTGGTTTCGGATAGGCCAACGTCAGCAAGTATCTTGCTGCTTTTGAATCAAACGCGGGt                                                                      <  1:2802295/100‑1 (MQ=255)
                                aTCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGACTCAAAAGTggggt                                                                     <  2:2983574/100‑3 (MQ=255)
                                   gCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGttt                                                                    <  1:4192133/98‑1 (MQ=255)
                                   gCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCGTGCTGATTTTGAATCAAAAGTGGGTTTaa                                                                  >  1:3458546/1‑100 (MQ=255)
                                    cTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTaaa                                                                 >  2:3083590/1‑100 (MQ=255)
                                    cTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTaaa                                                                 >  1:4099759/1‑100 (MQ=255)
                                    cTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTaaa                                                                 <  2:3311681/100‑1 (MQ=255)
                                     tGCCAGTGTAAATATTTGCTATGCTGCGATTGTACTTGGTATTGAAGAGGAAATCGTCTGCAAATATCTTGCTGATTTTGTATCAAAAGTGGGTTTa                                                                   <  1:3083893/97‑1 (MQ=255)
                                      gCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAg                                                               >  2:4240713/1‑100 (MQ=255)
                                        cAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTc                                                             >  2:2329492/1‑100 (MQ=255)
                                         aGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTGTGAATCAAAAGTGGGTTTAAAAGTCt                                                            >  2:2187282/1‑100 (MQ=255)
                                          gTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCtt                                                           >  2:182026/1‑100 (MQ=255)
                                           tGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGTTTTTGAATCAAAAGTGGGTTTAAAAGTCtt                                                           <  1:2891888/99‑1 (MQ=255)
                                           tGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCttt                                                          <  2:2749662/100‑1 (MQ=255)
                                            gAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTg                                                         <  1:4091087/100‑1 (MQ=255)
                                            gAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTg                                                         <  2:155108/100‑1 (MQ=255)
                                            gAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTg                                                         >  1:680602/1‑100 (MQ=255)
                                            gAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTg                                                         <  2:666000/100‑1 (MQ=255)
                                               aaTATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGAcc                                                      <  2:3255597/100‑1 (MQ=255)
                                                 tatCAGCCATGCCGCGACTGTACTTGGCATCGCACATGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGa                                                    >  2:3099311/1‑100 (MQ=255)
                                                   tCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGaan                                                  <  2:3874436/100‑2 (MQ=255)
                                                    cAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAacc                                                 <  1:2230303/100‑1 (MQ=255)
                                                    cAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAacc                                                 <  2:4382157/100‑1 (MQ=255)
                                                     aGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAacca                                                >  2:857235/1‑100 (MQ=255)
                                                     aGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAacca                                                >  1:3003774/1‑100 (MQ=255)
                                                      gCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAaccac                                               >  2:83485/1‑100 (MQ=255)
                                                        cATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCAccc                                             >  1:3847171/1‑100 (MQ=255)
                                                        cATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCAccc                                             >  2:130136/1‑100 (MQ=255)
                                                             cgcgACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCcagcag                                        >  1:886386/1‑100 (MQ=255)
                                                              gcgACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGc                                       >  2:4291398/1‑100 (MQ=255)
                                                              gcgACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGc                                       <  2:3803992/100‑1 (MQ=255)
                                                              gcgACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGc                                       >  2:3736908/1‑100 (MQ=255)
                                                                   tGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATg                                  <  2:3907694/100‑1 (MQ=255)
                                                                    gTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGc                                 <  2:3342455/100‑1 (MQ=255)
                                                                     tACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCt                                >  2:1733057/1‑100 (MQ=255)
                                                                      aCTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTc                               <  1:1133508/100‑1 (MQ=255)
                                                                      aCTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTc                               >  1:241217/1‑100 (MQ=255)
                                                                      aCTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTc                               >  2:2352400/1‑100 (MQ=255)
                                                                      aCTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTc                               <  2:1234309/100‑1 (MQ=255)
                                                                       cTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTca                              >  1:637860/1‑100 (MQ=255)
                                                                         tGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTcaca                            >  2:1614948/1‑100 (MQ=255)
                                                                            cATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCt                         >  2:966928/1‑100 (MQ=255)
                                                                               cGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCtttt                      <  1:4373230/100‑1 (MQ=255)
                                                                                 cacaGGCCAACGTCAGCACATATCTTGCTGCTTTTGACTCACAAGTGGGTTTAACAGTCTTTGTCCGAACCACCCGGCAGCTGATGCTCACGCCTTTTgg                    <  2:4126587/100‑1 (MQ=255)
                                                                                 cacaGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACTGAACCACCCGGCAGCTGATGCTCACACCTTTTgg                    <  2:3925579/100‑1 (MQ=255)
                                                                                  acaGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGc                   <  2:2459944/100‑1 (MQ=255)
                                                                                  acaGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGc                   <  1:391216/100‑1 (MQ=255)
                                                                                    aGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCAc                 <  2:3637207/100‑1 (MQ=255)
                                                                                     attcaaaatCAGCAAGATATTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGCCCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCAc                 >  2:40294/9‑99 (MQ=37)
                                                                                      ttcaaaatCAGCAAGATATTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACcg               >  1:773050/8‑100 (MQ=38)
                                                                                       tcaaaatCAGCAAGATATTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACcgc              >  1:2468643/7‑100 (MQ=255)
                                                                                       ccAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACcgc              >  1:3143515/1‑100 (MQ=255)
                                                                                           aatCAGCAAGATATTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGACACCGCGCTg          >  1:4123628/3‑100 (MQ=255)
                                                                                             tCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGtt        <  2:1664544/100‑1 (MQ=255)
                                                                                              cAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTa       >  2:3579465/1‑100 (MQ=255)
                                                                                               aGCAAATATCTTGCTGATTTTTAATCTAAAGTGTGTTTAAAAGTCTTTGATCGAACCAGCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTAc      <  1:3021115/100‑1 (MQ=255)
                                                                                               aGCAAATATCTTGCTGATTTTAAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTAc      <  1:2769675/100‑1 (MQ=255)
                                                                                                gCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTAcc     >  2:328486/1‑100 (MQ=255)
                                                                                                gCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTAcc     >  1:2678708/1‑100 (MQ=255)
                                                                                                gCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTAcc     >  1:3193177/1‑100 (MQ=255)
                                                                                                gCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTAcc     >  2:295573/1‑100 (MQ=255)
                                                                                                 cAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCg    <  2:2052534/100‑1 (MQ=255)
                                                                                                  aaaTATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGt   >  1:3850399/1‑100 (MQ=255)
                                                                                                  aaaTATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGt   >  1:916876/1‑100 (MQ=255)
                                                                                                   aaTATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTa  <  1:4031184/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
TTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTA  >  NC_000913/3752008‑3752206

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