Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 252,594 G→C 19.4% W98C (TGG→TGC yafO → mRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913252,5940GC19.4% 424.8 / 97.0 175W98C (TGG→TGCyafOmRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (59/82);  new base C (17/17);  total (76/99)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.88e-01

Only 100 of 178 total aligned reads displayed.

TGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGT  >  NC_000913/252495‑252693
                                                                                                    |                                                                                                   
tGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGACAAAGCATgg                                                                                                     >  2:3740051/1‑100 (MQ=255)
  agaATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGtt                                                                                                   >  1:3685738/1‑100 (MQ=255)
   gaATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTg                                                                                                  >  1:3023398/1‑100 (MQ=255)
   gaATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTg                                                                                                  >  2:2154209/1‑100 (MQ=255)
    aataatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGGCGTTTGATGAGCAAGCATGGTTg                                                                                                  <  2:794573/100‑1 (MQ=255)
    aataatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGc                                                                                                 <  2:3825154/100‑1 (MQ=255)
    aataatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGc                                                                                                 <  1:2347964/100‑1 (MQ=255)
      taatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTc                                                                                               <  1:1653982/100‑1 (MQ=255)
         tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAGGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                            <  2:2126473/100‑2 (MQ=255)
         tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                            <  1:2603216/100‑2 (MQ=255)
         tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                            <  2:2359787/100‑2 (MQ=255)
         tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                            <  2:3124003/100‑2 (MQ=255)
         tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                            <  1:1422946/100‑2 (MQ=255)
         tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAACGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCAt                                                                                             <  2:3068163/99‑1 (MQ=255)
          ccATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCAtt                                                                                            >  2:896018/1‑99 (MQ=255)
          ccATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTg                                                                                           <  1:3032127/100‑1 (MQ=255)
           cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTTGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaa                                                                                          <  2:164045/100‑4 (MQ=255)
           cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGc                                                                                          <  2:1540961/100‑1 (MQ=255)
           cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGc                                                                                          <  1:3410839/100‑1 (MQ=255)
           cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGc                                                                                          >  1:3361643/1‑100 (MQ=255)
           cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGc                                                                                          <  1:1060339/100‑1 (MQ=255)
           cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaa                                                                                          <  1:873663/100‑4 (MQ=255)
           cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaa                                                                                          <  1:1457464/100‑4 (MQ=255)
            aTTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGcc                                                                                         <  2:3716142/100‑1 (MQ=255)
            aTTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGcc                                                                                         >  1:3103015/1‑100 (MQ=255)
            aTTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGcc                                                                                         <  1:318922/100‑1 (MQ=255)
            aTTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaaa                                                                                         <  1:299291/100‑5 (MQ=39)
              ttCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCCTTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaaac                                                                                        <  2:1839538/99‑6 (MQ=38)
              ttCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCAt                                                                                       <  2:710678/100‑1 (MQ=255)
               tCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCAtt                                                                                      >  1:1620589/1‑100 (MQ=255)
                ccGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCAt                                                                                       <  1:2127800/98‑1 (MQ=255)
                  gCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCATGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTg                                                                                   >  1:1161425/1‑100 (MQ=255)
                     acaGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGaaa                                                                                <  2:3055547/100‑1 (MQ=255)
                      caGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAAc                                                                               <  1:537199/100‑1 (MQ=255)
                       aGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAAcc                                                                              <  2:2991503/100‑1 (MQ=255)
                       aGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAAcc                                                                              <  1:1725571/100‑1 (MQ=255)
                          tGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGa                                                                           <  1:4083027/100‑1 (MQ=255)
                           gcgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGaa                                                                          >  1:524987/1‑100 (MQ=255)
                            cgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAc                                                                         >  1:184978/1‑100 (MQ=255)
                            cgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAc                                                                         <  1:2125797/100‑1 (MQ=255)
                             gcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAcc                                                                        <  1:462585/100‑1 (MQ=255)
                              ccAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAc                                                                         <  2:121043/98‑1 (MQ=255)
                                aaTTCAGCAGAGCGGATGACGAAGCGCATTTGGTATATTGCCAGGGGCCGCTTGATGAGGAAGCTTGGTTGCTCATTGCCATGCGGaatcttgt                                                                           <  2:1486040/94‑7 (MQ=255)
                                aaTTCAGCAGAGCGAATGACGAAGCGCATTTGGTATATTGTCATGGGGCTTTTGATGAGCAAGCATGGTTGCTCATTACCATTCTGAAACCTGAACCTCa                                                                     <  2:1326298/100‑1 (MQ=255)
                                aaTTCAGCAGAACGAATGACTAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCa                                                                     >  2:2517111/1‑100 (MQ=255)
                                aaTTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCa                                                                     <  2:3347016/100‑1 (MQ=255)
                                    cagcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATaaa                                                                 >  1:3546037/1‑100 (MQ=255)
                                     ngcagAACGAATGACGAAGCGCATCTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAAc                                                                >  2:2853188/2‑100 (MQ=255)
                                      gcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACt                                                               <  1:1392796/100‑1 (MQ=255)
                                        agAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTgg                                                             <  2:2375346/100‑1 (MQ=255)
                                        agAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTgg                                                             <  1:3486153/100‑1 (MQ=255)
                                         gAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGc                                                            >  1:2730073/1‑100 (MQ=255)
                                          aaagAATCACGAAGCGCATTTGGTATCCTCTCAGCGGGGGTATGATCAGCAACCATGCTAGCTCATCGCCATTCTGAAACCTGACCCACATGAACTGGCt                                                           <  2:2249666/97‑1 (MQ=255)
                                             gAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCga                                                        >  1:859507/1‑100 (MQ=255)
                                             gAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCga                                                        <  1:585130/100‑1 (MQ=255)
                                                 gACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATa                                                    <  1:2538951/100‑1 (MQ=255)
                                                  aCGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGGGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaa                                                   >  2:1081011/1‑100 (MQ=255)
                                                  aCGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaa                                                   >  1:324106/1‑100 (MQ=255)
                                                   cGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaac                                                  >  1:2815806/1‑100 (MQ=255)
                                                    gAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaaca                                                 >  1:2411308/1‑100 (MQ=255)
                                                          cATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCaaa                                           >  2:1250165/1‑100 (MQ=255)
                                                             ttGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGc                                        <  1:3641020/100‑1 (MQ=255)
                                                                gTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATa                                     >  2:2769091/1‑100 (MQ=255)
                                                                  atatTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGATACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATaaa                                   <  2:2290041/100‑1 (MQ=255)
                                                                       gTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGGACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTgg                              <  1:2379688/100‑1 (MQ=255)
                                                                        tCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTggg                             <  2:3518176/100‑1 (MQ=255)
                                                                        tCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTggg                             <  1:1809877/100‑1 (MQ=255)
                                                                         cAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGa                            >  2:3413889/1‑100 (MQ=255)
                                                                         cAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGa                            >  1:289588/1‑100 (MQ=255)
                                                                          aGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGaa                           <  1:3574200/100‑1 (MQ=255)
                                                                             ggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGaaaa                         >  2:859855/1‑100 (MQ=255)
                                                                             ggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGaaaa                         >  2:1502395/1‑100 (MQ=255)
                                                                                 cGTTTGATGCGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGCCTCGAAATAACAACCAAATGCATAAAATTGGGAAAATGGc                     >  2:3979174/1‑100 (MQ=255)
                                                                                 cGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                     <  1:2987775/100‑1 (MQ=255)
                                                                                 cGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                     >  2:1418206/1‑100 (MQ=255)
                                                                                 cGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                     <  1:2517125/100‑1 (MQ=255)
                                                                                 cGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                     >  2:3166265/1‑100 (MQ=255)
                                                                                  gTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCa                    <  2:2815790/100‑1 (MQ=255)
                                                                                    gcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGa                  >  2:593808/4‑100 (MQ=255)
                                                                                     tgatgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGaa                 <  2:1985651/100‑1 (MQ=255)
                                                                                     caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAg                   >  2:3284086/3‑98 (MQ=255)
                                                                                     caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGaa                 >  1:1040737/3‑100 (MQ=255)
                                                                                     caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGaa                 >  2:2627460/3‑100 (MQ=255)
                                                                                     caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGaa                 >  1:420874/3‑100 (MQ=255)
                                                                                      gatgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAg                <  1:1526938/100‑1 (MQ=255)
                                                                                      aatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAg                >  1:3991472/2‑100 (MQ=255)
                                                                                       atgaGCAAGCATTGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGc               <  2:3248248/100‑1 (MQ=255)
                                                                                       atgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGc               >  1:3551239/1‑100 (MQ=255)
                                                                                       atgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGc               >  2:759081/1‑100 (MQ=255)
                                                                                       atgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGc               >  2:759101/1‑100 (MQ=255)
                                                                                         gaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGCTAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCg              >  2:3885694/1‑99 (MQ=255)
                                                                                         gaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCACATGCATAAAATTGGGAAAATGGCAGAAGCGt             >  2:569538/1‑100 (MQ=255)
                                                                                          agcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGtt            >  2:2840661/1‑100 (MQ=255)
                                                                                          agcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGtt            >  1:3851006/1‑100 (MQ=255)
                                                                                          agcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGtt            <  1:3982670/100‑1 (MQ=255)
                                                                                             aagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGGTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCg         <  2:10218/100‑1 (MQ=255)
                                                                                              agcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGc        >  1:269553/1‑100 (MQ=255)
                                                                                               gcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCa       <  2:3973276/100‑1 (MQ=255)
                                                                                                 aTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATg     <  1:3153127/100‑1 (MQ=255)
                                                                                                    gTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGt  >  1:1702477/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
TGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGT  >  NC_000913/252495‑252693

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: