Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,599,943 G→A A144T (GCC→ACC)  b1513 → putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,599,9430GA98.8% 543.1 / ‑6.0 163A144T (GCC→ACC) b1513putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (82/79);  total (82/81)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 171 total aligned reads displayed.

AAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCGCCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGAACCTACCGCCTCGCTTACCCCT  >  NC_000913/1599853‑1600032
                                                                                          |                                                                                         
  aCAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGa                                                                                       >  1:7071956/1‑93 (MQ=255)
    aGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATc                                                                                     <  1:7292206/93‑1 (MQ=255)
    aGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATc                                                                                     >  2:1153303/1‑93 (MQ=255)
     gCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCg                                                                                    >  2:7808220/1‑93 (MQ=255)
      ctctCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGc                                                                                   >  1:1091121/1‑93 (MQ=255)
      ctctCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGc                                                                                   <  1:95302/93‑1 (MQ=255)
      ctctCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGAGGTCACCGATCGc                                                                                   >  2:1225478/1‑93 (MQ=255)
         tcCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCaa                                                                                >  2:6053551/1‑93 (MQ=255)
           cATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAAt                                                                              <  1:5368656/93‑1 (MQ=255)
           cATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAAt                                                                              <  1:232763/93‑1 (MQ=255)
           cATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAAt                                                                              <  1:6837816/93‑1 (MQ=255)
           cATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAAt                                                                              <  2:5081593/93‑1 (MQ=255)
            aTGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAAtg                                                                             >  1:3747300/1‑93 (MQ=255)
             tGCNGAAAATGAAGAACTTGCTGGNGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAAtgg                                                                            <  2:5649781/93‑1 (MQ=255)
              gCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAAtggt                                                                           <  1:5616914/93‑1 (MQ=255)
               cAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAAtggtg                                                                          <  2:5306404/93‑1 (MQ=255)
               cAGAAAATGAAGAACTTCCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAAtggtg                                                                          <  2:1880365/93‑1 (MQ=255)
                 gAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGa                                                                        >  1:668023/1‑93 (MQ=255)
                 gAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGa                                                                        >  1:4209061/1‑93 (MQ=255)
                  aaaaTGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGGTCGCTGGATGTGACCGATCGCCAAATGGGGGaa                                                                       >  2:7109655/1‑93 (MQ=255)
                  aaaaTGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGaa                                                                       >  1:4706052/1‑93 (MQ=255)
                    aaTGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAAt                                                                     <  1:4205324/93‑1 (MQ=255)
                     aTGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAAtc                                                                    >  1:7587137/1‑93 (MQ=255)
                      tGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGGTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAAtcc                                                                   >  2:5906914/1‑93 (MQ=255)
                       gaagaaCTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAAtcct                                                                  <  1:2916207/93‑1 (MQ=255)
                        aagaaCTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAAtcctc                                                                 <  1:1968766/93‑1 (MQ=255)
                        aagaaCTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAAtcctc                                                                 <  2:6182971/93‑1 (MQ=255)
                         agaaCTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAAtcctcc                                                                <  2:1029694/93‑1 (MQ=255)
                         agaaCTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAAtcctcc                                                                <  2:2535548/93‑1 (MQ=255)
                          gaaCTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCGCCGATCGCCAAATGGTGGAAATCCTCcg                                                               <  2:7894240/93‑1 (MQ=255)
                           aaCTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCcgc                                                              <  1:1204246/93‑1 (MQ=255)
                           aaCTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCcgc                                                              >  1:6208500/1‑93 (MQ=255)
                            aCTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCcgcg                                                             >  1:5043440/1‑93 (MQ=255)
                                gCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCt                                                         <  2:3383859/93‑1 (MQ=255)
                                gCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCt                                                         <  1:4466631/93‑1 (MQ=255)
                                  tGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGa                                                       <  2:2717277/93‑1 (MQ=255)
                                  tGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGa                                                       <  1:6876890/93‑1 (MQ=255)
                                    gcggcgCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATg                                                     <  2:3697921/93‑1 (MQ=255)
                                      ggcgCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATgcg                                                   >  1:6974339/1‑93 (MQ=255)
                                       gcgcTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATgcgc                                                  >  1:7703092/1‑93 (MQ=255)
                                        cgcTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGcgcg                                                 >  1:3226811/1‑93 (MQ=255)
                                         gcTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGa                                                <  1:3818701/93‑1 (MQ=255)
                                          cTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGAc                                               <  1:603856/93‑1 (MQ=255)
                                          cTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGAc                                               <  2:3571951/93‑1 (MQ=255)
                                          cTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGAc                                               >  2:7653529/1‑93 (MQ=255)
                                           tGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACt                                              >  1:6446117/1‑93 (MQ=255)
                                           tGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACt                                              >  2:3857618/1‑93 (MQ=255)
                                            gggCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTc                                             >  2:1591435/1‑93 (MQ=255)
                                            gggCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTc                                             <  2:1780989/93‑1 (MQ=255)
                                            gggCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTc                                             <  1:397250/93‑1 (MQ=255)
                                               cTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTcgcg                                          >  2:4639469/1‑93 (MQ=255)
                                                tGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCgg                                         >  2:3207060/1‑93 (MQ=255)
                                                tGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCgg                                         >  1:162413/1‑93 (MQ=255)
                                                 gCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGa                                        >  1:1472228/1‑93 (MQ=255)
                                                  ccAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGAt                                       <  2:2070219/93‑1 (MQ=255)
                                                    aGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTc                                     <  1:2640709/93‑1 (MQ=255)
                                                    aGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTc                                     >  2:2754191/1‑93 (MQ=255)
                                                    aGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCAGCGGGCTGATGCGCGACTCGCGGATTc                                     >  1:4025640/1‑93 (MQ=255)
                                                     gTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCt                                    >  1:5895471/1‑93 (MQ=255)
                                                     gTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCt                                    <  1:191007/93‑1 (MQ=255)
                                                     gTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCt                                    >  1:4480010/1‑93 (MQ=255)
                                                     gTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCt                                    >  1:398520/1‑93 (MQ=255)
                                                       ttGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGa                                  <  1:2996317/93‑1 (MQ=255)
                                                       ttGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGa                                  >  2:6073932/1‑93 (MQ=255)
                                                       ttGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGa                                  <  1:6253069/93‑1 (MQ=255)
                                                       ttGATCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGa                                  <  1:7251813/93‑1 (MQ=255)
                                                          aTCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATcc                               <  1:3023224/93‑1 (MQ=255)
                                                           tCTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCt                              >  2:2616111/1‑93 (MQ=255)
                                                            cTGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTc                             <  1:6406174/93‑1 (MQ=255)
                                                             tGCATAGTCTGGCAGGATCGCTGGCTGTCACCGATCGCCACATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCg                            >  1:3526271/1‑93 (MQ=255)
                                                             tGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCg                            <  1:5254885/93‑1 (MQ=255)
                                                             tGCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCg                            <  1:5254886/93‑1 (MQ=255)
                                                              gCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGa                           <  2:3867667/93‑1 (MQ=255)
                                                              gCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGa                           >  1:4070488/1‑93 (MQ=255)
                                                              gCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGa                           >  1:6123603/1‑93 (MQ=255)
                                                              gCATAGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGa                           <  1:252133/93‑1 (MQ=255)
                                                                  aGTCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGaa                       >  1:6991830/1‑93 (MQ=255)
                                                                    tCTGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGAAcc                     <  2:1147637/93‑1 (MQ=255)
                                                                      tGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGAACCTa                   >  2:6855554/1‑93 (MQ=255)
                                                                      tGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGAACCTa                   >  2:5400744/1‑93 (MQ=255)
                                                                      tGGCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGAACCTa                   <  2:826459/93‑1 (MQ=255)
                                                                       ggCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGAACCTAc                  <  2:283752/93‑1 (MQ=255)
                                                                       ggCAGGATCGCTGGAGGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGAACCtgc                  <  2:2017198/93‑3 (MQ=255)
                                                                        gCAGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGAACCTAcc                 <  1:4085028/93‑1 (MQ=255)
                                                                          aGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGAACCTACCGc               >  1:3141937/1‑93 (MQ=255)
                                                                          aGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGAACCTACCGc               <  2:2552220/93‑1 (MQ=255)
                                                                          aGGATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGAACCTACCGc               <  2:6831024/93‑1 (MQ=255)
                                                                           ggATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGAACCTACCGcc              >  1:6463098/1‑93 (MQ=255)
                                                                           ggATCGCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGAACCTACCGcc              >  2:4317230/1‑93 (MQ=255)
                                                                            gaccGCTGGATGTCGCGGATCGCCAAATGGTGGAAAGCCTCCGCGGGCTGATGCGCGGCTCGCGGATTCTGATCCTCGATGAACCTACCGCCt             <  2:6728293/90‑1 (MQ=255)
                                                                               cGCTGTATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGAACCTACCGCCTCGc          <  2:4455969/93‑1 (MQ=255)
                                                                                gCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGAACCTACCGCCTCGCt         >  1:2247771/1‑93 (MQ=255)
                                                                                gCTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGAACCTACCGCCTCGCt         >  2:6279301/1‑93 (MQ=255)
                                                                                 cTGGATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGAACCTACCGCCTCGCtt        >  1:4909475/1‑93 (MQ=255)
                                                                                    gATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGTTGATGCGCGCCTCGCGGATTCTAATCCTCGATGAACCTACCGCCTCCCTTccc     >  2:4161681/1‑93 (MQ=255)
                                                                                    gATGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGAACCTACCGCCTCGCTTAcc     <  1:3267572/93‑1 (MQ=255)
                                                                                      tGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGAACCTACCGCCTCGCTTAcccc   >  2:2494519/1‑93 (MQ=255)
                                                                                      tGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGAACCTACCGCCTCGCTTAcccc   <  2:4844708/93‑1 (MQ=255)
                                                                                      tGTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTAGATGAACCTACCGCCTCGCTTAcccc   <  2:4475717/93‑1 (MQ=255)
                                                                                       gTCACCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGAACCTACCGCCTCGCTTACCCCt  >  2:1534502/1‑93 (MQ=255)
                                                                                          |                                                                                         
AAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGGCAGGATCGCTGGATGTCGCCGATCGCCAAATGGTGGAAATCCTCCGCGGGCTGATGCGCGACTCGCGGATTCTGATCCTCGATGAACCTACCGCCTCGCTTACCCCT  >  NC_000913/1599853‑1600032

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: