Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,777,363 T→C 28.7% F83F (TTT→TTC yfjZ → CP4‑57 prophage; antitoxin of the YpjF‑YfjZ toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,777,3630TC28.7% 148.2 / 76.7 108F83F (TTT→TTCyfjZCP4‑57 prophage; antitoxin of the YpjF‑YfjZ toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (39/38);  new base C (18/13);  total (57/51)
Fisher's exact test for biased strand distribution p-value = 5.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

Only 100 of 110 total aligned reads displayed.

TTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACAC  >  NC_000913/2777268‑2777460
                                                                                               |                                                                                                 
tttCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCt                                                                                               >  2:93615/1‑100 (MQ=255)
 ttCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTg                                                                                              <  1:128516/100‑1 (MQ=255)
 ttCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTg                                                                                              <  2:1127846/100‑1 (MQ=255)
  tCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGc                                                                                             >  2:658254/1‑100 (MQ=255)
  tCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGCTTTACCTGc                                                                                             >  2:2078969/1‑100 (MQ=255)
   ccACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCg                                                                                            >  2:255574/1‑100 (MQ=255)
   ccACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCg                                                                                            <  1:1542901/100‑1 (MQ=255)
    cacaTTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTAAATCACCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGa                                                                                           >  2:238022/1‑100 (MQ=255)
     acaTTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGaa                                                                                          <  2:1749203/100‑1 (MQ=255)
     acaTTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGaa                                                                                          >  1:1647742/1‑100 (MQ=255)
      caTTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAg                                                                                         >  2:312933/1‑100 (MQ=255)
      caTTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATACCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAg                                                                                         >  2:1998811/1‑100 (MQ=255)
        ttttATCAGCCAGATAGAGTCGATGCTGACCACTTGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGcc                                                                                       <  2:169540/100‑1 (MQ=255)
        ttttATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGcc                                                                                       <  1:2131529/100‑1 (MQ=255)
        ttttATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGcc                                                                                       <  2:1152923/100‑1 (MQ=255)
         tttATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCg                                                                                      <  2:323380/100‑1 (MQ=255)
          ttATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGa                                                                                     <  1:1540794/100‑1 (MQ=255)
            aTCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATa                                                                                   <  2:1349562/100‑1 (MQ=255)
              cAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATAct                                                                                 >  1:1867330/1‑100 (MQ=255)
              cAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATAct                                                                                 >  2:791420/1‑100 (MQ=255)
                 ccAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCtt                                                                              <  2:2365059/100‑1 (MQ=255)
                  cAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTg                                                                             >  1:2312662/1‑100 (MQ=255)
                   aGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTgg                                                                            <  2:739515/100‑1 (MQ=255)
                   aGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTgg                                                                            <  1:1273195/100‑1 (MQ=255)
                   aGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTgg                                                                            >  1:2345852/1‑100 (MQ=255)
                   aGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTgg                                                                            >  2:868915/1‑100 (MQ=255)
                   aGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTg                                                                             <  1:700163/99‑1 (MQ=255)
                     gtAGAGTCGATGCTGACCACTGGTGAACTGCATGCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGAGACTCTTGGTa                                                                          >  2:1578582/2‑100 (MQ=25)
                     aTAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTa                                                                          >  2:1995993/1‑100 (MQ=255)
                     aTAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTa                                                                          <  1:322230/100‑1 (MQ=255)
                       agagTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGt                                                                        <  2:83708/100‑1 (MQ=255)
                       agagTCGATGCTATCCACTGGTGAACTTAATCCCCGCCATGCCCAATGCGTCACCCTGTACTACAACGGTTTCACCTGCGAAGCCGATCCTCTTTGTa                                                                          >  2:759458/1‑98 (MQ=255)
                         agTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTg                                                                      >  2:88894/1‑100 (MQ=255)
                          gTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGc                                                                     >  2:399379/1‑100 (MQ=255)
                          gTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGc                                                                     <  2:210590/100‑1 (MQ=255)
                           tCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCTAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCg                                                                    <  2:1069378/100‑1 (MQ=255)
                             gATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGc                                                                  <  1:785999/100‑1 (MQ=255)
                             gATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGc                                                                  >  2:1573322/1‑100 (MQ=255)
                             gATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGc                                                                  >  2:1021419/1‑100 (MQ=255)
                              aTGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGc                                                                  <  1:789135/99‑1 (MQ=255)
                               tGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTa                                                                >  1:1003573/1‑100 (MQ=255)
                                 cTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACg                                                              <  2:1660616/100‑1 (MQ=255)
                                 cTGACCACGGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACg                                                              <  1:2350400/100‑1 (MQ=255)
                                  tGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGt                                                             >  2:681337/1‑100 (MQ=255)
                                   gACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGta                                                            <  2:590067/100‑1 (MQ=255)
                                     ccACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGtata                                                          >  1:2060721/1‑100 (MQ=255)
                                     ccACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGtata                                                          >  2:1948745/1‑100 (MQ=255)
                                       aCTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATTCTCTTGGTAGTTGCGGCTACGta                                                            <  2:268635/96‑1 (MQ=255)
                                            tGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCt                                                   >  1:204607/1‑100 (MQ=255)
                                               aCTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGtt                                                <  2:1787760/100‑1 (MQ=255)
                                                cTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGttt                                               <  2:1431814/100‑1 (MQ=255)
                                                 tGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGttt                                               <  2:2257486/99‑1 (MQ=255)
                                                 tGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTa                                              >  1:921818/1‑100 (MQ=255)
                                                  gAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAc                                             >  1:856500/1‑100 (MQ=255)
                                                   aaTCCCCGCCATGCGCAATGCGTCATCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAcc                                            <  2:1588251/100‑1 (MQ=255)
                                                   aaTCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAcc                                            <  1:2091571/100‑1 (MQ=255)
                                                   aaTCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAcc                                            >  2:588521/1‑100 (MQ=255)
                                                      ncccGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCa                                         >  2:328338/2‑100 (MQ=255)
                                                      ccccGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCa                                         >  2:1641844/1‑100 (MQ=255)
                                                      ccccGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCa                                         >  2:755154/1‑100 (MQ=255)
                                                         cGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCAc                                        <  2:813441/98‑1 (MQ=255)
                                                         cGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTc                                      <  1:533859/100‑1 (MQ=255)
                                                             aTGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACg                                  <  2:2379740/100‑1 (MQ=255)
                                                              tGCCCAGGGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTCGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGc                                 <  2:1170724/100‑1 (MQ=255)
                                                              tGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGc                                 <  1:1706616/100‑1 (MQ=255)
                                                              tGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGc                                 >  2:2391612/1‑100 (MQ=255)
                                                                cccAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGcta                               >  2:630297/1‑100 (MQ=255)
                                                                  cAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGctaac                             <  1:109784/100‑1 (MQ=255)
                                                                  cAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTTCGGCTACGTATACATCGCTGTTTACCCCACTCAACGctaac                             <  1:1188673/100‑1 (MQ=255)
                                                                       cGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAAtt                        <  1:952985/100‑1 (MQ=255)
                                                                        gTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCTCACTCAACGCTAACTAAttt                       >  2:575514/1‑100 (MQ=255)
                                                                        gTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAAtt                        <  1:333801/99‑1 (MQ=255)
                                                                        gTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAAttt                       >  2:1844465/1‑100 (MQ=255)
                                                                            cccTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACg                   >  2:839890/1‑100 (MQ=255)
                                                                            cagtGTACCGCAGCGTTTTCTGCTGGGAAGCTGGTACTCTTGGTAGGTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCAc                    <  2:2239003/96‑1 (MQ=21)
                                                                             ccTGTACCACAACGGTTTTACCTGCGAATCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTATCCCACTCAACGCTAACTAATTTCACg                   <  2:1313912/99‑1 (MQ=255)
                                                                             ccTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCActa                  >  2:1456147/1‑98 (MQ=255)
                                                                              cTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACgag                 <  1:1316100/100‑1 (MQ=255)
                                                                              cTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACgag                 >  1:2163616/1‑100 (MQ=255)
                                                                              cTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTGACCCCACTCAACGCTAACTAATTTCACga                  >  1:315303/1‑99 (MQ=255)
                                                                              cTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACgag                 >  1:1693156/1‑100 (MQ=255)
                                                                               tttACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACgaga                <  2:628276/98‑1 (MQ=255)
                                                                               tggaCCACAACGGTTTTACCTGCGAAGCCGATGCTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACgaga                <  1:1057740/97‑1 (MQ=255)
                                                                               tGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACgaga                <  1:340110/100‑1 (MQ=255)
                                                                               tGTACCACAACGGTTTCACCTGCGAAGCCTATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGGTTACCCCACTCAACGCTAACTAATTTCACgaga                <  1:345351/100‑1 (MQ=255)
                                                                                gTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGagag               <  2:1136486/100‑1 (MQ=255)
                                                                                 tACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagc              <  1:552980/100‑1 (MQ=255)
                                                                                  aCCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagca             >  2:1920510/1‑100 (MQ=255)
                                                                                  aCCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagca             >  2:472995/1‑100 (MQ=255)
                                                                                   ccACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaa            >  2:496979/1‑100 (MQ=255)
                                                                                    cacaACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaag           >  1:2039022/1‑100 (MQ=255)
                                                                                    cacaACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaag           >  1:565280/1‑100 (MQ=255)
                                                                                    cacaACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaag           >  2:1397077/1‑100 (MQ=255)
                                                                                     acaACGGTTTTACCTACGAAGCCGATACTCTTGATAGTTGCGGCTACGTATACATCGCTGTTTACCCAACTCAACGCTAACTAATTTCACGAGagcaagc          >  2:901478/1‑100 (MQ=255)
                                                                                     acaACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaagc          >  1:220459/1‑100 (MQ=255)
                                                                                      caACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaagca         <  1:1377873/100‑1 (MQ=255)
                                                                                       aaCGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaagca         >  1:280440/1‑99 (MQ=255)
                                                                                         cGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGa      >  1:677209/1‑100 (MQ=255)
                                                                                            tttcACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAac    <  1:864973/99‑1 (MQ=255)
                                                                                             tttACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAacac  <  2:2202674/100‑1 (MQ=255)
                                                                                               |                                                                                                 
TTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACAC  >  NC_000913/2777268‑2777460

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: