Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,777,363 T→C 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,3630TC98.4% 426.6 / ‑3.2 129F83F (TTT→TTCyfjZCP4‑57 prophage; antitoxin of the YpjF‑YfjZ toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (69/57);  total (72/57)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

Only 100 of 129 total aligned reads displayed.

TATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATG  >  NC_000913/2777266‑2777455
                                                                                                 |                                                                                            
tATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCAc                                                                                            >  1:365108/1‑100 (MQ=255)
 aTTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCAcc                                                                                           <  2:1341658/100‑1 (MQ=255)
 aTTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCAcc                                                                                           >  1:2863047/1‑100 (MQ=255)
    tCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCTCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGc                                                                                        >  1:2228521/1‑100 (MQ=255)
    tCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGc                                                                                        >  2:2420796/1‑100 (MQ=255)
    tCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGc                                                                                        >  1:3176144/1‑100 (MQ=255)
    tCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGc                                                                                        >  1:2398664/1‑100 (MQ=255)
    tCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGc                                                                                        >  1:2303075/1‑100 (MQ=255)
     ccACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCg                                                                                       >  1:452626/1‑100 (MQ=255)
     ccACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCg                                                                                       >  2:2078378/1‑100 (MQ=255)
     ccACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCg                                                                                       >  2:1938217/1‑100 (MQ=255)
     ccACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCg                                                                                       <  1:3097548/100‑1 (MQ=255)
     ccACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCg                                                                                       >  1:262194/1‑100 (MQ=255)
     ccACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCg                                                                                       >  1:1891877/1‑100 (MQ=255)
      cacaTTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGa                                                                                      >  2:858407/1‑100 (MQ=255)
        catgtgagcAGCCCGATAGAGTCGATGCTGACCGCTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAg                                                                                    <  1:2822903/92‑1 (MQ=21)
        caTTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAg                                                                                    >  2:575966/1‑100 (MQ=255)
         aTTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGc                                                                                   >  1:1647770/1‑100 (MQ=255)
          ttttATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGcc                                                                                  <  2:1538525/100‑1 (MQ=255)
            ttATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGa                                                                                >  2:1869685/1‑100 (MQ=255)
             tATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGAt                                                                               <  2:1203311/100‑1 (MQ=255)
             tATCAGCCAGATAGAGTCGATGCTGACCACTGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGAt                                                                               >  1:962525/1‑99 (MQ=255)
               tCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATAc                                                                             <  2:2319916/100‑1 (MQ=255)
               tCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATAc                                                                             <  1:2755397/100‑1 (MQ=255)
                 aGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATActc                                                                           >  1:1340657/1‑100 (MQ=255)
                 aGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATActc                                                                           >  1:3209356/1‑100 (MQ=255)
                    cAGTTAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTg                                                                        <  2:3443299/100‑1 (MQ=255)
                    cAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTg                                                                        <  2:1117099/100‑1 (MQ=255)
                      gATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGt                                                                      <  1:278428/100‑1 (MQ=255)
                      gATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGt                                                                      <  2:3176196/100‑1 (MQ=255)
                       tccgatctgATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACGACGGTTTCACCTGCGAAGCCGATACTCTTGGTa                                                                     <  2:453739/92‑1 (MQ=255)
                        tAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAg                                                                    >  1:2175681/1‑100 (MQ=255)
                        tAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAg                                                                    <  1:2803673/100‑1 (MQ=255)
                         agagTCGCTGCCGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAGCGGTTTCACCTGCGAAGCCGATACTCTTGGTAGt                                                                   >  1:1827545/1‑100 (MQ=255)
                         agagTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGt                                                                   >  1:3106684/1‑100 (MQ=255)
                          gagTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCATGTACCACAACGGTTTGACCTGCGAAGCCGATACTCTTGGTAGtt                                                                  >  1:1355347/1‑100 (MQ=255)
                            gTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGc                                                                >  2:1563806/1‑100 (MQ=255)
                             tCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCg                                                               <  2:2511240/100‑1 (MQ=255)
                             ccgatcttGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCg                                                               <  1:1225274/93‑1 (MQ=25)
                              cGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCgg                                                              >  2:2322707/1‑100 (MQ=255)
                               gATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACGACGGTTTCACCTGCGAAGCCGATACTCTTGGTaagatcgga                                                             >  1:453734/1‑92 (MQ=255)
                                 tGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTa                                                           <  1:355546/100‑1 (MQ=255)
                                   cTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACg                                                         <  2:1768120/100‑1 (MQ=255)
                                   cTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACg                                                         <  2:1034141/100‑1 (MQ=255)
                                    tGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCgagatcgg                                                        >  2:1225289/1‑93 (MQ=25)
                                     gACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGTCTACGta                                                       <  2:2649085/100‑1 (MQ=255)
                                      aCCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGtat                                                      >  1:3243540/1‑100 (MQ=255)
                                       ccACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTTCGAAGCCGATTCTCTTGGTAGTTGCGGCTACGtata                                                     >  1:1286/1‑100 (MQ=255)
                                       ccACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGtata                                                     <  2:3130798/100‑1 (MQ=255)
                                       ccACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGtata                                                     >  1:2586752/1‑100 (MQ=255)
                                        cACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATAc                                                    >  2:2136382/1‑100 (MQ=255)
                                        cACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATAc                                                    >  2:718241/1‑100 (MQ=255)
                                         aCTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACa                                                   >  2:2632199/1‑100 (MQ=255)
                                            ggTGAACTGAATCCCCGCCATGCCCAATGTGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCg                                                >  2:3157689/1‑100 (MQ=255)
                                             gTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGc                                               <  2:1910290/100‑1 (MQ=255)
                                               gAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTg                                             <  2:235520/100‑1 (MQ=255)
                                               gAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTg                                             >  2:1243245/1‑100 (MQ=255)
                                               gAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTg                                             >  1:2038610/1‑100 (MQ=255)
                                                aaCTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTg                                             <  1:317434/99‑1 (MQ=255)
                                                aaCTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGt                                            <  2:2070498/100‑1 (MQ=255)
                                                 aCTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGtt                                           <  1:2719310/100‑1 (MQ=255)
                                                   tGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTa                                         <  1:2853854/100‑1 (MQ=255)
                                                    gAATCCCCGCCATACCCAATGCGTCACATTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTATTTGCGGCTAGGTATACATCTCTGTTTAc                                        <  1:283247/100‑1 (MQ=255)
                                                     aaTCCCCGCCTTGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTCGTCGTTGCGGCTACGTATACATCGCTGTTTAcc                                       <  1:2285007/100‑1 (MQ=25)
                                                     aaTCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAcc                                       <  1:1563784/100‑1 (MQ=255)
                                                      aTCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAccc                                      <  1:1674446/100‑1 (MQ=255)
                                                      aTCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAccc                                      >  2:2879506/1‑100 (MQ=255)
                                                      aTCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAccc                                      >  2:2424278/1‑100 (MQ=255)
                                                      aTCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAccc                                      >  2:3362788/1‑100 (MQ=255)
                                                      aTCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAccc                                      >  2:198805/1‑100 (MQ=255)
                                                       tCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAcccc                                     >  1:408608/1‑100 (MQ=255)
                                                       tCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGGTTAcccc                                     <  2:1827570/100‑1 (MQ=255)
                                                        ccccGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCa                                    >  2:1900433/1‑100 (MQ=255)
                                                         cccGCCATGCCCAATGCGTCACCCTGTACCACAACGGCTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCAc                                   >  1:1608272/1‑100 (MQ=255)
                                                              cATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAAc                              >  2:3070211/1‑100 (MQ=255)
                                                               aTGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACg                             <  1:3465347/100‑1 (MQ=255)
                                                                 gCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGct                           >  2:3396710/1‑100 (MQ=255)
                                                                  ccccATGCGTGCCCCTTTCCCAAAGCGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAccc                                      <  1:1210145/88‑1 (MQ=18)
                                                                  cccAATGCGTCACCCTGTACCACAACGGCTTCACCTCCGAAGCCGCTACTCTCGCTAGTCGCGCCTACCTATACACCCCCGTTCACCCCCCCCAAccc                            >  1:3147144/1‑96 (MQ=12)
                                                                    cAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGctaac                        >  1:496932/1‑100 (MQ=255)
                                                                        gCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAAt                    <  1:2009675/100‑1 (MQ=255)
                                                                         cGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAAtt                   <  2:549227/100‑1 (MQ=255)
                                                                         cGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAAtt                   >  1:1340884/1‑100 (MQ=255)
                                                                           tCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTc                 <  1:718232/100‑1 (MQ=255)
                                                                            cACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCa                >  1:3379100/1‑100 (MQ=255)
                                                                             aCCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCAc               >  2:3060521/1‑100 (MQ=255)
                                                                                cTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACgag            <  1:376909/100‑1 (MQ=255)
                                                                                cTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCGTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACgag            >  1:1860242/1‑100 (MQ=255)
                                                                                  gTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGagag          <  1:606289/100‑1 (MQ=255)
                                                                                    aCCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagca        >  2:25813/1‑100 (MQ=255)
                                                                                     ccACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaa       <  1:2668837/100‑1 (MQ=255)
                                                                                     ccACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaa       >  2:1184678/1‑100 (MQ=255)
                                                                                      cacaACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaag      >  1:627585/1‑100 (MQ=255)
                                                                                      cacaACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaag      >  2:3394623/1‑100 (MQ=255)
                                                                                      cacaACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaag      >  2:2060261/1‑100 (MQ=255)
                                                                                       acaACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaagc     >  2:1389166/1‑100 (MQ=255)
                                                                                       acaACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaagc     >  2:650460/1‑100 (MQ=255)
                                                                                        caACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaagca    <  2:2812110/100‑1 (MQ=255)
                                                                                        caACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaagca    <  1:3158995/100‑1 (MQ=255)
                                                                                          aCGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATg  <  1:1218558/100‑1 (MQ=255)
                                                                                                 |                                                                                            
TATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATG  >  NC_000913/2777266‑2777455

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: