Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,777,363 T→C 100% 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,3630TC99.2% 426.1 / ‑4.8 127F83F (TTT→TTCyfjZCP4‑57 prophage; antitoxin of the YpjF‑YfjZ toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (63/63);  total (64/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 = 8.50e-01

Only 100 of 127 total aligned reads displayed.

TATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAAC  >  NC_000913/2777266‑2777458
                                                                                                 |                                                                                               
 aTTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCAcc                                                                                              >  1:317174/1‑100 (MQ=255)
  tttCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCAc                                                                                               >  1:1319807/1‑98 (MQ=255)
  tttCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCt                                                                                             <  2:1515686/100‑1 (MQ=255)
   ttCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTg                                                                                            <  2:2366730/100‑1 (MQ=255)
   ttCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAAATGAATCCCCTCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTg                                                                                            <  2:248775/100‑1 (MQ=255)
    tCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGc                                                                                           >  1:2665862/1‑100 (MQ=255)
    tCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGc                                                                                           <  2:2616464/100‑1 (MQ=255)
     ccACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCg                                                                                          >  2:503982/1‑100 (MQ=255)
     ccACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAACGCGTCACCCTGTACCACAACGGTTTCACCTGCg                                                                                          <  1:2258853/100‑1 (MQ=255)
      cacaTTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCg                                                                                          <  2:1399033/99‑1 (MQ=255)
      cacaTTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGa                                                                                         >  2:2082553/1‑100 (MQ=255)
      cacaTTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGa                                                                                         >  1:603889/1‑100 (MQ=255)
      cacaTTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGa                                                                                         <  2:2085913/100‑1 (MQ=255)
        caTTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAg                                                                                       >  2:351999/1‑100 (MQ=255)
           tttATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCg                                                                                    <  1:1922988/100‑1 (MQ=255)
            ttATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGATTCACCTGCGAAGCCGa                                                                                   >  1:1383884/1‑100 (MQ=255)
                 agacAGATAGATTCGATTCTAACCACTGGTGAAATAAATCCAATCCATGAAAAATGCGTCACCCTGTACAAAAACGGTTTCACCTGCGAAGCCGAtttt                                                                               >  2:346048/4‑96 (MQ=12)
                 aGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATActc                                                                              >  1:1880716/1‑100 (MQ=255)
                 aGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATActc                                                                              <  1:2211836/100‑1 (MQ=255)
                  gCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATActct                                                                             >  1:757529/1‑100 (MQ=255)
                     aGATAGAGTCGATGCTGACCACTGGTGAACTGAATCTCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTgg                                                                          <  1:1193878/100‑1 (MQ=255)
                      gATAGAGTCGATGCTGACCACTTGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGt                                                                         <  1:1219116/100‑1 (MQ=255)
                      gATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGt                                                                         <  1:906598/100‑1 (MQ=255)
                       aTAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTa                                                                        <  2:2087668/100‑1 (MQ=255)
                         agagTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGt                                                                      <  2:2552705/100‑1 (MQ=255)
                          gagTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCATTTGCGAAGCCGATACTCTTg                                                                           <  2:594109/94‑1 (MQ=255)
                          gagTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGt                                                                      <  2:757388/99‑1 (MQ=255)
                            gTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTg                                                                    >  1:388885/1‑99 (MQ=255)
                            gTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGc                                                                   >  2:2284035/1‑100 (MQ=255)
                            gTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGc                                                                   >  1:1398755/1‑100 (MQ=255)
                            gTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGc                                                                   <  2:571794/100‑1 (MQ=255)
                             tCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCg                                                                  <  1:1570605/100‑1 (MQ=255)
                             tCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCg                                                                  >  1:2451327/1‑100 (MQ=255)
                              cGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCgg                                                                 >  1:923664/1‑100 (MQ=255)
                               gATGCTGACCACTGGTGAGCTGAATCCCCGGCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGc                                                                <  2:809597/100‑1 (MQ=255)
                               gATGCTGACCACTGGTGAACTGAATCCCCGCCATGCTCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGc                                                                <  1:2136644/100‑1 (MQ=255)
                               gATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGc                                                                <  1:1314918/100‑1 (MQ=255)
                                 tGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTa                                                              >  2:2817691/1‑100 (MQ=255)
                                 tGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTa                                                              <  1:2455486/100‑1 (MQ=255)
                                   ttGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACg                                                            <  1:2834121/99‑1 (MQ=255)
                                   cTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACg                                                            >  1:2598552/1‑100 (MQ=255)
                                    tGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACga                                                           >  2:2834089/1‑99 (MQ=255)
                                    tGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGt                                                           >  2:2971442/1‑100 (MQ=255)
                                    tGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGt                                                           <  2:2564869/100‑1 (MQ=255)
                                      aCCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGtat                                                         <  1:52550/100‑1 (MQ=255)
                                      aCCAATGGTGAACTGAATCCCCGCCATGCCCATTGCGTCACCCTTTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGtat                                                         <  1:51845/100‑1 (MQ=21)
                                       ccACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGtata                                                        >  1:184279/1‑100 (MQ=255)
                                       ccACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGtat                                                         <  2:1855403/99‑1 (MQ=255)
                                            ggTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATc                                                    >  2:4284/1‑99 (MQ=255)
                                             gTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGc                                                  <  1:2941226/100‑1 (MQ=255)
                                               gAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAACCGATACTCTTGATAGATACAGCTACGTATACATc                                                    >  1:92484/1‑96 (MQ=21)
                                                aaCTGAATCCCCGCCATGCCCAATGCGTCACTCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGt                                               <  1:2501948/100‑1 (MQ=255)
                                                 aCTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGtt                                              >  1:607209/1‑100 (MQ=255)
                                                 aCTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGtt                                              <  2:2878421/100‑1 (MQ=255)
                                                   tGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGttt                                             <  2:2330238/99‑1 (MQ=255)
                                                    gAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAc                                           >  1:2174945/1‑100 (MQ=255)
                                                    gAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAc                                           >  1:2433595/1‑100 (MQ=255)
                                                    gAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCCACGTATACATCGCTGTTTAc                                           >  1:1798706/1‑100 (MQ=255)
                                                     aaTCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAcc                                          <  2:1217119/100‑1 (MQ=255)
                                                      aTCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAccc                                         <  1:2267474/100‑1 (MQ=255)
                                                      aTCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAccc                                         <  1:827043/100‑1 (MQ=255)
                                                      aTCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAccc                                         <  1:1133868/100‑1 (MQ=255)
                                                       tCCTCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAcccc                                        <  2:2625820/100‑1 (MQ=255)
                                                       tCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACTTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAcccc                                        <  1:2425320/100‑1 (MQ=255)
                                                       tCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAcccc                                        <  1:1687309/100‑1 (MQ=255)
                                                         cccGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCAc                                      >  1:1611413/1‑100 (MQ=255)
                                                          ccGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACt                                     <  1:2950927/100‑1 (MQ=255)
                                                          ccGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACt                                     >  2:1795849/1‑100 (MQ=255)
                                                            gCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCa                                   >  2:2820549/1‑100 (MQ=255)
                                                             ccATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGGAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCaa                                  >  2:3005712/1‑100 (MQ=255)
                                                              cATGCCCAATGCGTCACCCTGTCCCACAACGTTTTCAGCTGCGAAGCGGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAAc                                 <  2:2046332/100‑1 (MQ=18)
                                                               attcccAAAGCGTAACCCGGTACCACAACGTTTTCACCTGCGCTGCCGATCCTCTTGGTAGTTGCGGCTACGTATTCATCGCTGTTCACCCCACTGAACg                                <  2:1699027/97‑1 (MQ=17)
                                                                   ccAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGctaa                            >  1:1318916/1‑100 (MQ=255)
                                                                   ccAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGctaa                            >  2:84237/1‑100 (MQ=255)
                                                                   ccAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGctaa                            >  1:3001233/1‑100 (MQ=255)
                                                                   ccAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGctaa                            >  1:1617562/1‑100 (MQ=255)
                                                                         cGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAAtt                      >  2:872766/1‑100 (MQ=255)
                                                                         cGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAAtt                      <  1:224990/100‑1 (MQ=255)
                                                                          gTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAAttt                     >  2:2301675/1‑100 (MQ=255)
                                                                          gTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAAttt                     <  2:852205/100‑1 (MQ=255)
                                                                           tCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTc                    <  2:757525/100‑1 (MQ=255)
                                                                            cACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCa                   >  2:494977/1‑100 (MQ=255)
                                                                             aCCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCAc                  <  1:503985/100‑1 (MQ=255)
                                                                             aCCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCAc                  <  2:184277/100‑1 (MQ=255)
                                                                              cccTGTACCACAACGTTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACTCTAACTAATTTCACg                 <  2:1773809/100‑1 (MQ=255)
                                                                              cccTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACg                 >  2:1691633/1‑100 (MQ=255)
                                                                                cTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACgag               >  1:1284156/1‑100 (MQ=255)
                                                                                 tGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACgaga              <  2:2631813/100‑1 (MQ=255)
                                                                                 tGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACgaga              >  2:2231446/1‑100 (MQ=255)
                                                                                  gTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGagag             >  1:734791/1‑100 (MQ=255)
                                                                                    aCCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagca           >  1:757554/1‑100 (MQ=255)
                                                                                     ccACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaa          >  2:456437/1‑100 (MQ=255)
                                                                                      cacaACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaag         >  2:1296266/1‑100 (MQ=255)
                                                                                       acaACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaagc        >  1:2762478/1‑100 (MQ=255)
                                                                                         aaCGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCCAACTAATTTCACGAGagcaagc        <  2:1757574/98‑1 (MQ=255)
                                                                                          aCGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATg     <  2:594853/100‑1 (MQ=255)
                                                                                          aCGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATg     >  1:2124983/1‑100 (MQ=255)
                                                                                           cGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGa    <  2:734787/100‑1 (MQ=255)
                                                                                            ggTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACACTAACTAATTTCACGAGAGCAAGCATGaa   >  1:1197987/1‑100 (MQ=255)
                                                                                             gTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAac  <  1:1568568/100‑1 (MQ=255)
                                                                                                 |                                                                                               
TATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAAC  >  NC_000913/2777266‑2777458

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