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,2290AG98.9% 829.7 / ‑1.5 264I400V (ATT→GTT) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (141/119);  total (143/121)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

Only 100 of 267 total aligned reads displayed.

CTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGT  >  NC_000913/3965139‑3965318
                                                                                          |                                                                                         
 tgctCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTta                                                                                        >  2:1145883/1‑93 (MQ=255)
    tCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTtaata                                                                                     <  1:1950803/93‑1 (MQ=255)
    tCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTtaata                                                                                     <  1:2173018/93‑1 (MQ=255)
      aCGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATaaa                                                                                   <  1:1265416/93‑1 (MQ=255)
      aCGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATaaa                                                                                   <  1:1365033/93‑1 (MQ=255)
      aCGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATaaa                                                                                   <  2:4336437/93‑1 (MQ=255)
      aCGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATaaa                                                                                   <  2:2277463/93‑1 (MQ=255)
      aCGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATaaa                                                                                   <  1:3238906/93‑1 (MQ=255)
      aCGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATaaa                                                                                   <  2:1718829/93‑1 (MQ=255)
        gACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACt                                                                                 <  1:888507/93‑1 (MQ=255)
        gACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACt                                                                                 >  2:2851618/1‑93 (MQ=255)
         aCTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTg                                                                                >  1:3404849/1‑93 (MQ=255)
         aCTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTg                                                                                >  2:2164597/1‑93 (MQ=255)
           tCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGc                                                                              >  2:3135768/1‑93 (MQ=255)
           tCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGc                                                                              >  2:3645610/1‑93 (MQ=255)
              ggAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAAt                                                                           >  2:739924/1‑93 (MQ=255)
              ggAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAAt                                                                           >  1:530289/1‑93 (MQ=255)
              ggAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAAt                                                                           <  2:474913/93‑1 (MQ=255)
               gaagaaCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATg                                                                          >  2:1830935/1‑93 (MQ=255)
                 agaaCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgac                                                                        >  2:3344257/1‑93 (MQ=255)
                 agaaCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgac                                                                        >  2:133143/1‑93 (MQ=255)
                 agaaCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgac                                                                        >  2:3988059/1‑93 (MQ=255)
                  gaaCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgacc                                                                       >  1:1768096/1‑93 (MQ=255)
                    aCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaccaa                                                                     >  1:2982638/1‑93 (MQ=255)
                     cTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaccaag                                                                    >  2:2081824/1‑93 (MQ=255)
                     cTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaccaag                                                                    >  1:3872044/1‑93 (MQ=255)
                         aGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaccaagacca                                                                <  2:3458373/93‑1 (MQ=255)
                          gAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaccaagaccaa                                                               >  1:857603/1‑93 (MQ=255)
                          gAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaccaagaccaa                                                               >  2:883349/1‑93 (MQ=255)
                          gAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaccaagaccaa                                                               <  1:4655199/93‑1 (MQ=255)
                          gAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaccaagaccaa                                                               >  1:2925836/1‑93 (MQ=255)
                           aaaaTGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAAt                                                              <  2:1642196/93‑1 (MQ=255)
                           aaaaTGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAAt                                                              <  2:375444/93‑1 (MQ=255)
                             aaTGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGa                                                            >  2:3487980/1‑93 (MQ=255)
                              aTGTGGATCCTGCGGAAAATCAGTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGAc                                                           <  2:1560055/93‑1 (MQ=255)
                              aTGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATCACCAAGACCAATGAc                                                           <  2:1581684/93‑1 (MQ=255)
                               tgtgGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACg                                                          >  2:3376933/1‑93 (MQ=255)
                                gtgGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGa                                                         >  2:618580/1‑93 (MQ=255)
                                 tgGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGAt                                                        >  2:4079744/1‑93 (MQ=255)
                                  ggATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGAtt                                                       <  2:2381630/93‑1 (MQ=255)
                                   gATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGAttt                                                      <  1:4431056/93‑1 (MQ=255)
                                   gATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGAttt                                                      >  1:658749/1‑93 (MQ=255)
                                       cTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATttcttc                                                  >  2:814982/1‑93 (MQ=255)
                                        tGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCg                                                 >  2:4507537/1‑93 (MQ=255)
                                         gcgcAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGa                                                <  1:206858/93‑1 (MQ=255)
                                           gcAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGaaa                                              <  2:2454603/93‑1 (MQ=255)
                                           gcAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGaaa                                              <  2:2691135/93‑1 (MQ=255)
                                            cAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAat                                             <  1:1569835/93‑1 (MQ=255)
                                            cAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAat                                             >  1:2701922/1‑93 (MQ=255)
                                            cAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAat                                             <  1:1074458/93‑1 (MQ=255)
                                             aaaaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatg                                            >  1:2040402/1‑93 (MQ=255)
                                              aaaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatga                                           >  1:2812282/1‑93 (MQ=255)
                                              aaaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatga                                           >  2:2602614/1‑93 (MQ=255)
                                              aaaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatga                                           >  2:1210488/1‑93 (MQ=255)
                                               aaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatgat                                          <  1:4489067/93‑1 (MQ=255)
                                                aTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatgatg                                         >  1:4609498/1‑93 (MQ=255)
                                                 tCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAAtgatga                                        >  1:4267858/1‑93 (MQ=255)
                                                 tCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAAtgatga                                        <  2:4289583/93‑1 (MQ=255)
                                                  cATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGaa                                       >  1:4629953/1‑93 (MQ=255)
                                                   aTTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGaaa                                      <  1:4077470/93‑1 (MQ=255)
                                                   aTTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGaaa                                      <  1:3741221/93‑1 (MQ=255)
                                                    ttCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAAc                                     >  2:3916013/1‑93 (MQ=255)
                                                         ccGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCa                                <  1:400601/93‑1 (MQ=255)
                                                         ccGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCa                                <  1:3585464/93‑1 (MQ=255)
                                                         ccGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCa                                >  2:3566369/1‑93 (MQ=255)
                                                             tGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATaaa                            >  1:2634775/1‑93 (MQ=255)
                                                             tGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATaaa                            <  1:2895780/93‑1 (MQ=255)
                                                             tGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATaaa                            >  1:4191259/1‑93 (MQ=255)
                                                             tGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATaaa                            >  2:3763207/1‑93 (MQ=255)
                                                               ggCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAAtt                          <  2:3689626/93‑1 (MQ=255)
                                                                 cGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTg                        <  1:3329640/93‑1 (MQ=255)
                                                                  gAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGt                       >  2:573794/1‑93 (MQ=255)
                                                                   aaaTCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTc                      <  2:3086643/93‑1 (MQ=255)
                                                                    aaTCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGACACGCTCATAAATTTGTCt                     >  1:3503860/1‑93 (MQ=255)
                                                                      tCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTa                   <  1:4322376/93‑1 (MQ=255)
                                                                       ataTGCAATGGAATACCTCGTTAATAAACTGGGAATGACCAAGACCAATGACGTTTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTAt                  <  2:709546/91‑1 (MQ=255)
                                                                        gATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATg                 <  1:1064502/93‑1 (MQ=255)
                                                                         aTGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGc                >  2:3060858/1‑93 (MQ=255)
                                                                          tGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGcc               >  2:1235476/1‑93 (MQ=255)
                                                                           gCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCa              <  2:3598577/93‑1 (MQ=255)
                                                                             aaTGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCaaa            <  1:2574776/93‑1 (MQ=255)
                                                                              aTGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCaaaa           >  2:1258054/1‑93 (MQ=255)
                                                                              aTGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCaaaa           <  2:200221/93‑1 (MQ=255)
                                                                               tGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCaaaca          >  2:1087109/1‑91 (MQ=255)
                                                                               tGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCaaaaa          >  2:4201901/1‑93 (MQ=255)
                                                                                ggAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAAc         <  2:2482117/93‑1 (MQ=255)
                                                                                  aaTTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGc       <  1:3840367/93‑1 (MQ=255)
                                                                                   aTTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGcc      >  1:144058/1‑93 (MQ=255)
                                                                                   aTTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGcc      <  1:2490124/93‑1 (MQ=255)
                                                                                    ttCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCTAAAAACGCCa     >  1:2438729/1‑93 (MQ=255)
                                                                                    ttCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCa     >  2:177237/1‑93 (MQ=255)
                                                                                    ttCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCa     <  1:3979970/93‑1 (MQ=255)
                                                                                    ttCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCa     >  2:3248635/1‑93 (MQ=255)
                                                                                    ttCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCa     >  2:1937773/1‑93 (MQ=255)
                                                                                    ttCCTCGTTAATAAACAGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCa     <  2:3806902/93‑1 (MQ=255)
                                                                                     tCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGACATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACACCag    >  2:4621544/1‑92 (MQ=255)
                                                                                     tCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCAc    >  1:2324411/1‑93 (MQ=255)
                                                                                       cTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACgt  >  2:2743349/1‑93 (MQ=255)
                                                                                       cTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACgt  >  2:1116204/1‑93 (MQ=255)
                                                                                       cTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACgt  <  2:4541258/93‑1 (MQ=255)
                                                                                          |                                                                                         
CTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGT  >  NC_000913/3965139‑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: