Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,965,229 A→G I400V (ATT→GTT)  rho → transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,965,2290AG97.4% 949.6 / 8.4 303I400V (ATT→GTT) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base A (6/2);  new base G (193/102);  total (199/104)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-01

Only 100 of 307 total aligned reads displayed.

TGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGT  >  NC_000913/3965140‑3965318
                                                                                         |                                                                                         
tgctCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTta                                                                                        >  2:1340183/1‑93 (MQ=255)
 gctCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTtaa                                                                                       <  1:1902080/93‑1 (MQ=255)
 gctCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTtaa                                                                                       >  1:4333992/1‑93 (MQ=255)
    cACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTtaataa                                                                                    >  2:3542219/1‑93 (MQ=255)
     aCGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATaaa                                                                                   <  2:635829/93‑1 (MQ=255)
      cGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAAc                                                                                  <  1:4670840/93‑1 (MQ=255)
        aCTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTg                                                                                <  2:1606799/93‑1 (MQ=255)
        aCTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTg                                                                                >  1:1002831/1‑93 (MQ=255)
        aCTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTg                                                                                >  1:2785608/1‑93 (MQ=255)
         cTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTgg                                                                               >  1:3777702/1‑93 (MQ=255)
         cTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTgg                                                                               >  1:4208720/1‑93 (MQ=255)
         cTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTgg                                                                               <  1:4022012/93‑1 (MQ=255)
          tCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGc                                                                              >  2:2670074/1‑93 (MQ=255)
          tCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGc                                                                              <  2:3716741/93‑1 (MQ=255)
             ggAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTTGTTAATAAACTGGCAAt                                                                           <  1:2292408/93‑1 (MQ=255)
             ggAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAAt                                                                           <  2:2169600/93‑1 (MQ=255)
             ggAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAAt                                                                           >  1:3516627/1‑93 (MQ=255)
              gaagaaCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATg                                                                          >  2:2438505/1‑93 (MQ=255)
              gaagaaCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATg                                                                          >  2:2376080/1‑93 (MQ=255)
              gaagaaCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATg                                                                          >  2:666908/1‑93 (MQ=255)
              gaagaaCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATg                                                                          >  1:4103572/1‑93 (MQ=255)
                 gaaCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgacc                                                                       >  1:2557106/1‑93 (MQ=255)
                 gaaCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgacc                                                                       >  1:2955803/1‑93 (MQ=255)
                 gaaCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaca                                                                       >  1:1457843/1‑92 (MQ=255)
                   aCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaccaa                                                                     >  1:3089141/1‑93 (MQ=255)
                     tGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaccaaga                                                                   >  1:3385980/1‑93 (MQ=255)
                       cAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaccaagacc                                                                 >  2:3890652/1‑93 (MQ=255)
                       cAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaccaagacc                                                                 >  2:2458275/1‑93 (MQ=255)
                         gAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaccaagaccaa                                                               <  2:4368762/93‑1 (MQ=255)
                          aaaaTGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAAt                                                              <  1:3430287/93‑1 (MQ=255)
                            aaTGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACAAAGACCAATGa                                                            >  1:935641/1‑93 (MQ=255)
                             aTGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGAc                                                           <  1:349049/93‑1 (MQ=255)
                              tgtgGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACg                                                          >  1:2630691/1‑93 (MQ=255)
                               gtgGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGa                                                         >  2:2532701/1‑93 (MQ=255)
                               gtgGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGa                                                         >  1:3755158/1‑93 (MQ=255)
                               gtgGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGa                                                         <  1:2391605/93‑1 (MQ=255)
                               gtgGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGa                                                         <  1:1400496/93‑1 (MQ=255)
                               gtgGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGa                                                         >  2:1636958/1‑93 (MQ=255)
                                 ggATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGAtt                                                       >  1:3642911/1‑93 (MQ=255)
                                  gATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGAttt                                                      <  1:3826439/93‑1 (MQ=255)
                                  gATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGAttt                                                      <  1:4273729/93‑1 (MQ=255)
                                      cTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATttcttc                                                  >  1:78375/1‑93 (MQ=255)
                                          gcAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGaaa                                              <  2:1561792/93‑1 (MQ=255)
                                          gcAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGaaa                                              <  2:4609986/93‑1 (MQ=255)
                                           cAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGATCAATGACGATTTCTTCGAAat                                             <  2:4672280/93‑1 (MQ=255)
                                           cAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAat                                             <  2:4531672/93‑1 (MQ=255)
                                           cAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAat                                             <  2:3119768/93‑1 (MQ=255)
                                            aaaaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatg                                            >  2:3397940/1‑93 (MQ=255)
                                            aaaaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatg                                            >  2:3335011/1‑93 (MQ=255)
                                            aaaaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGAGTTCTTCGAAatg                                            >  2:2997903/1‑93 (MQ=255)
                                             aaaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatga                                           >  1:1318384/1‑93 (MQ=255)
                                              aaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatgat                                          >  1:1809375/1‑93 (MQ=255)
                                              aaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatgat                                          >  1:4392511/1‑93 (MQ=255)
                                              aaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatgat                                          >  1:672887/1‑93 (MQ=255)
                                                tCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAAtgatga                                        <  1:1134133/93‑1 (MQ=255)
                                                tCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAAtgatga                                        >  2:2720431/1‑93 (MQ=255)
                                                 cATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGaa                                       >  1:3335796/1‑93 (MQ=255)
                                                  aTTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGaaa                                      >  1:4312928/1‑93 (MQ=255)
                                                      aCCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCt                                  <  2:1816386/93‑1 (MQ=255)
                                                       cccGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTc                                 >  2:442364/1‑93 (MQ=255)
                                                       cccGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTc                                 <  1:323098/93‑1 (MQ=255)
                                                        ccGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCa                                >  1:2256372/1‑93 (MQ=255)
                                                        ccGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCa                                >  1:3904663/1‑93 (MQ=255)
                                                         cGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCAt                               <  2:1478772/93‑1 (MQ=255)
                                                         cGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCAt                               >  2:3984083/1‑93 (MQ=255)
                                                         cGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCAt                               >  2:3189479/1‑93 (MQ=255)
                                                         cGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCAt                               <  1:617133/93‑1 (MQ=255)
                                                            tGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATaaa                            >  2:3158592/1‑93 (MQ=255)
                                                               gCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAAttt                         <  2:765064/93‑1 (MQ=255)
                                                                 gAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATACATTTGt                       >  1:716290/1‑93 (MQ=255)
                                                                  aaaTCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTc                      >  2:1774684/1‑93 (MQ=255)
                                                                  aaaTCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTc                      >  2:3628788/1‑93 (MQ=255)
                                                                   aaTCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCt                     >  2:2221706/1‑93 (MQ=255)
                                                                    aTCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCtt                    >  2:191775/1‑93 (MQ=255)
                                                                    aTCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCtt                    >  2:864877/1‑93 (MQ=255)
                                                                    aTCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCtt                    >  1:1934493/1‑93 (MQ=255)
                                                                    aTCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCtt                    >  2:1169835/1‑93 (MQ=255)
                                                                     tCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTa                   <  1:4549563/93‑1 (MQ=255)
                                                                     tCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTa                   >  2:1516543/1‑93 (MQ=255)
                                                                       gATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATg                 >  1:1955522/1‑93 (MQ=255)
                                                                       gATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATg                 >  2:360771/1‑93 (MQ=255)
                                                                        attcAATGGAATTCCTCGTTAATAAACTGGCAAGGACCAAGACCAATGACGATTTCTTCGAAATAATGAAACGCTCATAAATTTGGCTTATGc                <  2:4119798/90‑1 (MQ=255)
                                                                        aTGNAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATNNAACGCTCATAAATTTGTCTTNTGc                >  2:4720467/1‑93 (MQ=255)
                                                                        aTGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGCAATGATGAAACGCTCATAAATTTGTCTTATGc                >  1:1968701/1‑93 (MQ=255)
                                                                        aTGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGc                >  1:2722715/1‑93 (MQ=255)
                                                                            aaTGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCaaa            >  1:2030532/1‑93 (MQ=255)
                                                                             aTGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCaaaa           >  1:1203809/1‑93 (MQ=255)
                                                                             aTGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCaaaa           >  2:3513603/1‑93 (MQ=255)
                                                                               ggAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAAc         <  1:842474/93‑1 (MQ=255)
                                                                               ggAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAAc         <  2:2084189/93‑1 (MQ=255)
                                                                               ggAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAAc         >  1:3176309/1‑93 (MQ=255)
                                                                                gAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACg        <  2:3676350/93‑1 (MQ=255)
                                                                                gAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACg        >  2:1423364/1‑93 (MQ=255)
                                                                                  aTTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGcc      >  1:2878635/1‑93 (MQ=255)
                                                                                  aTTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGcc      >  2:972885/1‑93 (MQ=255)
                                                                                    tCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCAc    <  2:3828301/93‑1 (MQ=255)
                                                                                    tCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCAc    >  1:3202917/1‑93 (MQ=255)
                                                                                    tCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCAc    >  2:798706/1‑93 (MQ=255)
                                                                                      cTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACgt  >  1:3059834/1‑93 (MQ=255)
                                                                                      cTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACgt  >  1:1745736/1‑93 (MQ=255)
                                                                                         |                                                                                         
TGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGT  >  NC_000913/3965140‑3965318

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