Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,777,363 T→C 21.5% 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,3630TC21.5% 293.0 / 83.8 158F83F (TTT→TTCyfjZCP4‑57 prophage; antitoxin of the YpjF‑YfjZ toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (60/64);  new base C (22/12);  total (82/76)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

Only 100 of 163 total aligned reads displayed.

GGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTC  >  NC_000913/2777264‑2777462
                                                                                                   |                                                                                                   
ggTATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGtttt                                                                                                     >  2:3095832/1‑100 (MQ=255)
ggTATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGtttt                                                                                                     >  1:3484738/1‑100 (MQ=255)
 gTATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTa                                                                                                    <  1:1335280/100‑1 (MQ=255)
 gTATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTa                                                                                                    <  1:1425822/100‑1 (MQ=255)
   aTTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTAcc                                                                                                  <  1:2110302/100‑1 (MQ=255)
    tttCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCt                                                                                                 >  1:3887934/1‑100 (MQ=255)
    tttCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCt                                                                                                 <  2:3077534/100‑1 (MQ=255)
    tttCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCt                                                                                                 >  2:1297020/1‑100 (MQ=255)
     ttCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTg                                                                                                >  2:3609647/1‑100 (MQ=255)
     ttCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTg                                                                                                >  2:532010/1‑100 (MQ=255)
      tCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCt                                                                                                 >  1:372553/1‑98 (MQ=255)
      tCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCATGTACCACAACGGTTTTACCTGc                                                                                               <  2:170827/100‑1 (MQ=255)
        cacaTTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGa                                                                                             >  2:3197318/1‑100 (MQ=255)
         acaTTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGaa                                                                                            >  1:1193533/1‑100 (MQ=255)
          caTTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAg                                                                                           <  2:3963539/100‑1 (MQ=255)
               tATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGAt                                                                                      <  1:924073/100‑1 (MQ=255)
                  naGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATAct                                                                                   >  2:844227/2‑100 (MQ=255)
                  cAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCGATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGAtcc                                                                                    >  2:2036265/1‑97 (MQ=255)
                   aGCAAGATAGAGTCGATACCGAGCACTGGTGAACTGCATCCCCGTCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGAt                                                                                      <  2:1253596/96‑1 (MQ=25)
                    gCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATActct                                                                                 >  2:961902/1‑100 (MQ=255)
                    gCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATActct                                                                                 >  2:3848386/1‑100 (MQ=255)
                    gCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATActct                                                                                 >  1:1019091/1‑100 (MQ=255)
                       aGATAGAGTGGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTgg                                                                              <  1:371695/100‑1 (MQ=255)
                        gATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGt                                                                             <  1:3911607/100‑1 (MQ=255)
                         aTAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGt                                                                             <  2:1493301/99‑1 (MQ=255)
                          tAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAg                                                                           <  2:3155292/100‑1 (MQ=255)
                             agTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTg                                                                        >  2:2379413/1‑100 (MQ=255)
                              gTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGc                                                                       >  2:794336/1‑100 (MQ=255)
                               tCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCg                                                                      <  2:3090057/100‑1 (MQ=255)
                               tCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGGCACCCTGTACCACAAGGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCg                                                                      <  2:1565049/100‑1 (MQ=255)
                                 gATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGc                                                                    <  2:4049162/100‑1 (MQ=255)
                                 gATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGc                                                                    <  1:1883559/100‑1 (MQ=255)
                                 gATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGc                                                                    >  2:1998488/1‑100 (MQ=255)
                                 gATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGc                                                                    >  2:420862/1‑100 (MQ=255)
                                     cTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACg                                                                >  2:3831589/1‑100 (MQ=255)
                                      tGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACACCGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGt                                                               >  2:2855542/1‑100 (MQ=255)
                                       gACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGta                                                              <  1:2184650/100‑1 (MQ=255)
                                       gACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGta                                                              <  1:867603/100‑1 (MQ=255)
                                        aCCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGtat                                                             <  1:1523998/100‑1 (MQ=255)
                                            cTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACAt                                                         <  1:444603/100‑1 (MQ=255)
                                              ggTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCg                                                       >  2:637189/1‑100 (MQ=255)
                                               gTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGc                                                      <  1:2564015/100‑1 (MQ=255)
                                               gTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGc                                                      <  1:338307/100‑1 (MQ=255)
                                               gTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGc                                                      <  2:3013076/100‑1 (MQ=255)
                                                tGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCt                                                     >  2:3268946/1‑100 (MQ=255)
                                                   aCTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATCCTCTTGGTAGTTGCGGCTACGTATACATCGCTGt                                                   >  1:737786/1‑99 (MQ=255)
                                                   aCTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGtt                                                  >  1:1204047/1‑100 (MQ=255)
                                                    cTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACGACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGttt                                                 <  1:1925789/100‑1 (MQ=255)
                                                    cTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGttt                                                 >  1:133597/1‑100 (MQ=255)
                                                     tGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATAGCTGTTTa                                                >  1:1850136/1‑100 (MQ=255)
                                                      gAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAc                                               >  1:1020501/1‑100 (MQ=255)
                                                      gAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAc                                               <  1:3075833/100‑1 (MQ=255)
                                                      gAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAc                                               >  1:455631/1‑100 (MQ=255)
                                                        aacaccGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAccc                                             <  2:359581/96‑1 (MQ=255)
                                                        aTCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAccc                                             >  2:30693/1‑100 (MQ=255)
                                                        aTCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAccc                                             >  1:196005/1‑100 (MQ=255)
                                                         tCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAcccc                                            >  2:114457/1‑100 (MQ=255)
                                                         tCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAcccc                                            >  2:3564100/1‑100 (MQ=255)
                                                          ccccGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCa                                           >  1:665736/1‑100 (MQ=255)
                                                           cccGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCAc                                          <  2:2356484/100‑1 (MQ=255)
                                                           cccGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCAc                                          >  2:1529856/1‑100 (MQ=255)
                                                            ccGCGATGCGCAATGCGGCACCCGGTAGCACATCGGGCTTACCTGCGAAGCCGTTTCTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAccc                                             <  2:634750/96‑1 (MQ=21)
                                                             cGCCATGCCCAATGCGTCACCCTGTACCGCAACGGTTTTACCTGCGAAGCCGTTACTCTTGGTAGTGGCGGCTACGTATACATCGCTGTTTACCCCACTc                                        <  1:481758/100‑1 (MQ=255)
                                                             cGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTc                                        <  2:2697845/100‑1 (MQ=255)
                                                             cGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTc                                        >  2:1921949/1‑100 (MQ=255)
                                                               ccATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCaa                                      >  1:270853/1‑100 (MQ=255)
                                                                  tGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGc                                   <  2:1485029/100‑1 (MQ=255)
                                                                   gCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGCTACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGct                                  >  2:2296843/1‑100 (MQ=255)
                                                                    cccAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGcta                                 <  2:2886746/100‑1 (MQ=255)
                                                                    cccAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGcta                                 >  2:3528487/1‑100 (MQ=255)
                                                                    cccAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGcta                                 >  1:3608959/1‑100 (MQ=255)
                                                                    cccAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGcta                                 >  2:124087/1‑100 (MQ=255)
                                                                    cccAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGcta                                 >  1:4054755/1‑100 (MQ=255)
                                                                     ccAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGctaa                                >  1:1819483/1‑100 (MQ=255)
                                                                      cAATGCCTCACCCTGTACCACAACTGTTTTACCTTCTAATCCCATACTCTTGGTAGTTCCGGCTACGTATACATCGCTGTTTACCCCACTCAACGctaac                               <  1:2377765/100‑1 (MQ=21)
                                                                       aaTGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGctaact                              <  2:49709/100‑1 (MQ=255)
                                                                       aaTGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGctaact                              <  1:2199114/100‑1 (MQ=255)
                                                                            gTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAAttt                         <  1:2948832/100‑1 (MQ=255)
                                                                            gTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAAttt                         >  2:1673717/1‑100 (MQ=255)
                                                                             tCACTCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTc                        <  2:1766771/100‑1 (MQ=255)
                                                                             tCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTc                        <  1:1181448/100‑1 (MQ=255)
                                                                             tCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTc                        >  2:3011252/1‑100 (MQ=255)
                                                                             tCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTc                        <  1:86710/100‑1 (MQ=255)
                                                                               aCCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCAc                      <  1:3770396/100‑1 (MQ=255)
                                                                                cccTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGCATACATCGCTGTTTACCCCACTAAACGCTAACTAATTTCACg                     >  2:2761165/1‑100 (MQ=255)
                                                                                 ccTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACg                     <  1:2318614/99‑1 (MQ=255)
                                                                                    gTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACgaga                  >  2:3693297/1‑99 (MQ=255)
                                                                                    gTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGagag                 <  1:161517/100‑1 (MQ=255)
                                                                                     tACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagc                >  1:83921/1‑100 (MQ=255)
                                                                                      aCCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagca               <  2:2867328/100‑1 (MQ=255)
                                                                                        cacaACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaag             >  1:1923167/1‑100 (MQ=255)
                                                                                        cacaACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaag             >  2:860176/1‑100 (MQ=255)
                                                                                        cacaACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaag             >  2:1704881/1‑100 (MQ=255)
                                                                                        caaaaCGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaag             <  2:2588307/97‑1 (MQ=255)
                                                                                         acaACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaagc            >  2:2382359/1‑100 (MQ=255)
                                                                                         acaACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaagc            <  1:2556323/100‑1 (MQ=255)
                                                                                          caACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagaaagca           >  2:510673/1‑100 (MQ=255)
                                                                                          caACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaagca           <  1:488068/100‑1 (MQ=255)
                                                                                                  ttACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACAct   >  1:955909/1‑100 (MQ=255)
                                                                                                   tACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACActc  <  1:579479/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTC  >  NC_000913/2777264‑2777462

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: