Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,597,272 T→A 18.9% pseudogene (716/3861 nt) yneO ← pseudogene, AidA homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,597,2720TA18.9% 279.6 / 59.7 126pseudogene (716/3861 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base T (57/45);  new base A (12/12);  total (69/57)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-02

Only 100 of 128 total aligned reads displayed.

ATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAA  >  NC_000913/1597173‑1597370
                                                                                                   |                                                                                                  
atatTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAt                                                                                                    <  2:3789059/100‑1 (MQ=255)
  atTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTg                                                                                                  <  1:1392393/100‑1 (MQ=255)
  atTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACAGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTg                                                                                                  <  1:1281297/100‑1 (MQ=255)
  atTGACTGTGCCATTACCGAAATATCAATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTg                                                                                                  >  2:294584/1‑100 (MQ=255)
    tGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGcc                                                                                                >  1:1413798/1‑100 (MQ=255)
     gACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCCGGGTTCCGATACCGCTGCCATTGcc                                                                                                >  2:1513210/1‑99 (MQ=255)
      aCTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCAc                                                                                              <  1:2695654/100‑1 (MQ=255)
      aCTGAGCCATTACCGAAATATCCATTATAGAGTCTGAGTACATTGATGACAGAGTCTTGATCCTGCACATTCAGTGTTACGATACCGGTGGTATTg                                                                                                  <  2:1021839/96‑1 (MQ=255)
        tgtgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCg                                                                                            >  2:3540347/1‑100 (MQ=255)
          tgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCggt                                                                                          <  1:845764/100‑3 (MQ=255)
          tgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCggt                                                                                          <  1:2962107/100‑3 (MQ=255)
          tgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCggt                                                                                          <  1:2932377/100‑3 (MQ=255)
           gCCATTACCGAATTATCCATTATAGAGTCTGCGTAATGTGATGATTTAGTCTTGATCCTGCACATTCAGGGTGCCGATACCGGTGGCAATGCCACCGGTa                                                                                         <  1:694212/100‑1 (MQ=255)
           gCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATATCGGTGGCAATGCCACCGGTa                                                                                         <  1:1035559/100‑1 (MQ=255)
           gCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATa                                                                                         <  2:1830240/100‑1 (MQ=255)
            ccATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCTACCGATaa                                                                                        >  2:73095/1‑100 (MQ=255)
            ccATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATaa                                                                                        <  2:2069667/100‑1 (MQ=255)
            ccATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATaa                                                                                        <  2:1726276/100‑1 (MQ=255)
             cATTANCTAAATATCCATTATAGAGTTTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGGTCCGATACCGGTGGCAATGCCACCGGTat                                                                                        <  1:3212245/99‑2 (MQ=255)
             cATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCGGTatc                                                                                       <  2:1178575/100‑3 (MQ=255)
              aTTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATTACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCGGTatc                                                                                       <  1:1342421/99‑3 (MQ=255)
              aTTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCggtatcn                                                                                      <  2:528063/100‑7 (MQ=255)
              aTTACCGAAATACCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGTTTCCGATACCGGTGGCAATGCCACCGGTatcg                                                                                      <  2:709634/100‑4 (MQ=255)
               ttACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCg                                                                                     <  2:1292752/100‑1 (MQ=255)
               ttACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCg                                                                                     <  2:3329969/100‑1 (MQ=255)
               ttACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCg                                                                                     <  2:1191753/100‑1 (MQ=255)
                tacctaataatCCATTAAAAAGTCTGCGTAACTTGATGACAGAGTCTTTATCCTGCACATTCCGGGTTTCGATTCCTGTGGCATTGCTACCGATAATCg                                                                                     <  2:1825617/90‑1 (MQ=255)
                tACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCga                                                                                    >  1:1292766/1‑99 (MQ=255)
                tACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGt                                                                                    >  2:1098324/1‑100 (MQ=255)
                 aCCGAAATATCCATTATAGAGTCTGCGTACCATGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTa                                                                                   >  2:3054221/1‑100 (MQ=255)
                 aCCAAAATATCCATTATAGATTCTCCGTACCGTGATGACACAGTCTTTATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCAACCGGTatcgg                                                                                     <  2:224430/98‑5 (MQ=255)
                  ccGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTTATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAt                                                                                  <  2:1194968/100‑1 (MQ=255)
                  ccGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTa                                                                                   >  2:1460855/1‑99 (MQ=255)
                  ccGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATAGTAt                                                                                  >  1:758923/1‑100 (MQ=255)
                    gAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAttt                                                                                >  1:2919067/1‑100 (MQ=255)
                     aaaTATCCATTAGAGAGTCTTCGGACCGTTATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGCTACACGTGGCCTTGCTCCTGATCATCGtccttt                                                                               <  2:2679849/100‑6 (MQ=255)
                      aaTATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTc                                                                              >  2:952490/1‑100 (MQ=255)
                         atCGATTATAGAGTTTGCGTACCGTGCTGACAGAGTCTTGATCCTGCACATTTAGGGTTCCGATACCGGTGGCATTGCTACCGATAATCGTATTTTCAGc                                                                           <  2:3186164/100‑1 (MQ=255)
                         atCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGc                                                                           <  1:3682009/100‑1 (MQ=255)
                         atCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGc                                                                           >  2:1462933/1‑100 (MQ=255)
                           ccATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCg                                                                         >  1:1099320/1‑100 (MQ=255)
                             aTTATAGAGTCTGCGTACCGTGATGAGAGAGTCTTGATCATAAACATTAAAGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGACg                                                                         >  2:3593763/1‑98 (MQ=255)
                             aTTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCAGTGGCATTGCCACCGATAATAGTATTTTCAGcc                                                                          >  2:3593411/1‑97 (MQ=255)
                             aTTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATACTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTa                                                                       >  2:3137424/1‑100 (MQ=255)
                                 taGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaacca                                                                   >  1:86042/1‑100 (MQ=255)
                                     gTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAccc                                                               <  1:1716043/100‑1 (MQ=255)
                                          cGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCt                                                          >  1:1174270/1‑100 (MQ=255)
                                            tACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTcg                                                        >  1:2002511/1‑100 (MQ=255)
                                              ccGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTcgcg                                                      <  1:2187752/100‑1 (MQ=255)
                                               cGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGa                                                     >  2:2811839/1‑100 (MQ=255)
                                                 tgatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAAt                                                   >  1:3780301/1‑100 (MQ=255)
                                                    tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGt                                                <  2:1099309/100‑1 (MQ=255)
                                                    tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCAATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGt                                                >  1:102781/1‑100 (MQ=255)
                                                       cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAATCCCCCCTCGCGAATAGTCGc                                             >  1:890251/1‑100 (MQ=255)
                                                       cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGc                                             >  1:1027593/1‑100 (MQ=255)
                                                       cAGAGTCTTGATCCTGCACACTCAGGGTTCCGATACCGGTGGCATTGCCAACGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGaa                                                    >  2:2827124/1‑93 (MQ=255)
                                                        agagTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGcc                                            <  1:449066/100‑1 (MQ=255)
                                                         gagTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCt                                           >  2:587443/1‑100 (MQ=255)
                                                         gagTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCCCGAATAGTCGc                                             >  1:2307316/1‑98 (MQ=255)
                                                           gTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTcc                                         <  2:751914/100‑1 (MQ=255)
                                                             cTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTcccc                                       <  2:2609853/100‑1 (MQ=255)
                                                              ttGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCTCCCCTCGCGAATAGTCGCCTcccc                                       <  2:791070/99‑1 (MQ=255)
                                                               tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAg                                     >  2:3784639/1‑100 (MQ=255)
                                                                   ccTGCACATTCAGGTTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCTAATAGTCGCCTCCCCAGTTcc                                 <  2:1501696/100‑1 (MQ=255)
                                                                   ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTcc                                 >  1:1144264/1‑100 (MQ=255)
                                                                   ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTcc                                 >  1:360906/1‑100 (MQ=255)
                                                                   ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAATTc                                  >  1:782002/1‑99 (MQ=255)
                                                                    cTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCt                                >  1:3845103/1‑100 (MQ=255)
                                                                     tGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCtt                               >  2:682837/1‑100 (MQ=255)
                                                                       cacaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGa                             <  1:2910000/100‑1 (MQ=255)
                                                                         caTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                           <  1:740195/100‑1 (MQ=255)
                                                                         caTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCCTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                           <  2:2160546/100‑1 (MQ=255)
                                                                          aTTCTGGGTTCCGATACCGGTGGCATTGCCATCGATAATCGTATTTTCTTCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTa                          <  1:986272/100‑1 (MQ=255)
                                                                          aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTa                          <  2:2373692/100‑1 (MQ=255)
                                                                             catggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCa                       <  1:1322892/97‑1 (MQ=255)
                                                                             cAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCa                       >  2:9761/1‑100 (MQ=255)
                                                                              aGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCaa                      >  1:3479688/1‑100 (MQ=255)
                                                                                ggTTCCGATACCGGTGTCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAtt                    <  2:2582175/100‑1 (MQ=255)
                                                                                  tacgattatcGGTGGCAATGCCACCGGTAATCGTATTTTCAGCAGTAACCAACCCCCCCTCGCGAATAGTCTCCTCCCCAGTTCCTTGATTAACAAtttt                  >  1:1439329/10‑99 (MQ=37)
                                                                                   tCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGa                 <  2:2943412/100‑1 (MQ=255)
                                                                                   acgattatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAATAACCTCCCCAGTTCCTTGATTACCAATTTGa                 >  1:3372608/9‑100 (MQ=255)
                                                                                    ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCCCCCCAGTTCCTTGATTACCAATTTGaa                >  2:507794/1‑100 (MQ=255)
                                                                                     gattatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa               >  2:1533010/7‑100 (MQ=255)
                                                                                      gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt              <  1:2138018/100‑1 (MQ=255)
                                                                                      attatcGGTGGCAATGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCAAATAATCAAATCCCCAATTCCTTGATTACCAATTTGAAAt              >  1:2160787/6‑100 (MQ=255)
                                                                                      attatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt              >  1:1165576/6‑100 (MQ=255)
                                                                                        tatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCACCTCCCCAGTTCCTTGATTACCAATTTGAAAtt             >  1:474659/4‑99 (MQ=255)
                                                                                         atcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa           >  1:2736096/3‑100 (MQ=255)
                                                                                         aCCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa           >  1:1889332/1‑100 (MQ=255)
                                                                                          ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCaa          >  2:1094753/1‑100 (MQ=255)
                                                                                           cGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCaa          <  1:3166399/99‑1 (MQ=255)
                                                                                           cGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAt         <  1:3912830/100‑1 (MQ=255)
                                                                                            ngTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTAATTACCAATTTGAAATTCAAtt        >  1:644799/2‑100 (MQ=255)
                                                                                            ggTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt        >  2:1201708/1‑100 (MQ=255)
                                                                                            ggTGGCAATGCCACTGATAATCGTATTTTCAGCAGTAATCAACCCCAACTCGCGAATAGTCGCCTCCCCACTTCCTTGATTACCAATTTGAAATTCaaat        >  2:213452/1‑98 (MQ=255)
                                                                                            ggTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCACATCCCCAGTTCCTTGATTACCAATTCGAAATACAAtt        >  2:3818973/1‑100 (MQ=255)
                                                                                             gTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt        >  2:468836/1‑99 (MQ=255)
                                                                                              tGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTg       <  1:935611/99‑1 (MQ=255)
                                                                                                 cATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGa   >  1:2675511/1‑100 (MQ=255)
                                                                                                  aTTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGaa  <  1:1048059/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
ATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAA  >  NC_000913/1597173‑1597370

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