Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 252,589 G→C 20.2% A97P (GCA→CCA)  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,5890GC20.2% 351.9 / 92.7 154A97P (GCA→CCA) yafOmRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (73/50);  new base C (18/13);  total (91/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-01

Only 100 of 174 total aligned reads displayed.

GCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCA  >  NC_000913/252490‑252688
                                                                                                   |                                                                                                   
gCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCTAATTCAGCAGAACGAACGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGagcaag                                                                                                     <  1:1825492/100‑1 (MQ=255)
gCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGagcaag                                                                                                     <  1:3144811/100‑1 (MQ=255)
 cAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTGGGTATATTGTCAGGGGGCGTTTGATGagcaag                                                                                                     <  2:2916933/99‑1 (MQ=255)
     tGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATgg                                                                                                >  1:2782690/1‑100 (MQ=255)
     tGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATgg                                                                                                >  1:1305987/1‑100 (MQ=255)
     tGAGAATAATCAATTTCCGCCACAGTTGCGCCAATTCAGTAGAAAGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATgg                                                                                                >  2:816838/1‑100 (MQ=255)
        gaATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTg                                                                                             <  2:3971585/100‑1 (MQ=255)
         aataatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGc                                                                                            <  2:1676467/100‑1 (MQ=255)
           taatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTc                                                                                          >  2:2095901/1‑100 (MQ=255)
            aatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCa                                                                                         >  1:4166532/1‑100 (MQ=255)
            aatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCa                                                                                         >  1:3084579/1‑100 (MQ=255)
             atCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCAt                                                                                        >  2:1907786/1‑100 (MQ=255)
             atCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCAt                                                                                        >  1:20971/1‑100 (MQ=255)
             atCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAt                                                                                        <  1:2380048/100‑1 (MQ=255)
              tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                       <  1:859703/100‑2 (MQ=255)
              tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                       <  1:2576608/100‑2 (MQ=255)
              tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                       <  2:3231801/100‑2 (MQ=255)
              tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                       <  1:1479897/100‑2 (MQ=255)
               ccATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtca                                                                                      <  1:1691361/100‑3 (MQ=255)
                cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGc                                                                                     <  1:1976057/100‑1 (MQ=255)
                 aTTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaaa                                                                                    <  1:111650/100‑5 (MQ=39)
                 aTTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaaa                                                                                    <  1:1389873/100‑5 (MQ=39)
                     cgtccACACTTGCGCCAATTCATCCGATCGAATAACTGAGTGCATTTTGTATATTGTCAGGTGGCGTTTGCTGATCAAGCATGGTTGCTCATTGCCATTc                                                                                <  2:2594751/97‑1 (MQ=255)
                      cGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGTAAGCATGGTTGCTCATTGCCa                                                                                   >  2:1906557/1‑96 (MQ=255)
                       gCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTg                                                                              <  1:2053754/100‑1 (MQ=255)
                           caGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAAc                                                                          <  1:963869/100‑1 (MQ=255)
                           caGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAAc                                                                          <  1:3253689/100‑1 (MQ=255)
                            aGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAAcc                                                                         <  1:1541196/100‑1 (MQ=255)
                              ttGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTg                                                                       <  1:3489674/100‑1 (MQ=255)
                                gcgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGaa                                                                     >  1:3866440/1‑100 (MQ=255)
                                 cgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAc                                                                    <  1:318671/100‑1 (MQ=255)
                                   ccAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCt                                                                  >  2:4001793/1‑100 (MQ=255)
                                   ccAATTCAGCAGAACGAATGACGAAGAGCATTTGGTATATTGTCAGGGAGCGTTTGATGAGCAAGCATGGTTGCACATTGCCATTCTGAAACCTGAAc                                                                    >  1:3417142/1‑98 (MQ=255)
                                      aTTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCAt                                                               >  1:3830699/1‑100 (MQ=255)
                                       ttCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATa                                                              <  2:2162130/100‑1 (MQ=255)
                                          aggagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAAGAAGTATGGCTGATCATTGCCATACTGAAACGTGaatat                                                                  >  2:3471414/4‑90 (MQ=255)
                                          agcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAAc                                                           >  1:1386258/1‑100 (MQ=255)
                                           gcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACt                                                          >  2:2958093/1‑100 (MQ=255)
                                           gcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACt                                                          <  2:3295394/100‑1 (MQ=255)
                                               aaCGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGTGTTTGATGAGCAAGCATGGTTGCTCATTGCCCTTCTGAAACCTGAACCTCATAAACTGGCt                                                      >  2:2267437/1‑100 (MQ=255)
                                                 cGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCg                                                    <  2:3637445/100‑1 (MQ=255)
                                                  gAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTc                                                     <  1:190475/98‑1 (MQ=255)
                                                  gAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCg                                                    >  2:1313398/1‑99 (MQ=255)
                                                   aaTGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgag                                                  >  2:3839723/1‑100 (MQ=255)
                                                      gACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATa                                               <  2:714064/100‑1 (MQ=255)
                                                       aCGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaa                                              >  2:1859236/1‑100 (MQ=255)
                                                        cGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaac                                             >  2:2700345/1‑100 (MQ=255)
                                                        cGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaac                                             >  2:1650573/1‑100 (MQ=255)
                                                         gAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaaca                                            <  2:3139319/100‑1 (MQ=255)
                                                         gAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaaca                                            <  2:625544/100‑1 (MQ=255)
                                                           aGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaacaac                                          >  2:739148/1‑100 (MQ=255)
                                                            gcgcATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAAcc                                         <  1:3470087/100‑1 (MQ=255)
                                                              gcATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCaa                                       >  2:2475909/1‑100 (MQ=255)
                                                               cATTTGGTATATTGTCAGGGGGAGTTTGATGAGCAAGCATGGTTGCTCCTTGCCATTCTGCAACCTGAACCTCATAAACTGGCTCGAGATAACAACCaaa                                      >  2:196415/1‑100 (MQ=255)
                                                                 tttGGTAGATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCCCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGTTAACAACCAAAt                                     >  2:1393439/1‑99 (MQ=255)
                                                                  ttGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGc                                   >  2:3552667/1‑100 (MQ=255)
                                                                    ggTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCAt                                 >  2:3587147/1‑100 (MQ=255)
                                                                     gTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCaga                                >  2:2549754/1‑98 (MQ=255)
                                                                     gTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATa                                >  2:2688290/1‑100 (MQ=255)
                                                                      tataTTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATaa                               >  1:1455564/1‑100 (MQ=255)
                                                                      tataTTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATaa                               >  2:2243499/1‑100 (MQ=255)
                                                                       atatTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATaaa                              >  2:1320314/1‑100 (MQ=255)
                                                                         atTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAAt                            <  2:3290879/100‑1 (MQ=255)
                                                                          ttGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAAtt                           >  1:553353/1‑100 (MQ=255)
                                                                           tGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATGGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTg                          <  1:816838/100‑1 (MQ=255)
                                                                            gTCAGGGCGCGTTTGATGAGCAAGCATGGGTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTgg                         <  2:3858101/100‑1 (MQ=255)
                                                                             tCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTggg                        <  2:2993259/100‑1 (MQ=255)
                                                                             tCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTggg                        <  1:3902466/100‑1 (MQ=255)
                                                                                gggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGaaa                     >  2:1115031/1‑100 (MQ=255)
                                                                                   ggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATg                  >  1:3812884/1‑100 (MQ=255)
                                                                                   ggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATg                  >  1:2809174/1‑100 (MQ=255)
                                                                                   ggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATg                  >  1:2928773/1‑100 (MQ=255)
                                                                                   ggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACATGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATg                  <  2:2609677/100‑1 (MQ=255)
                                                                                   ggCGTTTGATGAGCAAGCATGGTTGCTCATGGCCATTCTGAAACTTGAACCCCATAAACTGGCTCGAGATAACATCCAAATGCATAAAATTGGGAAAATg                  <  2:1925290/100‑1 (MQ=255)
                                                                                    gCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATgg                 <  1:322254/100‑1 (MQ=255)
                                                                                     cGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                <  2:941871/100‑1 (MQ=255)
                                                                                     cGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                <  1:1435259/100‑1 (MQ=255)
                                                                                     cGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                >  2:680297/1‑100 (MQ=255)
                                                                                       tttGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGAATAAAATTGGGAAAATGGCa               >  2:3757012/1‑99 (MQ=255)
                                                                                       ggcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAg              >  1:1873092/5‑100 (MQ=39)
                                                                                        gcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGa             >  2:252585/4‑100 (MQ=255)
                                                                                        gcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGAAAACAACCAAATGCATAAAATTGGGAAAATGGCAGa             >  1:376783/4‑100 (MQ=255)
                                                                                         tgatgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGaa            >  1:2681424/1‑100 (MQ=255)
                                                                                         caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGaa            >  1:2490793/3‑100 (MQ=255)
                                                                                         caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGaa            >  2:2016900/3‑100 (MQ=255)
                                                                                         caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGaa            >  1:1741657/3‑100 (MQ=255)
                                                                                         caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGaa            >  1:4111512/3‑100 (MQ=255)
                                                                                          aatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAACTGCATAAAATTGGGAAAATGGCAGAAg           >  2:724086/2‑100 (MQ=255)
                                                                                          aatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAg           >  1:2487222/2‑100 (MQ=255)
                                                                                          aatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAg           >  2:2771950/2‑100 (MQ=255)
                                                                                              agcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGtt       >  1:1017240/1‑100 (MQ=255)
                                                                                               gcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGttt      <  1:3669099/100‑1 (MQ=255)
                                                                                                caagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTc     >  2:4125266/1‑100 (MQ=255)
                                                                                                 aagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCg    <  2:2029218/100‑1 (MQ=255)
                                                                                                 aagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCg    <  1:764180/100‑1 (MQ=255)
                                                                                                 aagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCg    >  1:3090616/1‑100 (MQ=255)
                                                                                                  agcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGc   >  2:524452/1‑100 (MQ=255)
                                                                                                  agcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGc   >  1:1619744/1‑100 (MQ=255)
                                                                                                  agcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGc   >  1:3565035/1‑100 (MQ=255)
                                                                                                   gcaTGGTTGCCAATTGCCATTCTGAAACCTGAACCTCATAAGCTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGATAGAAGCGTTTCGCa  <  2:1771509/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCA  >  NC_000913/252490‑252688

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: