Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,966,721 G→A 20.5% R102H (CGC→CAC)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,966,7210GA20.5% 244.5 / 54.6 117R102H (CGC→CAC) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base G (51/41);  new base A (10/14);  total (62/55)
Fisher's exact test for biased strand distribution p-value = 2.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 128 total aligned reads displayed.

CAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAA  >  NC_000913/3966623‑3966820
                                                                                                  |                                                                                                   
cAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAACTCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAACATTcgc                                                                                                    <  2:3106810/100‑1 (MQ=255)
cAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgc                                                                                                    <  1:1177307/100‑1 (MQ=255)
 aGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcc                                                                                                   <  2:2143837/100‑1 (MQ=255)
  gCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGCTTcgccc                                                                                                  >  2:1594445/1‑100 (MQ=255)
  gCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGCTTcgccc                                                                                                  >  2:1594581/1‑100 (MQ=255)
  gCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgccc                                                                                                  >  1:2651751/1‑100 (MQ=255)
     ccTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcccgcc                                                                                               <  2:1278031/100‑1 (MQ=255)
     ccTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcccgcc                                                                                               >  2:1888924/1‑100 (MQ=255)
      cTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCccgccg                                                                                              >  1:798718/1‑100 (MQ=255)
         ccTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCAGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGaa                                                                                            >  2:622413/1‑99 (MQ=255)
          cTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGCTTCGCCCGCCGAaag                                                                                          <  1:671925/100‑1 (MQ=255)
          cTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAaag                                                                                          >  2:482850/1‑100 (MQ=255)
           tCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGCTTCGCCCGCCGAagga                                                                                         <  1:1710099/100‑1 (MQ=255)
           tCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGCTTCGCCCGCCGAaaga                                                                                         <  1:1710223/100‑1 (MQ=255)
            cGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAaagaa                                                                                        >  2:1697767/1‑100 (MQ=255)
                ggTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTg                                                                                    >  1:1339602/1‑100 (MQ=255)
                ggTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTg                                                                                    <  1:1927275/100‑1 (MQ=255)
                    cTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACg                                                                                >  1:2231225/1‑100 (MQ=255)
                    cTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACg                                                                                >  1:1012596/1‑100 (MQ=255)
                     tgatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGCTTCGCCCGCCGAAAGAAGGTGAACGc                                                                               >  2:2539298/1‑100 (MQ=255)
                     tgatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGc                                                                               <  1:1111579/100‑1 (MQ=255)
                     tgatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGc                                                                               >  1:2190452/1‑100 (MQ=255)
                     tgatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGc                                                                               <  1:800433/100‑1 (MQ=255)
                      gatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTAGCCCGCCGAAAAAAGGTGAACGc                                                                               >  2:183458/1‑99 (MQ=255)
                          aCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAttt                                                                          <  1:1645764/100‑1 (MQ=255)
                           cATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAttt                                                                          >  2:1375057/1‑99 (MQ=255)
                             tCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGCTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgc                                                                       >  1:2625357/1‑100 (MQ=255)
                             tCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgc                                                                       >  2:894390/1‑100 (MQ=255)
                                aCGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGCTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgct                                                                    <  2:1455117/100‑1 (MQ=255)
                                 cGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAGAGAAGGTGAACGCTATTTTGCgctg                                                                   >  1:3163911/1‑100 (MQ=255)
                                 cGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCACCACTGGTAAGCTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTtcg                                                                      >  2:851986/1‑97 (MQ=255)
                                  gTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGcgccgc                                                                  <  1:622415/100‑1 (MQ=255)
                                  gTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgc                                                                  <  1:79088/100‑1 (MQ=255)
                                   tttCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTCTTTTGCgctgct                                                                 <  1:2151805/100‑1 (MQ=255)
                                   tttCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgct                                                                 >  1:2650806/1‑100 (MQ=255)
                                    ttCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGctgctg                                                                <  2:600974/100‑1 (MQ=255)
                                     tCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGa                                                               >  2:1362432/1‑100 (MQ=255)
                                      ccccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaa                                                              <  1:1755036/100‑1 (MQ=255)
                                       cccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaaa                                                             >  2:755637/1‑100 (MQ=255)
                                       cccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaaa                                                             >  1:3268631/1‑100 (MQ=255)
                                        ccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAg                                                            <  2:149811/100‑1 (MQ=255)
                                        ccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAg                                                            >  1:55903/1‑100 (MQ=255)
                                        ccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAg                                                            >  2:2504325/1‑100 (MQ=255)
                                            gCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGCTTCGCCCGCCGAACGAAGGTGCAAGCTATTTTGAGCTGCTGCCAGTTaa                                                        >  1:2731735/1‑100 (MQ=255)
                                            gCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGCAAGTTaa                                                        >  2:266206/1‑100 (MQ=255)
                                            gCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTaa                                                        >  1:2724938/1‑100 (MQ=255)
                                            gCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTaa                                                        >  1:2020979/1‑100 (MQ=255)
                                                 aTCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAg                                                   <  2:3137003/100‑1 (MQ=255)
                                                 aTCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAg                                                   >  2:1169217/1‑100 (MQ=255)
                                                  tCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGt                                                  <  2:2136581/100‑1 (MQ=255)
                                                  tCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGt                                                  <  1:2065071/100‑1 (MQ=255)
                                                   ccgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGtt                                                 <  1:3032544/100‑1 (MQ=255)
                                                    cgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGCTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGtt                                                 >  2:857912/1‑99 (MQ=255)
                                                    cgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGCTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTa                                                >  1:2959222/1‑100 (MQ=255)
                                                    cgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTa                                                >  1:596503/1‑100 (MQ=255)
                                                       cgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACt                                             >  2:2651715/1‑100 (MQ=255)
                                                        gTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACtt                                            >  1:1678880/1‑100 (MQ=255)
                                                         tttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTc                                           <  1:1012554/100‑1 (MQ=255)
                                                         tttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTc                                           <  2:1384278/100‑1 (MQ=255)
                                                          ttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCg                                          >  2:487992/1‑100 (MQ=255)
                                                          ttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCg                                          <  1:954241/100‑1 (MQ=255)
                                                           tCAACCTCCGCACTGGTGATACCATCTCTGGTAAGCTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGa                                         >  2:2331577/1‑100 (MQ=255)
                                                           tCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGa                                         >  1:2769166/1‑100 (MQ=255)
                                                           tCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGa                                         <  1:1392461/100‑1 (MQ=255)
                                                            cAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                        >  2:1680313/1‑100 (MQ=255)
                                                            cAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                        >  2:2850788/1‑100 (MQ=255)
                                                             aaCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACa                                       >  2:64786/1‑100 (MQ=255)
                                                             aaCCTCCGCACTGGTGATACCAACTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGa                                         <  2:1625231/98‑1 (MQ=255)
                                                               ccTCCGCACTGGTGATACCATCTCTGGTAAGCTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaaa                                     >  2:173918/1‑100 (MQ=255)
                                                               ccTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaaa                                     >  1:3184351/1‑100 (MQ=255)
                                                                cTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAc                                    <  1:2488825/100‑1 (MQ=255)
                                                                 tCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGGAAGTTAACTTCGTTAACtttcagcagcg                                   <  2:2445071/100‑10 (MQ=37)
                                                                 tCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAcc                                   <  1:2099704/100‑1 (MQ=255)
                                                                   cGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTg                                 >  1:2466481/1‑100 (MQ=255)
                                                                   cGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTg                                 <  1:1637709/100‑1 (MQ=255)
                                                                   cGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTg                                 <  2:2318620/100‑1 (MQ=255)
                                                                    gCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGa                                <  1:2619959/100‑1 (MQ=255)
                                                                      accgctGATACCATCACTGGCAAGATTCTCCATCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAATTAACGAAGTTAACTTCGACAAACCTGaaa                              >  2:3059389/6‑100 (MQ=255)
                                                                      aCTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaa                              <  1:2504336/100‑1 (MQ=255)
                                                                         ggTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGACGGTGAACGCTATTTTTCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaaa                             >  2:3060974/1‑98 (MQ=255)
                                                                          gTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGc                          >  1:2910889/1‑100 (MQ=255)
                                                                            gATACCATCTCTGGTAAGATTCGCCCGCCGTAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGccc                        >  2:2789597/1‑100 (MQ=255)
                                                                              tACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                      >  1:1637496/1‑100 (MQ=255)
                                                                                ccATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaa                    <  2:1861963/100‑1 (MQ=255)
                                                                                 cATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaac                   <  2:1339599/100‑1 (MQ=255)
                                                                                 cATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaac                   <  1:482852/100‑1 (MQ=255)
                                                                                   tctcTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaacaa                 <  2:362920/100‑1 (MQ=255)
                                                                                     tctGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACaaaa               >  2:1618344/1‑100 (MQ=255)
                                                                                      ctGGTAAGCTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAt              >  2:1038451/1‑100 (MQ=255)
                                                                                        ggTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc            <  2:142155/100‑1 (MQ=255)
                                                                                        ggTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc            >  1:709957/1‑100 (MQ=255)
                                                                                        ggTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTATGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc            >  2:1753340/1‑100 (MQ=255)
                                                                                        ggTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc            >  1:1410769/1‑100 (MQ=255)
                                                                                           aaGCTTCGCCCGCCGAAAGAAGGTGAACGGTATTTTGCGCTGCTGAAAGTTAGCGAAGTTACCTTCGACAAACATGAAAACGCCCGTAACAAACTcatat         >  2:2778726/1‑96 (MQ=255)
                                                                                            agcttCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCtt        >  1:1241214/4‑100 (MQ=255)
                                                                                            aGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCtt        >  1:3022166/1‑100 (MQ=255)
                                                                                             gcttCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTCAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt       <  2:2295791/98‑1 (MQ=255)
                                                                                              aTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg      >  1:3130576/1‑100 (MQ=255)
                                                                                               ttCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTga     >  2:2075271/1‑100 (MQ=255)
                                                                                                 cgcccgccGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTgaga   >  1:561192/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAA  >  NC_000913/3966623‑3966820

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