Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,966,598 G→A 100% G61E (GGA→GAA)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,966,5980GA99.2% 422.6 / ‑4.9 118G61E (GGA→GAA) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (58/59);  total (58/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 122 total aligned reads displayed.

ATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCA  >  NC_000913/3966501‑3966693
                                                                                                 |                                                                                               
aTGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggatgaat                                                                                               <  2:3054102/100‑1 (MQ=255)
  gCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAAttt                                                                                             >  2:3266307/1‑100 (MQ=255)
   cGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTg                                                                                            <  2:2866675/100‑1 (MQ=255)
   cGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTg                                                                                            >  2:2961868/1‑100 (MQ=255)
    gTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTgg                                                                                           <  2:302004/100‑1 (MQ=255)
    gTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTgg                                                                                           <  2:1649821/100‑1 (MQ=255)
     tAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGt                                                                                          <  2:1067207/100‑1 (MQ=255)
       aGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGttt                                                                                        >  1:3192712/1‑100 (MQ=255)
         cAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTtcc                                                                                      >  1:1693612/1‑100 (MQ=255)
          aGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTtcct                                                                                     <  1:2112185/100‑1 (MQ=255)
           ggACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTtcctc                                                                                    >  1:3049812/1‑100 (MQ=255)
            gACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTtcctcc                                                                                   <  2:118165/100‑1 (MQ=255)
              cATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGt                                                                                 >  2:1415037/1‑100 (MQ=255)
              cATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGt                                                                                 <  2:1093331/100‑1 (MQ=255)
                ttattTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTc                                                                               <  1:962669/100‑1 (MQ=255)
                ttattTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTc                                                                               <  2:1019743/100‑1 (MQ=255)
                 tattTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTcc                                                                              >  2:3084681/1‑100 (MQ=255)
                     tttGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAg                                                                          <  1:1992642/100‑1 (MQ=255)
                       tGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGAc                                                                        >  1:799997/1‑100 (MQ=255)
                       tGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGAc                                                                        <  2:2833530/100‑1 (MQ=255)
                        gCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACa                                                                       >  1:2894300/1‑100 (MQ=255)
                        gCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACa                                                                       >  2:2697817/1‑100 (MQ=255)
                         ccATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATGTGGTTTCCTCCGTTCCGCAGACAg                                                                      >  1:1503695/1‑100 (MQ=255)
                          cATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGc                                                                     >  2:714558/1‑100 (MQ=255)
                           aTCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCt                                                                    <  2:316038/100‑1 (MQ=255)
                             ccTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTcc                                                                  >  2:2150605/1‑100 (MQ=255)
                             ccTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTcc                                                                  <  2:511981/100‑1 (MQ=255)
                             ccTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTcc                                                                  >  2:2945833/1‑100 (MQ=255)
                             ccTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTcc                                                                  <  1:1647428/100‑1 (MQ=255)
                               tGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTa                                                                >  1:1579722/1‑100 (MQ=255)
                                gAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTAc                                                               >  2:1939716/1‑100 (MQ=255)
                                 aaGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTAcc                                                              <  2:782732/100‑1 (MQ=255)
                                  agcagcACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCt                                                             <  2:2568095/100‑1 (MQ=255)
                                  agcagcACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCt                                                             >  2:1375070/1‑100 (MQ=255)
                                     agcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGc                                                          >  2:1698955/1‑100 (MQ=255)
                                      gcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGcc                                                         <  1:981486/100‑1 (MQ=255)
                                      gcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGcc                                                         <  2:3224674/100‑1 (MQ=255)
                                        aCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCgg                                                       <  1:327234/100‑1 (MQ=255)
                                        aCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCgg                                                       >  1:1513885/1‑100 (MQ=255)
                                         cGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGt                                                      <  2:1889970/100‑1 (MQ=255)
                                           cAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTcc                                                    <  1:3119359/100‑1 (MQ=255)
                                           cAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTcc                                                    >  2:3067645/1‑100 (MQ=255)
                                            aaaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCt                                                   <  2:1492719/100‑1 (MQ=255)
                                            aaaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCt                                                   >  1:910964/1‑100 (MQ=255)
                                                agTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgatg                                               <  1:2360551/100‑1 (MQ=255)
                                                agTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgatg                                               <  2:344326/100‑1 (MQ=255)
                                                agTGGCGAAGATATCGTTGGTGATGGCGTACTGGAGATATTCCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgat                                                <  1:1157452/99‑1 (MQ=255)
                                                 gTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgatga                                              >  2:2158558/1‑100 (MQ=255)
                                                 gTGGCGAAGATATCTTTGGTGATGGCGTACGGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCTGTCCtgatg                                               >  1:2192682/1‑99 (MQ=255)
                                                  tGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCGTGATGAc                                             <  2:1762025/100‑1 (MQ=255)
                                                  tGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGAc                                             <  2:89545/100‑1 (MQ=255)
                                                   ggCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACa                                            >  1:2077541/1‑100 (MQ=255)
                                                   ggCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACa                                            >  2:2623336/1‑100 (MQ=255)
                                                   ggCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACa                                            <  2:580580/100‑1 (MQ=255)
                                                   ggCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACa                                            >  1:3180259/1‑100 (MQ=255)
                                                   ggCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACa                                            >  2:1081956/1‑100 (MQ=255)
                                                   ggCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACa                                            <  1:2753666/100‑1 (MQ=255)
                                                    gCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACAt                                           <  2:2114998/100‑1 (MQ=255)
                                                     cGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACAtn                                          <  1:3257978/100‑2 (MQ=255)
                                                     cGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATc                                          <  2:1379194/100‑1 (MQ=255)
                                                      gAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCt                                         >  1:923302/1‑100 (MQ=255)
                                                            atCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGttt                                   >  2:2228297/1‑100 (MQ=255)
                                                              cTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTcc                                 >  2:1579109/1‑100 (MQ=255)
                                                              cTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTcc                                 >  2:1569843/1‑100 (MQ=255)
                                                              cTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTcc                                 >  2:1569426/1‑100 (MQ=255)
                                                               tttGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTccc                                <  2:3049838/100‑1 (MQ=255)
                                                                ttGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTcccc                               <  1:1415032/100‑1 (MQ=255)
                                                                  ggTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTa                             <  2:799999/100‑1 (MQ=255)
                                                                  ggTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTa                             >  2:2950960/1‑100 (MQ=255)
                                                                   gTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAg                            >  2:303636/1‑100 (MQ=255)
                                                                    tGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGc                           <  2:1205806/100‑1 (MQ=255)
                                                                    tGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGc                           <  1:2815258/100‑1 (MQ=255)
                                                                    tGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTACCCTAGc                           >  1:2500336/1‑100 (MQ=255)
                                                                     gATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGcc                          <  1:3266281/100‑1 (MQ=255)
                                                                       tGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCaa                        >  2:670574/1‑100 (MQ=255)
                                                                         gCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAAt                      <  2:810250/100‑1 (MQ=255)
                                                                          cGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATc                     >  1:2618966/1‑100 (MQ=255)
                                                                            tACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccg                   <  2:1565398/100‑1 (MQ=255)
                                                                             aCTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccgc                  >  1:647749/1‑100 (MQ=255)
                                                                             aCTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccgc                  <  2:1198894/100‑1 (MQ=255)
                                                                              cTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccgcc                 >  1:828531/1‑100 (MQ=255)
                                                                              cTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccgcc                 <  1:2925547/100‑1 (MQ=255)
                                                                              cTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGGTTCCCCTAGCCAAATccgcc                 <  1:2053849/100‑1 (MQ=255)
                                                                                   gaTATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTc            >  2:2457301/1‑100 (MQ=255)
                                                                                   gaTATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTc            >  1:1402453/1‑100 (MQ=255)
                                                                                     tatTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCaa          >  1:3322400/1‑100 (MQ=255)
                                                                                     tatTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCaa          >  1:1385721/1‑100 (MQ=255)
                                                                                      atTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCTCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTc            <  1:134203/97‑1 (MQ=255)
                                                                                      atTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAAc         >  2:2553518/1‑100 (MQ=255)
                                                                                      atTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAAc         <  2:1402458/100‑1 (MQ=255)
                                                                                       cgatctgatgaatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAAcc        <  1:2242730/94‑1 (MQ=255)
                                                                                        tGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCt       >  1:3052363/1‑100 (MQ=255)
                                                                                        tGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCt       >  2:68039/1‑100 (MQ=255)
                                                                                        tGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCt       >  2:1785292/1‑100 (MQ=255)
                                                                                          cAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTcc     >  2:2718315/1‑100 (MQ=255)
                                                                                          cAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTcc     <  2:191583/100‑1 (MQ=255)
                                                                                          cAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTcc     >  2:1744720/1‑100 (MQ=255)
                                                                                           aGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCg    >  2:2514234/1‑100 (MQ=255)
                                                                                            ggatgaatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGc   >  2:1667903/1‑100 (MQ=255)
                                                                                             gatgaatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAAccagatcg  >  2:2242747/1‑94 (MQ=255)
                                                                                                 |                                                                                               
ATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCA  >  NC_000913/3966501‑3966693

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