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,8750TG98.2% 765.9 / 2.5 225S153S (TCT→TCG) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base T (0/4);  new base G (97/123);  total (97/128)
Fisher's exact test for biased strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

Only 100 of 225 total aligned reads displayed.

TCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAAC  >  NC_000913/3966778‑3966970
                                                                                                 |                                                                                               
tCGACAAACCGGAAAACGCCCGGATCAAAATCCTCGGTGTGATCTTAACCCCGCAGCACGCAAGCTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGAc                                                                                               <  2:688237/100‑1 (MQ=255)
 cGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACt                                                                                              <  2:396489/100‑1 (MQ=255)
 cGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACt                                                                                              <  2:4982712/100‑1 (MQ=255)
 cGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACt                                                                                              <  2:3548645/100‑1 (MQ=255)
 cGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACt                                                                                              <  2:1126092/100‑1 (MQ=255)
   aCAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGa                                                                                            >  1:307483/1‑100 (MQ=255)
   aCAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGa                                                                                            >  2:1517415/1‑100 (MQ=255)
    cAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGaa                                                                                           >  1:69429/1‑100 (MQ=255)
    cAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGAGTATGGAACGTGGTAACGGTTCGACTGaa                                                                                           >  1:5083303/1‑100 (MQ=255)
     aaaCCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAg                                                                                          <  1:1017984/100‑1 (MQ=255)
      aaCCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGa                                                                                         >  1:4610079/1‑100 (MQ=255)
        ccTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGAtt                                                                                       <  2:369592/100‑1 (MQ=255)
          tGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTACCCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTa                                                                                     >  2:5373475/1‑100 (MQ=255)
          tGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTa                                                                                     >  2:2388041/1‑100 (MQ=255)
            aaaaCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAAc                                                                                   <  1:3643884/100‑1 (MQ=255)
            aaaaCGCCCGCAACAAAATCCTCGTTGAGTACTTAACCGCGCTGCACGCAAACTCTCGTCTGCGTATGTAACGTGGTAACGGTTCGACTGAAGATTTAAc                                                                                   <  2:3606269/100‑1 (MQ=255)
             caaCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCTGCTGTACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACt                                                                                  <  2:2838565/99‑1 (MQ=255)
              aaCGCCCGCAACAAAATCCACTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGAGTGAAGATTTAACTg                                                                                 >  2:5032581/1‑100 (MQ=255)
                   ccGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTcgc                                                                            >  2:3849638/1‑100 (MQ=255)
                     gCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGt                                                                          <  2:5333465/100‑1 (MQ=255)
                     gCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGt                                                                          <  1:2113105/100‑1 (MQ=255)
                     gCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGt                                                                          >  1:2678519/1‑100 (MQ=255)
                      caacaaAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTa                                                                         >  2:4550616/1‑100 (MQ=255)
                      caacaaAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTa                                                                         >  1:1735871/1‑100 (MQ=255)
                        acaaAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACt                                                                       >  2:2430184/1‑100 (MQ=255)
                         caaAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTg                                                                      <  2:4859833/100‑1 (MQ=255)
                         caaAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTg                                                                      >  1:3247527/1‑100 (MQ=255)
                          aaaaTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTgg                                                                     <  1:732788/100‑1 (MQ=255)
                             aTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATc                                                                  <  1:2562514/100‑1 (MQ=255)
                             aTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATc                                                                  >  1:3896016/1‑100 (MQ=255)
                              tCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCt                                                                 <  1:1536645/100‑1 (MQ=255)
                              tCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCt                                                                 >  2:3750755/1‑100 (MQ=255)
                              tCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCt                                                                 >  1:2642109/1‑100 (MQ=255)
                               ccTCTTTGAGATCTTAACACCGCTGTACGCAACCCTTCGCCTGCCTATGGAACGTGTTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATTTg                                                                <  2:1353557/100‑1 (MQ=255)
                                 tctTTGAGAACTTAGCCACGCTGCACGCAAACTCGCGTTTGCGTATGGAACGTGGTAACGGTGCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGc                                                              <  2:5131246/100‑1 (MQ=255)
                                 tctTTGAGAACTTAACCCCGCTGCACGCAAACTTTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGc                                                              <  1:5507991/100‑1 (MQ=255)
                                  ctTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGTTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCa                                                             <  1:3567883/100‑1 (MQ=255)
                                  ctTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCa                                                             >  2:2541582/1‑100 (MQ=255)
                                  ctTTGAGAACTGAACCCCGCTGCACGCAAACTCTCGTCTTCGTCTGGAACGTGGTAACGGTACGACTAAAGATTGAACTGCTCGCGTACTGGATCTGGCa                                                             >  2:5009964/1‑100 (MQ=255)
                                   tttGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCAt                                                            <  2:4883530/100‑1 (MQ=255)
                                   tttGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCAt                                                            <  2:304600/100‑1 (MQ=255)
                                    ttGAGCACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATc                                                           <  2:5500567/100‑1 (MQ=255)
                                      gagaACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCan                                                         <  1:4435407/100‑2 (MQ=255)
                                      gagaACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCAc                                                         <  2:4168477/100‑1 (MQ=255)
                                       agaACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCAcc                                                        <  1:1445693/100‑1 (MQ=255)
                                        gaACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCt                                                       <  2:4665784/100‑1 (MQ=255)
                                          aCTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTAt                                                     <  2:4442802/100‑1 (MQ=255)
                                           cTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATc                                                    <  2:1275745/100‑1 (MQ=255)
                                             tAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCgg                                                  >  1:807369/1‑100 (MQ=255)
                                               aCCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTc                                                >  2:3322453/1‑100 (MQ=255)
                                                ccccGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCg                                               >  1:1246372/1‑100 (MQ=255)
                                                ccccGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCg                                               >  2:2903513/1‑100 (MQ=255)
                                                 cccGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGt                                              <  1:4962575/100‑1 (MQ=255)
                                                 cccGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGAGTTAACTGCTCGCGTACTGGATCTGGCGTCACCTATCGGTCGt                                              <  2:1862801/100‑1 (MQ=255)
                                                   cGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTgg                                            <  1:5484552/100‑1 (MQ=255)
                                                    gCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGt                                           >  2:2186420/1‑100 (MQ=255)
                                                     gatctcGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTc                                          <  1:1377424/95‑1 (MQ=255)
                                                        cACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGc                                       >  1:5131708/1‑100 (MQ=255)
                                                        cACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGc                                       <  1:5103024/100‑1 (MQ=255)
                                                          cGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGt                                     >  2:2737781/1‑100 (MQ=255)
                                                           gCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTg                                    <  1:4266587/100‑1 (MQ=255)
                                                            caagttctcGTCTTCGTATAGAACGTGCGAACGGTTCGACTGAAGATTTAACTGGTCGCGCACTGGATCTGGCATCCCCTATCGGTCGTGGTCAGCGTgg                                   <  2:945819/95‑1 (MQ=255)
                                                            cAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTgg                                   >  1:2000857/1‑100 (MQ=255)
                                                            cAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTgg                                   >  2:316002/1‑100 (MQ=255)
                                                             aaaCTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGt                                  >  1:2847383/1‑100 (MQ=255)
                                                               aCTCTCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCt                                <  1:4765402/100‑1 (MQ=255)
                                                                ctctCGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTg                               >  1:2063615/1‑100 (MQ=255)
                                                                  ctcGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGAt                             >  1:1087546/1‑100 (MQ=255)
                                                                   tcGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGAtt                            <  2:3079639/100‑1 (MQ=255)
                                                                    cGTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtg                           <  1:1895665/100‑1 (MQ=255)
                                                                     gTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtgt                          >  1:3702282/1‑100 (MQ=255)
                                                                     gTCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtgt                          >  2:3261478/1‑100 (MQ=255)
                                                                      tCTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtgtg                         <  2:5471549/100‑1 (MQ=255)
                                                                       cTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTgg                        >  2:1166788/1‑100 (MQ=255)
                                                                       cTGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTgg                        <  1:2993003/100‑1 (MQ=255)
                                                                       cTGCGTATGGAACGTCGTAACGGTTCCACTGAAGATTTAACTGCCTGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTgg                        <  2:5239328/100‑1 (MQ=255)
                                                                        tGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGc                       <  1:2282586/100‑1 (MQ=255)
                                                                        tGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGc                       <  1:2709056/100‑1 (MQ=255)
                                                                        tGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGc                       <  1:4509305/100‑1 (MQ=255)
                                                                        tGCGTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGc                       <  1:1401396/100‑1 (MQ=255)
                                                                           gTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAcc                    >  1:1603602/1‑100 (MQ=255)
                                                                           gTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAcc                    >  2:1611058/1‑100 (MQ=255)
                                                                           gTATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAcc                    <  1:2170235/100‑1 (MQ=255)
                                                                            tATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccg                   <  1:1298861/100‑1 (MQ=255)
                                                                            tATGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccg                   >  2:4324214/1‑100 (MQ=255)
                                                                             aTGGAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccgc                  >  1:3146360/1‑100 (MQ=255)
                                                                               ggAACGTGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccgccg                <  1:4115260/100‑1 (MQ=255)
                                                                                  aCGTGATAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGACGTGGTCAGTGTGTTCTGATTGTGGCACCGCCGaaa             <  2:877123/100‑1 (MQ=255)
                                                                                     tGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGcc          >  2:1747258/1‑100 (MQ=255)
                                                                                     tGGTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGcc          <  1:5410481/100‑1 (MQ=255)
                                                                                      ggTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCg         >  2:2051382/1‑100 (MQ=255)
                                                                                      ggTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCg         >  2:4003253/1‑100 (MQ=255)
                                                                                       gTAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGccag        >  2:1424221/1‑98 (MQ=255)
                                                                                        tAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGt       >  2:1585561/1‑100 (MQ=255)
                                                                                        tAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGt       >  1:3740999/1‑100 (MQ=255)
                                                                                        tAACGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGt       >  1:989384/1‑100 (MQ=255)
                                                                                         aaCGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTa      <  1:3353655/100‑1 (MQ=255)
                                                                                         aaCGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTa      <  1:678479/100‑1 (MQ=255)
                                                                                          aCGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTaa     >  2:3424628/1‑100 (MQ=255)
                                                                                           cGGTTCGACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTaaa    >  1:921509/1‑100 (MQ=255)
                                                                                                 |                                                                                               
TCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAAC  >  NC_000913/3966778‑3966970

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