Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,597,272 T→A 17.3% 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,2720TA17.3% 249.7 / 43.2 104pseudogene (716/3861 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base T (51/35);  new base A (11/7);  total (62/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

Only 100 of 107 total aligned reads displayed.

ATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAA  >  NC_000913/1597173‑1597370
                                                                                                   |                                                                                                   
atatTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAt                                                                                                     <  1:1218274/100‑1 (MQ=255)
atatTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAt                                                                                                     <  1:1703397/100‑1 (MQ=255)
  atTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTg                                                                                                   >  2:1094275/1‑100 (MQ=255)
    tGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTgtc                                                                                                 <  1:1281621/100‑3 (MQ=255)
     gACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCTTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCa                                                                                                >  2:510475/1‑100 (MQ=255)
     gACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCa                                                                                                >  1:343744/1‑100 (MQ=255)
      aCTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCAc                                                                                               <  1:728235/100‑1 (MQ=255)
         gtgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGa                                                                                            <  1:1031082/100‑1 (MQ=255)
         gtgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGa                                                                                            >  1:2493158/1‑100 (MQ=255)
         gtgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCgg                                                                                            <  1:636632/100‑2 (MQ=255)
          tgCCATTACCTAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCggt                                                                                           <  1:628618/100‑3 (MQ=255)
          tgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCggt                                                                                           <  1:346422/100‑3 (MQ=255)
          tgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCggt                                                                                           <  1:199027/100‑3 (MQ=255)
          tgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCggt                                                                                           <  1:1763610/100‑3 (MQ=255)
            ccATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATaa                                                                                         >  1:1917041/1‑100 (MQ=255)
            ccATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATaa                                                                                         <  1:1683018/100‑1 (MQ=255)
             cATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCGGTatc                                                                                        <  1:2200819/100‑3 (MQ=255)
              aTTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCggtatcg                                                                                       <  2:410964/100‑7 (MQ=255)
               ttACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCg                                                                                      <  2:1990517/100‑1 (MQ=255)
                  ccGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAt                                                                                   >  1:277815/1‑100 (MQ=255)
                   cGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAtt                                                                                  <  1:1236633/100‑1 (MQ=255)
                    gAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAttt                                                                                 >  1:252764/1‑100 (MQ=255)
                      aaTATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTc                                                                               >  1:1598365/1‑100 (MQ=255)
                      aaTATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTc                                                                               <  1:825381/100‑1 (MQ=255)
                        tatCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATGGTATTTTCAg                                                                             >  2:2442301/1‑100 (MQ=255)
                        tatCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAg                                                                             <  1:1027253/100‑1 (MQ=255)
                         atCCATTATCGAGTCTGCGTACCATGCTGACAGAGCCTTCATCCTGCCCATTCAGGGTTCCGATACCGGTGGCATGGCCACCGATAATCGTCtttt                                                                                >  2:365470/1‑96 (MQ=255)
                         atCCATTATACGGTCTGCGTACCGTTATGACAGATTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGc                                                                            <  2:1225411/100‑1 (MQ=255)
                           ccATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCg                                                                          >  1:1843438/1‑100 (MQ=255)
                           ccATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCg                                                                          >  1:1125298/1‑100 (MQ=255)
                           ccATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCg                                                                          <  1:1293130/100‑1 (MQ=255)
                               tataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaac                                                                      <  1:1568422/100‑1 (MQ=255)
                                   gagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaccaac                                                                 >  2:800325/1‑100 (MQ=255)
                                    agTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaccaacc                                                                >  2:427554/1‑100 (MQ=255)
                                         cgatctcGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTGCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAccccccc                                                           <  2:1793550/94‑1 (MQ=255)
                                          cGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCt                                                          >  1:1680533/1‑100 (MQ=255)
                                            tACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTcg                                                        >  1:1067257/1‑100 (MQ=255)
                                              ccGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTcgcg                                                      >  2:1762462/1‑100 (MQ=255)
                                              ccGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTcgcg                                                      <  2:179218/100‑1 (MQ=255)
                                               cGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAcccccccagatcg                                                     >  1:1793564/1‑94 (MQ=255)
                                               cGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGa                                                     >  2:1987019/1‑100 (MQ=255)
                                                gcgatgaCAGAGTCTTGATCCTGCTCATTCAGGGGTCCGATACCGGTGGCATTGCCACCGATAATCGTCTTTTCAGCCGTAACCAACCCCCCCTCGCGa                                                     >  2:434590/3‑99 (MQ=255)
                                                gTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGaa                                                    <  2:1446180/100‑1 (MQ=255)
                                                  gatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATa                                                  >  2:49333/1‑100 (MQ=255)
                                                  gatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATa                                                  <  1:1479231/100‑1 (MQ=255)
                                                   atgaCAGAGTCTTGATCCTGCACATTCAGGCTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAg                                                 <  2:327473/100‑1 (MQ=255)
                                                    tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAg                                                 >  2:761280/1‑99 (MQ=255)
                                                    tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGt                                                >  2:1548155/1‑100 (MQ=255)
                                                    tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGt                                                >  1:703941/1‑100 (MQ=255)
                                                       cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGc                                             <  1:1487325/100‑1 (MQ=255)
                                                               tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAg                                     >  1:2368212/1‑100 (MQ=255)
                                                               tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAg                                     >  1:628977/1‑100 (MQ=255)
                                                               tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAg                                     >  1:177725/1‑100 (MQ=255)
                                                                gATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGt                                    <  1:1565204/100‑1 (MQ=255)
                                                                 aTCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGtt                                   >  1:2164254/1‑100 (MQ=255)
                                                                   ccTGCACATTCAGGGTTCCGATATCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTcc                                 <  1:48480/100‑1 (MQ=255)
                                                                   ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTcc                                 <  1:2514353/100‑1 (MQ=255)
                                                                   ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTcc                                 <  2:1257481/100‑1 (MQ=255)
                                                                   ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTcc                                 >  1:1747358/1‑100 (MQ=255)
                                                                   ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTcc                                 >  1:1996205/1‑100 (MQ=255)
                                                                     tGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCtt                               <  2:2185715/100‑1 (MQ=255)
                                                                      gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                              >  1:2485162/1‑100 (MQ=255)
                                                                       cacaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGa                             >  2:956708/1‑100 (MQ=255)
                                                                        acaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAt                            <  1:200875/100‑1 (MQ=255)
                                                                         caTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                           <  2:705167/100‑1 (MQ=255)
                                                                          aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTa                          >  2:518032/1‑100 (MQ=255)
                                                                          aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTa                          >  1:123843/1‑100 (MQ=255)
                                                                            tCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAcc                        <  1:298051/100‑1 (MQ=255)
                                                                               gggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAt                     >  1:758743/1‑100 (MQ=255)
                                                                                ggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAtt                    >  1:541549/1‑100 (MQ=255)
                                                                                ggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAtt                    <  2:1747346/100‑1 (MQ=255)
                                                                                   tCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGa                 <  2:2049190/100‑1 (MQ=255)
                                                                                   acgattatcGGTGGCAATGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCCGTTCCTTGATTACCAATTTGa                 >  2:733329/9‑100 (MQ=38)
                                                                                    cgattatcGGTGGCAATGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaa                >  2:357349/8‑100 (MQ=38)
                                                                                    ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaa                <  1:1924642/100‑1 (MQ=255)
                                                                                    ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaa                >  2:2286194/1‑100 (MQ=255)
                                                                                    ccGACACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAATCAGTCCACCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAACTTGaa                >  2:1155472/1‑100 (MQ=255)
                                                                                     gattatcGGTGGCAATGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa               >  2:334909/7‑100 (MQ=39)
                                                                                     cGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa               <  2:334713/100‑1 (MQ=255)
                                                                                     cGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa               >  2:675373/1‑100 (MQ=255)
                                                                                      gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCGGTTCCTTGATTACCAATTTGAAAt              >  1:894321/1‑100 (MQ=255)
                                                                                      gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt              <  2:737857/100‑1 (MQ=255)
                                                                                      attatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAGTAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt              >  2:1390717/6‑100 (MQ=255)
                                                                                      attatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt              >  1:1412033/6‑100 (MQ=255)
                                                                                       ttatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAtt             >  2:262162/5‑100 (MQ=255)
                                                                                         atcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa           >  2:711071/3‑100 (MQ=255)
                                                                                         atcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa           >  2:1604808/3‑100 (MQ=255)
                                                                                         atcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTCAAATTc            >  2:1094280/3‑99 (MQ=255)
                                                                                         atcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCAAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa           >  2:1636594/3‑100 (MQ=255)
                                                                                         aCCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTc            >  2:2202503/1‑99 (MQ=255)
                                                                                          ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa           <  1:342115/100‑1 (MQ=255)
                                                                                          ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCaa          >  1:686197/1‑100 (MQ=255)
                                                                                          ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCaa          >  1:837494/1‑100 (MQ=255)
                                                                                            ggTGGCATTGCAACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTGCCTTGATTACCAATTTGAAATTCAAtt        >  2:26155/1‑100 (MQ=255)
                                                                                             gTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTg       <  2:1000809/100‑1 (MQ=255)
                                                                                                gCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg    <  1:49334/100‑1 (MQ=255)
                                                                                                 cATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCACTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGa   >  1:686932/1‑100 (MQ=255)
                                                                                                  aTTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGaa  >  1:183387/1‑100 (MQ=255)
                                                                                                  aTTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGaa  >  1:656910/1‑100 (MQ=255)
                                                                                                   ttGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGaa  <  1:913423/99‑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 < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: