Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,471,300 T→A 18.1% Q73L (CAG→CTG)  yfdK ← CPS‑53 (KpLE1) prophage; conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,471,3000TA18.1% 214.4 / 37.2 90Q73L (CAG→CTG) yfdKCPS‑53 (KpLE1) prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base T (34/39);  new base A (11/5);  total (46/44)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

AGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGCGGTAATCCGGAAAACTCAA  >  NC_000913/2471201‑2471399
                                                                                                   |                                                                                                   
aGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCt                                                                                                     >  1:1755300/1‑100 (MQ=255)
aGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCt                                                                                                     >  2:809405/1‑100 (MQ=255)
   cgacCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGtt                                                                                                  <  2:1818222/100‑1 (MQ=255)
    gacCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTc                                                                                                 >  2:833087/1‑100 (MQ=255)
     acCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTcc                                                                                                <  1:1105481/100‑1 (MQ=255)
     acCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTcc                                                                                                >  1:1859189/1‑100 (MQ=255)
     acCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTcc                                                                                                >  1:2032806/1‑100 (MQ=255)
     acCAATAGCCGCTTTAAAAACCCAAAAATTAATGTACATATAATCGTAGGCCTGATTAATCAATTACTGCTTTTTAAGTTCGGCTGCAGCAAACTGTTcc                                                                                                >  2:955794/1‑100 (MQ=255)
      ccAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTcat                                                                                               >  1:185547/1‑98 (MQ=255)
        aaTAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCa                                                                                             <  1:1188942/100‑1 (MQ=255)
        aaTAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCa                                                                                             <  1:943801/100‑1 (MQ=255)
        aaTAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCa                                                                                             >  1:68703/1‑100 (MQ=255)
           aGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCAtgt                                                                                          <  1:648221/100‑1 (MQ=255)
            gCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCAtgt                                                                                          <  1:648231/99‑1 (MQ=255)
               gCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCAGTTCCTCATGTGTtg                                                                                      >  1:829014/1‑100 (MQ=255)
                 tctACCAGCCCATTGTTTACTGTTCAGATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTCCGGCTGCAGCAATCTGTTCCTCATGTGTtggt                                                                                    <  2:827150/98‑1 (MQ=255)
                  ttACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTtggtg                                                                                   <  2:1364946/100‑1 (MQ=255)
                   tACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTtggtgg                                                                                  >  2:511800/1‑100 (MQ=255)
                   tACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTtggtgg                                                                                  <  2:155442/100‑1 (MQ=255)
                   tACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCAGTTCCTCATGTGTtggtgg                                                                                  <  2:464717/100‑1 (MQ=255)
                    aCCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTtggtggt                                                                                 <  1:1298468/100‑1 (MQ=255)
                     ccAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTtggtggtg                                                                                <  1:495894/100‑1 (MQ=255)
                      cAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTtggtggtgg                                                                               <  2:1628948/100‑1 (MQ=255)
                        gCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGaa                                                                             >  2:2189813/1‑100 (MQ=255)
                         cccATTGTTTACTGTTCATATAATCGTTGGCCTGATTCATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAAt                                                                            >  2:1276318/1‑100 (MQ=255)
                          ccATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCAGTTCCTCATGTGTTGGTGGTGGAAtt                                                                           >  2:1483484/1‑100 (MQ=255)
                           cATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAAttt                                                                          <  2:401140/100‑1 (MQ=255)
                            aTTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCAGTTCCTCATGTGTTGGTGGTGGAATTTc                                                                         >  2:1882743/1‑100 (MQ=255)
                             ttGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCa                                                                        <  1:353996/100‑1 (MQ=255)
                             ttGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCa                                                                        <  2:499221/100‑1 (MQ=255)
                             ttGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCAGTTCCTCATGTGTTGGTGGTGGAATTTCa                                                                        >  1:54429/1‑100 (MQ=255)
                              tctttACTGTTCACATAATCGTTGGCCTGATTAATCAATTGCCGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAg                                                                       >  1:1622543/3‑100 (MQ=255)
                              tGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCa                                                                        >  2:574940/1‑99 (MQ=255)
                              tGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAg                                                                       <  1:777834/100‑1 (MQ=255)
                               gTTTACTGTTCATATTATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCAGTTCCTCATGTGTTGGTGGTGGAATTTCAGa                                                                      <  1:2370200/100‑1 (MQ=255)
                                   aCTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCAt                                                                  >  1:1307833/1‑100 (MQ=255)
                                   aCTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCAt                                                                  >  1:1747687/1‑100 (MQ=255)
                                     tGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGc                                                                <  1:2044494/100‑1 (MQ=255)
                                     tGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGc                                                                >  1:414768/1‑100 (MQ=255)
                                     tGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGc                                                                <  1:252911/100‑1 (MQ=255)
                                      nttCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCAGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCa                                                               >  1:750596/2‑100 (MQ=255)
                                        tCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTTGTGGAATTTCAGACCATGCAgg                                                             <  2:341790/100‑1 (MQ=255)
                                        tCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAgg                                                             <  2:1629596/100‑1 (MQ=255)
                                        tCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAgg                                                             <  1:1200217/100‑1 (MQ=255)
                                         cATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCAGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGa                                                            >  2:1504038/1‑100 (MQ=255)
                                           tataATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGaaa                                                          >  1:1273511/1‑100 (MQ=255)
                                            ataATCGTTGGCCTGATTAATCAATTGCTGCTTTTTTTTTTTTTTTTTAGCAATCTGTTCCTAATGTGTTGTTAATGGAATTTCAGTCCATGCAGGata                                                          >  2:327360/1‑97 (MQ=255)
                                             tcaTCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTACTCATGTGTTGGTGTTTGAATTTCAGACCATGCACGAAAAc                                                        >  2:1866518/3‑100 (MQ=255)
                                             taATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGCCCATGCAGGAAAAc                                                        >  1:2200200/1‑100 (MQ=255)
                                              aaTCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCAGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAAcc                                                       >  1:1536304/1‑100 (MQ=255)
                                                    tGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTtc                                                 <  2:544061/100‑1 (MQ=255)
                                                     ggCCTGATTAATCAATTGCTGCTTTTTCAGTTCTGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGACTTTCAGACCATGCAGGAAAACCATTTtct                                                >  2:1434784/1‑100 (MQ=255)
                                                      gCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTtctc                                               >  1:105723/1‑100 (MQ=255)
                                                      gCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCAGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTtctc                                               >  1:531658/1‑100 (MQ=255)
                                                       ccTGATTAATCAATTGCTGCTTTTTTTTTTTTTTTTCAGCAATAAATTTCTCATATTTTTGTGGTGGAATTTCAGACCATGCAGGAAAACTATTTTCTcc                                              >  2:104900/1‑100 (MQ=255)
                                                       ccTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTcc                                              <  1:833096/100‑1 (MQ=255)
                                                         tGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAg                                            >  2:2191196/1‑100 (MQ=255)
                                                          gATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGc                                           >  1:235485/1‑100 (MQ=255)
                                                           aTTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTACCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCg                                          <  2:1536035/100‑1 (MQ=255)
                                                               aTCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATAc                                      <  2:909420/100‑1 (MQ=255)
                                                               aTCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATAc                                      <  1:746923/100‑1 (MQ=255)
                                                                  aaTTGCTGCTTTTTCAGTTCGGCTGCAGCAATCAGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGa                                   <  2:54428/100‑1 (MQ=255)
                                                                    ttGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGAtt                                 >  2:1262744/1‑100 (MQ=255)
                                                                    ttGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGAtt                                 <  1:619332/100‑1 (MQ=255)
                                                                    ttGCTGCTTGTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGTAATTTCAGACCATGCAGTAAAATCATTTTCTCCAGCGATACGGAt                                  <  2:1919767/99‑1 (MQ=255)
                                                                       ctgctTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGAttttt                              >  1:1388221/1‑100 (MQ=255)
                                                                       ctgctTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGAttttt                              <  1:790955/100‑1 (MQ=255)
                                                                       ctgctTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGAttttt                              >  2:725467/1‑100 (MQ=255)
                                                                       ctgctTTTTCAGTTCGGCTGCAGCAATCAGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGAAttttt                              >  1:1080142/1‑100 (MQ=255)
                                                                        tgctTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTc                             <  2:1690107/100‑1 (MQ=255)
                                                                         gctTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTcc                            >  2:1084637/1‑100 (MQ=255)
                                                                         gctTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTcc                            >  2:184475/1‑100 (MQ=255)
                                                                          ctTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCTGCGATACGGATTTTTCCt                           <  2:947861/100‑1 (MQ=255)
                                                                           tttttCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCtt                          <  2:1860553/100‑1 (MQ=255)
                                                                            ttttCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCttt                         >  1:1122714/1‑100 (MQ=255)
                                                                             tttCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTc                        <  1:353873/100‑1 (MQ=255)
                                                                               tCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTcgg                      <  1:1375543/100‑1 (MQ=255)
                                                                                cAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTcggc                     <  2:105721/100‑1 (MQ=255)
                                                                                cAGTTCGGCTGCAGCAATCCGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTcgg                      >  1:603927/1‑99 (MQ=255)
                                                                                cAGTTCGGCTGCAGCAATCAGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTcggc                     <  1:1504055/100‑1 (MQ=255)
                                                                                  gTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTcggcgg                   >  1:1050161/1‑100 (MQ=255)
                                                                                  gTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTcggcgg                   <  2:414765/100‑1 (MQ=255)
                                                                                  gTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTcggcgg                   <  2:2079172/100‑1 (MQ=255)
                                                                                      ggCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGCGGTAAt               >  1:1914319/1‑100 (MQ=255)
                                                                                      ggCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGCGGTAAt               >  2:1761330/1‑100 (MQ=255)
                                                                                       gCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGCGGTAAt               >  2:61383/1‑99 (MQ=255)
                                                                                          gcagcaATCAGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGCGGTAATCCgg           >  1:804928/1‑100 (MQ=255)
                                                                                             gcaATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGGGGTAATCCGGaca        <  1:1626681/100‑3 (MQ=255)
                                                                                             gcaATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGCGGTAATCCGGaaa        >  1:1846171/1‑100 (MQ=255)
                                                                                             gcaATCAGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGCGGTAATCCGGaaa        <  2:1539835/100‑1 (MQ=255)
                                                                                                aTCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGCGGTAATCCGGAAAACt     <  2:2346144/100‑1 (MQ=255)
                                                                                                 tCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGCAGTAATCCGGAAAAc      >  2:2086657/1‑98 (MQ=255)
                                                                                                   tGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGCGGTAATCCGGAAAACTCaa  >  1:234734/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGCGGTAATCCGGAAAACTCAA  >  NC_000913/2471201‑2471399

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: