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→TGC) yafOmRNA 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.

TGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGT  >  NC_000913/252495‑252693
                                                                                                    |                                                                                                   
tGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGACAAAGCATgg                                                                                                     >  2:3740051/1‑100 (MQ=255)
  agaATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGtt                                                                                                   >  1:3685738/1‑100 (MQ=255)
   gaATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTg                                                                                                  >  1:3023398/1‑100 (MQ=255)
   gaATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTg                                                                                                  >  2:2154209/1‑100 (MQ=255)
    aataatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGGCGTTTGATGAGCAAGCATGGTTg                                                                                                  <  2:794573/100‑1 (MQ=255)
    aataatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGc                                                                                                 <  2:3825154/100‑1 (MQ=255)
    aataatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGc                                                                                                 <  1:2347964/100‑1 (MQ=255)
      taatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTc                                                                                               <  1:1653982/100‑1 (MQ=255)
         tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAGGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                            <  2:2126473/100‑2 (MQ=255)
         tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                            <  1:2603216/100‑2 (MQ=255)
         tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                            <  2:2359787/100‑2 (MQ=255)
         tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                            <  2:3124003/100‑2 (MQ=255)
         tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                            <  1:1422946/100‑2 (MQ=255)
         tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAACGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCAt                                                                                             <  2:3068163/99‑1 (MQ=255)
          ccATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCAtt                                                                                            >  2:896018/1‑99 (MQ=255)
          ccATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTg                                                                                           <  1:3032127/100‑1 (MQ=255)
           cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTTGTATATTGTCA‑GGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaa                                                                                          <  2:164045/100‑4 (MQ=255)
           cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGc                                                                                          <  2:1540961/100‑1 (MQ=255)
           cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGc                                                                                          <  1:3410839/100‑1 (MQ=255)
           cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGc                                                                                          >  1:3361643/1‑100 (MQ=255)
           cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGc                                                                                          <  1:1060339/100‑1 (MQ=255)
           cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaa                                                                                          <  1:873663/100‑4 (MQ=255)
           cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaa                                                                                          <  1:1457464/100‑4 (MQ=255)
            aTTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGcc                                                                                         <  2:3716142/100‑1 (MQ=255)
            aTTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGcc                                                                                         >  1:3103015/1‑100 (MQ=255)
            aTTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGcc                                                                                         <  1:318922/100‑1 (MQ=255)
            aTTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaaa                                                                                         <  1:299291/100‑5 (MQ=39)
              ttCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCCTTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaaac                                                                                        <  2:1839538/99‑6 (MQ=38)
              ttCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCAt                                                                                       <  2:710678/100‑1 (MQ=255)
               tCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCAtt                                                                                      >  1:1620589/1‑100 (MQ=255)
                ccGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCAt                                                                                       <  1:2127800/98‑1 (MQ=255)
                  gCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑TGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTg                                                                                   >  1:1161425/1‑100 (MQ=255)
                     acaGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGaaa                                                                                <  2:3055547/100‑1 (MQ=255)
                      caGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAAc                                                                               <  1:537199/100‑1 (MQ=255)
                       aGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAAcc                                                                              <  2:2991503/100‑1 (MQ=255)
                       aGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAAcc                                                                              <  1:1725571/100‑1 (MQ=255)
                          tGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGa                                                                           <  1:4083027/100‑1 (MQ=255)
                           gcgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGaa                                                                          >  1:524987/1‑100 (MQ=255)
                            cgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAc                                                                         >  1:184978/1‑100 (MQ=255)
                            cgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAc                                                                         <  1:2125797/100‑1 (MQ=255)
                             gcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAcc                                                                        <  1:462585/100‑1 (MQ=255)
                              ccAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAc                                                                         <  2:121043/98‑1 (MQ=255)
                                aaTTCAGCAGAGCGGATGACGAAGCGCATTTGGTATATTGCCA‑GGGGCCGCTTGATGAGGAAGCTTGGTTGCTCATTGCCATGCGGaatcttgt                                                                           <  2:1486040/94‑7 (MQ=255)
                                aaTTCAGCAGAGCGAATGACGAAGCGCATTTGGTATATTGTCA‑TGGGGCTTTTGATGAGCAAGCATGGTTGCTCATTACCATTCTGAAACCTGAACCTCa                                                                     <  2:1326298/100‑1 (MQ=255)
                                aaTTCAGCAGAACGAATGACTAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCa                                                                     >  2:2517111/1‑100 (MQ=255)
                                aaTTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCa                                                                     <  2:3347016/100‑1 (MQ=255)
                                    cagcagAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATaaa                                                                 >  1:3546037/1‑100 (MQ=255)
                                     ngcagAACGAATGACGAAGCGCATCTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAAc                                                                >  2:2853188/2‑100 (MQ=255)
                                      gcagAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACt                                                               <  1:1392796/100‑1 (MQ=255)
                                        agAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTgg                                                             <  2:2375346/100‑1 (MQ=255)
                                        agAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTgg                                                             <  1:3486153/100‑1 (MQ=255)
                                         gAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGc                                                            >  1:2730073/1‑100 (MQ=255)
                                          aaagAATCACGAAGCGCATTTGGTATCCTCTCA‑GCGGGGGTATGATCAGCAACCATGCTAGCTCATCGCCATTCTGAAACCTGACCCACATGAACTGGCt                                                           <  2:2249666/97‑1 (MQ=255)
                                             gAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCga                                                        >  1:859507/1‑100 (MQ=255)
                                             gAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCga                                                        <  1:585130/100‑1 (MQ=255)
                                                 gACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATa                                                    <  1:2538951/100‑1 (MQ=255)
                                                  aCGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGGGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaa                                                   >  2:1081011/1‑100 (MQ=255)
                                                  aCGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaa                                                   >  1:324106/1‑100 (MQ=255)
                                                   cGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaac                                                  >  1:2815806/1‑100 (MQ=255)
                                                    gAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaaca                                                 >  1:2411308/1‑100 (MQ=255)
                                                          cATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCaaa                                           >  2:1250165/1‑100 (MQ=255)
                                                             ttGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGc                                        <  1:3641020/100‑1 (MQ=255)
                                                                gTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATa                                     >  2:2769091/1‑100 (MQ=255)
                                                                  atatTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGATACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATaaa                                   <  2:2290041/100‑1 (MQ=255)
                                                                       gTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGGACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTgg                              <  1:2379688/100‑1 (MQ=255)
                                                                        tCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTggg                             <  2:3518176/100‑1 (MQ=255)
                                                                        tCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTggg                             <  1:1809877/100‑1 (MQ=255)
                                                                         cA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGa                            >  2:3413889/1‑100 (MQ=255)
                                                                         cA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGa                            >  1:289588/1‑100 (MQ=255)
                                                                          a‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGaa                           <  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)
                                                                                                    |                                                                                                   
TGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCA‑GGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGT  >  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: ‑