Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,597,272 T→A 15.8% 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,2720TA15.8% 416.2 / 67.6 171pseudogene (716/3861 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base T (72/71);  new base A (11/16);  total (83/88)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-02

Only 100 of 173 total aligned reads displayed.

TATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAA  >  NC_000913/1597174‑1597370
                                                                                                  |                                                                                                  
tatTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAtt                                                                                                   >  2:3978393/1‑100 (MQ=255)
 atTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTg                                                                                                  <  2:1495761/100‑1 (MQ=255)
 atTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTg                                                                                                  <  2:3347514/100‑1 (MQ=255)
 atTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTg                                                                                                  <  2:2999259/100‑1 (MQ=255)
  ttGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGc                                                                                                 <  2:1372549/100‑1 (MQ=255)
    gACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCa                                                                                               <  2:3299786/100‑1 (MQ=255)
      cTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCAcc                                                                                             <  1:871587/100‑1 (MQ=255)
       tgtgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCg                                                                                            <  2:4314578/100‑1 (MQ=255)
         tgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGAt                                                                                          <  1:3137734/100‑1 (MQ=255)
         tgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCggt                                                                                          <  2:2671782/100‑3 (MQ=255)
         tgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTCCCGATACCGGTGGCAATGCCACCggt                                                                                          <  2:2671233/100‑3 (MQ=255)
          gCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTGCGATATCAGTGGCATTGCCACCGATa                                                                                         >  1:474286/1‑100 (MQ=255)
          gCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCCATGCCACCGGTa                                                                                         <  2:2887553/100‑1 (MQ=255)
          gCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCGGTa                                                                                         <  1:976064/100‑1 (MQ=255)
          gCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCGGTa                                                                                         <  2:2140071/100‑1 (MQ=255)
            cATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCGGTatc                                                                                       <  2:322471/100‑3 (MQ=255)
            cATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCGGTatc                                                                                       <  1:144011/100‑3 (MQ=255)
             aTTACTGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGCTCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATc                                                                                      <  1:1021073/100‑1 (MQ=255)
             aTTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCggtatcg                                                                                      <  2:1301092/100‑7 (MQ=255)
             aTTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCggtatcg                                                                                      <  1:2447610/100‑7 (MQ=255)
              ttACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCg                                                                                     <  1:2187866/100‑1 (MQ=255)
              ttACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCggtatcg                                                                                      <  1:3294165/99‑7 (MQ=255)
               naCCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGt                                                                                    >  1:2725438/2‑100 (MQ=255)
                aCCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCATGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTa                                                                                   <  1:2353685/100‑1 (MQ=255)
                 ccGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAt                                                                                  <  1:574687/100‑1 (MQ=255)
                  cGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAt                                                                                  >  2:3694605/1‑99 (MQ=255)
                   gAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAttt                                                                                <  1:3347958/100‑1 (MQ=255)
                   gAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAttt                                                                                >  2:1478562/1‑100 (MQ=255)
                   gAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAttt                                                                                >  1:771300/1‑100 (MQ=255)
                   gAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAttt                                                                                >  2:4108975/1‑100 (MQ=255)
                    aaaTATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAtttt                                                                               >  1:1104515/1‑100 (MQ=255)
                        atCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGc                                                                           >  2:1262176/1‑100 (MQ=255)
                        atCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGc                                                                           >  1:1550114/1‑100 (MQ=255)
                        atCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCATGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGc                                                                           >  2:3766231/1‑100 (MQ=255)
                         tCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGcc                                                                          <  2:1921120/100‑1 (MQ=255)
                          ccATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGcc                                                                          >  1:1804126/1‑99 (MQ=255)
                          ccATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCg                                                                         >  2:1923294/1‑100 (MQ=255)
                          ccATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCg                                                                         >  1:3115157/1‑100 (MQ=255)
                           cATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGt                                                                        >  1:2490490/1‑100 (MQ=255)
                           cATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGt                                                                        >  1:3823431/1‑100 (MQ=255)
                           cATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAGCGTATTTTCAGCCGt                                                                        >  2:4339370/1‑100 (MQ=255)
                             ttATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGACGTaa                                                                      >  2:2385868/1‑100 (MQ=255)
                              tataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTTAGCCGTaac                                                                     <  2:708215/100‑1 (MQ=255)
                              tataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaac                                                                     >  2:541403/1‑100 (MQ=255)
                              tataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaac                                                                     >  2:2034473/1‑100 (MQ=255)
                                 agagccgacGTAACGTGATGTCAGAGTCTTGATCCTGCACATCCAGGGTTCCGATACCGGTTGCAGTGCCACCGACAATCGTATTTCCAGCCGTa                                                                       <  1:3604246/87‑1 (MQ=255)
                                   attCTGCGTACCGAGATGACAGAGTCTTGATCCTGCACATTCGGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaccaacc                                                                >  2:2979409/3‑100 (MQ=255)
                                        gCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCAGCGATAATCGTATTTTCAGCCGTAACCAAccccccc                                                           >  2:4292083/1‑100 (MQ=255)
                                            aCCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTcgc                                                       <  2:2442541/100‑1 (MQ=255)
                                             ccGTGATGACATAGTCTTGATCCTGCACATTTAGTTTTCTGATACCTGTCGCATTGCCACCGATAATCGTATTTTCGGCTGTAAGCAACCCCCCcgcgcg                                                      <  1:293372/100‑1 (MQ=255)
                                              cGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTcgcg                                                      <  1:2804253/99‑1 (MQ=255)
                                               gTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGaa                                                    >  1:3987782/1‑100 (MQ=255)
                                               gTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGaa                                                    <  2:2904376/100‑1 (MQ=255)
                                               gTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGaa                                                    <  2:2189454/100‑1 (MQ=255)
                                                tgatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCTATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAAt                                                   <  1:3619219/100‑1 (MQ=255)
                                                 gatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATa                                                  >  1:639899/1‑100 (MQ=255)
                                                 gatgaCAGAGTCTTGATCCTGCACAGTCAGGGTTCCGATACCGGTTGCTTTGCCACCGATAATCGTATTTTCAGCCGTTACGAACTCCCCCTCTCGaag                                                   >  1:2850395/1‑98 (MQ=255)
                                                  atgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAg                                                 <  2:1550119/100‑1 (MQ=255)
                                                   tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGt                                                >  2:3092454/1‑100 (MQ=255)
                                                      cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGc                                             >  1:1763892/1‑100 (MQ=255)
                                                      cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGc                                             >  1:2549791/1‑100 (MQ=255)
                                                      cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCAACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCg                                              >  2:907702/1‑99 (MQ=255)
                                                           tCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTccc                                        >  2:2103717/1‑100 (MQ=255)
                                                            cTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTcccc                                       >  2:1398738/1‑100 (MQ=255)
                                                              tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAg                                     >  2:2640664/1‑100 (MQ=255)
                                                              tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAg                                     >  1:4069928/1‑100 (MQ=255)
                                                                  ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTcc                                 >  1:3202601/1‑100 (MQ=255)
                                                                  ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTcc                                 >  1:4266125/1‑100 (MQ=255)
                                                                    tGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCtt                               <  1:1268709/100‑1 (MQ=255)
                                                                    tGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCtt                               <  1:3725097/100‑1 (MQ=255)
                                                                     gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                              <  1:393451/100‑1 (MQ=255)
                                                                      cacaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCtt                               <  1:54716/98‑1 (MQ=255)
                                                                       acaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGCTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAt                            <  2:3095461/100‑1 (MQ=255)
                                                                        caTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                           <  2:1122548/100‑1 (MQ=255)
                                                                         aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                           >  1:2681055/1‑99 (MQ=255)
                                                                            cAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGACTAGTCGCCTCCCCAGTTCCTTGATTACCa                       >  2:3503492/1‑100 (MQ=255)
                                                                             aGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATCACCaa                      <  1:2942598/100‑1 (MQ=255)
                                                                               ggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAtt                    <  1:1010463/100‑1 (MQ=255)
                                                                                gTTCCGATACCGGTGCCATTGCCACCGAGAACCGTATTTTCAGCCGGAACTAACCGACCCTCGCGAATAGCCGCCGACGCAGTTCCTGGATTACCAAtt                    <  1:3779074/99‑1 (MQ=255)
                                                                                  tCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGa                 <  2:2855617/100‑1 (MQ=255)
                                                                                  acgattatcGGTGGCAATGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTg                  >  1:2283990/9‑99 (MQ=38)
                                                                                   ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaa                <  2:2490502/100‑1 (MQ=255)
                                                                                   ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaa                <  1:2281167/100‑1 (MQ=255)
                                                                                    gattatcGGTGGCAATGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa               >  1:1243189/7‑100 (MQ=39)
                                                                                     gATACCGGTGGCATTGCCACCGATAATCGTATTGTCAGCCGGAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCCGTTCCTTGATTACCACTTTGaaa               >  1:1681340/1‑99 (MQ=255)
                                                                                       tatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTc            >  2:2050045/4‑100 (MQ=255)
                                                                                       tatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTc            >  2:3404591/4‑100 (MQ=255)
                                                                                       tatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTc            >  2:231735/4‑100 (MQ=255)
                                                                                       tatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTc            >  2:2048635/4‑100 (MQ=255)
                                                                                         ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCaa          <  1:3501352/100‑1 (MQ=255)
                                                                                         ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCaa          >  2:3943830/1‑100 (MQ=255)
                                                                                         ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCaa          >  2:2355160/1‑100 (MQ=255)
                                                                                           ggTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt        <  1:3194322/100‑1 (MQ=255)
                                                                                           ggTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt        >  2:245851/1‑100 (MQ=255)
                                                                                           ggTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAt         <  1:3023955/99‑1 (MQ=255)
                                                                                            gTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTg       <  1:3696857/100‑1 (MQ=255)
                                                                                              ggCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAg     >  2:472514/1‑100 (MQ=255)
                                                                                              ggCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAg     <  2:1613872/100‑1 (MQ=255)
                                                                                              ggCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCCCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGa      >  1:2595642/1‑99 (MQ=255)
                                                                                                cATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGATCCTTGATTACAAATTTAAAATTCAAtt        >  1:4046541/1‑95 (MQ=255)
                                                                                                  |                                                                                                  
TATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAA  >  NC_000913/1597174‑1597370

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: ‑