Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,966,875 T→G S153S (TCT→TCG rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,966,8750TG97.7% 556.6 / ‑1.4 171S153S (TCT→TCGrhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base G (77/89);  total (80/91)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 171 total aligned reads displayed.

TCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACC  >  NC_000913/3966778‑3966971
                                                                                                 |                                                                                                
tCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGAc                                                                                                <  1:1061750/100‑1 (MQ=255)
tCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGAc                                                                                                <  2:3314528/100‑1 (MQ=255)
tCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGAc                                                                                                <  2:1978367/100‑1 (MQ=255)
 cGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACt                                                                                               >  1:81176/1‑100 (MQ=255)
 cGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACt                                                                                               >  1:3637462/1‑100 (MQ=255)
 cGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACt                                                                                               <  1:1946138/100‑1 (MQ=255)
 cGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACt                                                                                               >  1:1584724/1‑100 (MQ=255)
  gACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTg                                                                                              >  2:2667186/1‑100 (MQ=255)
    cAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGaa                                                                                            >  2:2067354/1‑100 (MQ=255)
     aaaCCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAg                                                                                           <  2:2671037/100‑1 (MQ=255)
      aaCCTGAAAACGCCCGCAACAAAATCCTCTTTGCGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTACCGGTTCGACTAAAGa                                                                                          >  2:3169190/1‑100 (MQ=255)
      aaCCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGa                                                                                          >  1:2759531/1‑100 (MQ=255)
        ccTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGAtt                                                                                        <  2:3435630/100‑1 (MQ=255)
        ccTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGAtt                                                                                        >  2:662070/1‑100 (MQ=255)
         cTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGAttt                                                                                       >  2:115871/1‑100 (MQ=255)
          tGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTa                                                                                      <  1:4027535/100‑1 (MQ=255)
          tGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTa                                                                                      >  1:1685786/1‑100 (MQ=255)
          tGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTa                                                                                      >  1:286206/1‑100 (MQ=255)
           gAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTaa                                                                                     >  2:189340/1‑100 (MQ=255)
            aaaaCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGTGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAAc                                                                                    <  1:2640959/100‑1 (MQ=255)
            aaaaCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAAc                                                                                    <  1:3819867/100‑1 (MQ=255)
            aaaaCGCCCGCAACAAAACCCTCCTTGAGAACTTAACCCCGCTGCACGCAACCTCTCCTCTGCGTATGGACCGTGGTAACGGTTCGACTGAAGATTTAAc                                                                                    <  2:1513879/100‑1 (MQ=255)
              aaCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTg                                                                                  >  2:2548963/1‑100 (MQ=255)
                cGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCt                                                                                >  1:2253985/1‑100 (MQ=255)
                cGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCt                                                                                <  2:1377379/100‑1 (MQ=255)
                 gCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTc                                                                               >  2:2531989/1‑100 (MQ=255)
                  cccGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTcg                                                                              >  1:3840558/1‑100 (MQ=255)
                   ccGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTcgc                                                                             >  2:264546/1‑100 (MQ=255)
                      caacaaAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTa                                                                          >  2:3183350/1‑100 (MQ=255)
                         caaAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTg                                                                       >  1:652518/1‑100 (MQ=255)
                          aaaaTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTgg                                                                      <  2:1557895/100‑1 (MQ=255)
                          aaaaTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTgg                                                                      <  2:79800/100‑1 (MQ=255)
                          aaaaTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTgg                                                                      <  1:1960925/100‑1 (MQ=255)
                           aaaTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGa                                                                     >  1:1669614/1‑100 (MQ=255)
                            aaTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGAt                                                                    >  1:2731312/1‑100 (MQ=255)
                             aTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATc                                                                   <  2:905474/100‑1 (MQ=255)
                              tCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCt                                                                  <  1:2276422/100‑1 (MQ=255)
                               ccTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCtn                                                                 <  1:1279251/100‑2 (MQ=255)
                               ccTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTg                                                                 <  2:3221650/100‑1 (MQ=255)
                                ctctTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTgg                                                                >  1:362349/1‑100 (MQ=255)
                                ctctTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTgg                                                                >  2:1150584/1‑100 (MQ=255)
                                 tctTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGc                                                               >  1:208312/1‑100 (MQ=255)
                                 tctTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGc                                                               <  1:3345485/100‑1 (MQ=255)
                                    ttGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATc                                                            <  1:2905362/100‑1 (MQ=255)
                                     tGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCa                                                           <  1:3613225/100‑1 (MQ=255)
                                       agaACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCAcc                                                         <  2:2036917/100‑1 (MQ=255)
                                       agaACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCAcc                                                         <  1:300267/100‑1 (MQ=255)
                                           cTTACCCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATc                                                     <  1:3659574/100‑1 (MQ=255)
                                           cTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTAtn                                                     <  1:3253529/100‑2 (MQ=255)
                                           cTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATc                                                     <  2:1101677/100‑1 (MQ=255)
                                           cTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATc                                                     >  1:4026407/1‑100 (MQ=255)
                                              aaCCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGt                                                  <  2:367947/100‑1 (MQ=255)
                                                ccccGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCg                                                >  1:3142267/1‑100 (MQ=255)
                                                ccccGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCg                                                <  2:3858982/100‑1 (MQ=255)
                                                ccccGCTGCACGCAAACTCTCGTCTGCGTATGGAACGCGGTAACGGTTCAACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGtag                                                >  2:1827463/1‑98 (MQ=255)
                                                 cccGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGt                                               <  1:3703092/100‑1 (MQ=255)
                                                  ccGCGGCACGCAAACTCTCGTCTGCGTATGGAACGAGGTAACGGTTCGGCTGAACATTTCACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTg                                              <  2:758400/100‑1 (MQ=255)
                                                   cGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTgg                                             <  2:777784/100‑1 (MQ=255)
                                                   cGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTgg                                             <  1:2786847/100‑1 (MQ=255)
                                                          cGCAAAGTCTCGTCGGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGt                                      <  2:555532/100‑1 (MQ=255)
                                                          cGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGt                                      >  2:1296976/1‑100 (MQ=255)
                                                          cGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGt                                      <  1:1692282/100‑1 (MQ=255)
                                                          cGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGt                                      >  1:2201509/1‑100 (MQ=255)
                                                           gCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTg                                     >  1:3227391/1‑100 (MQ=255)
                                                           gCAAACTCTCGTCAGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCACGCGTACTGGATCTGGCATCACCTATCGATCGTGGTCAGCGta                                     >  2:203300/1‑99 (MQ=255)
                                                            cAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTgg                                    >  2:3492421/1‑100 (MQ=255)
                                                              agctctCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTc                                  <  1:46057/98‑1 (MQ=255)
                                                                ctctCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTg                                >  1:1042592/1‑100 (MQ=255)
                                                                ctctCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTg                                <  1:1784533/100‑1 (MQ=255)
                                                                 tctcGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGa                               >  1:791287/1‑100 (MQ=255)
                                                                 tctcGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGa                               >  1:790635/1‑100 (MQ=255)
                                                                  ctcGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGAt                              >  2:1009674/1‑100 (MQ=255)
                                                                  ctcGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGAt                              <  1:923622/100‑1 (MQ=255)
                                                                   tcGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGAtt                             <  1:2514438/100‑1 (MQ=255)
                                                                   gagtttGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGAtt                             <  1:1321059/95‑1 (MQ=255)
                                                                    cGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtg                            >  2:1317108/1‑100 (MQ=255)
                                                                     gTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtgt                           <  2:2924312/100‑1 (MQ=255)
                                                                       cTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTgg                         >  2:2883290/1‑100 (MQ=255)
                                                                        tGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGc                        <  2:2432093/100‑1 (MQ=255)
                                                                         gCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCa                       <  1:1341353/100‑1 (MQ=255)
                                                                          cGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAc                      >  1:929630/1‑100 (MQ=255)
                                                                           gTATGGAACGGGGTAACGGTTCGACTGAAGATTTAAGTGCTCGCGTACTGGATCTGGCGTCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAcc                     >  2:141875/1‑100 (MQ=255)
                                                                            tATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccg                    <  2:2717687/100‑1 (MQ=255)
                                                                            tATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGATCTGATTGTGGCAccg                    >  2:3352668/1‑100 (MQ=255)
                                                                             acggAACGTGGTACCGGTTCGAATCAAGCTTTAACTGCTCGAGTTCTGAATCTTGCATCACCTTTCGCTCGTGGTCAGCGTGGTCAGACTGTGGCAccgc                   <  2:1578186/98‑1 (MQ=255)
                                                                             aTGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccgc                   >  1:348332/1‑100 (MQ=255)
                                                                               ggAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccgccg                 >  2:3005760/1‑100 (MQ=255)
                                                                               ggAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccgccg                 >  1:1524895/1‑100 (MQ=255)
                                                                                   cGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAg             <  2:3045630/100‑1 (MQ=255)
                                                                                      ggTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCg          >  2:1713701/1‑100 (MQ=255)
                                                                                      ggTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCg          >  1:3226809/1‑100 (MQ=255)
                                                                                       gTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCgg         <  1:1634519/100‑1 (MQ=255)
                                                                                        tAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGt        <  1:2243074/100‑1 (MQ=255)
                                                                                         aaCGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTa       <  1:329368/100‑1 (MQ=255)
                                                                                          aCGGTTCGACTGAAGAGTTAACTGCTTGAGGACTGGATCTGGCATCACCTATCGATCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTaa      <  2:3512116/100‑1 (MQ=255)
                                                                                           cGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTaaa     >  2:3534465/1‑100 (MQ=255)
                                                                                           cGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTaaa     <  1:1895725/100‑1 (MQ=255)
                                                                                            ggTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTaaaa    <  1:2137203/100‑1 (MQ=255)
                                                                                              ttCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAAcc  >  1:2157359/1‑100 (MQ=255)
                                                                                              ttCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAAcc  >  2:3248536/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACC  >  NC_000913/3966778‑3966971

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: