Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,964,245 T→C Y72H (TAC→CAC)  rho → transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,964,2450TC99.4% 516.6 / ‑4.9 164Y72H (TAC→CAC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (77/86);  total (78/86)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01

Only 100 of 174 total aligned reads displayed.

ACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATT  >  NC_000913/3964156‑3964334
                                                                                         |                                                                                         
aCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCAcc                                                                                        >  2:3084362/1‑93 (MQ=255)
    aaaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGc                                                                                    >  1:2474627/1‑93 (MQ=255)
     aaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGca                                                                                   <  1:3059618/93‑2 (MQ=255)
       gagTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCgg                                                                                 >  2:385543/1‑93 (MQ=255)
         gTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTc                                                                               >  1:2579178/1‑93 (MQ=255)
           ggCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCt                                                                             <  2:775990/93‑1 (MQ=255)
            gCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCCGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCtg                                                                            <  1:1598710/93‑1 (MQ=255)
             cGAATATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCtga                                                                           <  1:4061706/93‑1 (MQ=255)
             cGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCtga                                                                           <  2:3305387/93‑1 (MQ=255)
              gAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCtgat                                                                          <  2:1949908/93‑1 (MQ=255)
              gAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCtgat                                                                          >  1:4154429/1‑93 (MQ=255)
              gAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCtgat                                                                          <  1:2644602/93‑1 (MQ=255)
               aaGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCtgatg                                                                         >  1:2083023/1‑93 (MQ=255)
                aGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCtgatga                                                                        >  2:2633847/1‑93 (MQ=255)
                  atatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACa                                                                      <  1:3895890/93‑1 (MQ=255)
                  atatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACa                                                                      <  1:3174445/93‑1 (MQ=255)
                   tatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACAt                                                                     <  1:1704120/93‑1 (MQ=255)
                    atCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATc                                                                    <  2:2006677/93‑1 (MQ=255)
                    atCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATc                                                                    >  2:4631505/1‑93 (MQ=255)
                    atCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATc                                                                    >  2:4628857/1‑93 (MQ=255)
                     tCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCt                                                                   <  2:311146/93‑1 (MQ=255)
                     tCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCt                                                                   <  1:366456/93‑1 (MQ=255)
                      cTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTa                                                                  <  2:3962921/93‑1 (MQ=255)
                      cTTTGGTGATGGCGTACTGGAGATATTGAAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTa                                                                  <  2:2707681/93‑1 (MQ=255)
                       tttGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTAc                                                                 >  1:3222849/1‑93 (MQ=255)
                        ttGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACg                                                                <  1:2997773/93‑1 (MQ=255)
                        ttGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACg                                                                <  1:17477/93‑1 (MQ=255)
                         tGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGt                                                               <  1:2897234/93‑1 (MQ=255)
                          ggTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGtt                                                              <  2:3323018/93‑1 (MQ=255)
                          ggTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGtt                                                              <  2:1692449/93‑1 (MQ=255)
                           gTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGttt                                                             >  2:3889307/1‑93 (MQ=255)
                            tGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTc                                                            >  1:1566505/1‑93 (MQ=255)
                            tGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTc                                                            >  2:1984868/1‑93 (MQ=255)
                            tGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTc                                                            <  1:1797127/93‑1 (MQ=255)
                             gATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTcc                                                           <  2:2659771/93‑1 (MQ=255)
                             gATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTcc                                                           <  1:3056494/93‑1 (MQ=255)
                                ggCGTACTGGATATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCt                                                        <  2:1652168/93‑1 (MQ=255)
                                ggCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCt                                                        >  2:4655545/1‑93 (MQ=255)
                                ggCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCt                                                        >  2:481842/1‑93 (MQ=255)
                                   gTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGc                                                     <  1:2347277/93‑1 (MQ=255)
                                   gTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGc                                                     >  1:38177/1‑93 (MQ=255)
                                   gTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGc                                                     >  1:1773156/1‑93 (MQ=255)
                                    tACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGcc                                                    <  2:962986/93‑1 (MQ=255)
                                    tACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGcc                                                    >  2:2957847/1‑93 (MQ=255)
                                       tGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCaaa                                                 >  1:2975420/1‑93 (MQ=255)
                                        ggAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAAt                                                <  1:75526/93‑1 (MQ=255)
                                        ggAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAAt                                                <  1:3317064/93‑1 (MQ=255)
                                         gagaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATc                                               >  2:2301531/1‑93 (MQ=255)
                                         gagaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATc                                               <  2:3874118/93‑1 (MQ=255)
                                         gagaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATc                                               <  2:1884389/93‑1 (MQ=255)
                                          agaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATcc                                              >  2:1432277/1‑93 (MQ=255)
                                            atatTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccgc                                            >  2:1835039/1‑93 (MQ=255)
                                             tatTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccgcc                                           <  1:1452485/93‑1 (MQ=255)
                                                tGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGtt                                        <  2:2111441/93‑1 (MQ=255)
                                                  cAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTc                                      <  1:1357342/93‑1 (MQ=255)
                                                    ggagggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCaa                                    <  2:2779399/93‑1 (MQ=255)
                                                     gatggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATACGACGTTTCAAc                                   >  1:2979715/1‑93 (MQ=255)
                                                      atggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAAcc                                  <  2:956577/93‑1 (MQ=255)
                                                        ggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCACCTAGCCAAATCCGCCGTTTCAAACTc                                >  2:44745/1‑93 (MQ=255)
                                                         gatTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTcc                               <  2:3590301/93‑1 (MQ=255)
                                                         gatTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTcc                               >  1:1878719/1‑93 (MQ=255)
                                                         gatTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTcc                               >  2:659548/1‑93 (MQ=255)
                                                          atTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTccc                              >  2:1165903/1‑92 (MQ=255)
                                                          atTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCg                              >  2:889709/1‑93 (MQ=255)
                                                          atTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCg                              >  1:236752/1‑93 (MQ=255)
                                                           tttGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGc                             <  2:63290/93‑1 (MQ=255)
                                                           tttGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGc                             >  2:4510857/1‑93 (MQ=255)
                                                            ttGGTTTCCTCCGTGCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCa                            <  2:848102/93‑1 (MQ=255)
                                                               gTTTCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTg                         >  2:3728693/1‑93 (MQ=255)
                                                                tttCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTgg                        <  2:723863/93‑1 (MQ=255)
                                                                tttCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTgg                        <  2:2347665/93‑1 (MQ=255)
                                                                tttCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGACGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTgg                        <  2:806698/93‑1 (MQ=255)
                                                                 ttCCTCCGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGt                       >  1:3360142/1‑93 (MQ=255)
                                                                  tcctccGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTg                      <  1:1096877/93‑1 (MQ=255)
                                                                   cctccGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGa                     <  2:1354538/93‑1 (MQ=255)
                                                                   cctccGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGa                     >  2:2402636/1‑93 (MQ=255)
                                                                   cctccGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGa                     >  1:479748/1‑93 (MQ=255)
                                                                     tccGTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATa                   >  1:3884297/1‑93 (MQ=255)
                                                                      ccGTTCCGCAGACAGCTCCCCCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATAc                  <  2:2569852/93‑1 (MQ=255)
                                                                       cGTTCCGCAAACAGCTCCCACCTGGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATAcc                 <  2:2599796/93‑1 (MQ=255)
                                                                        gTTCCGCAGACAGTTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCa                >  1:4210571/1‑93 (MQ=255)
                                                                        gTTCCGCAGACAGCTCCCACCTTGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCa                >  2:36150/1‑93 (MQ=255)
                                                                        gTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCa                <  2:814506/93‑1 (MQ=255)
                                                                        gTTCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCa                >  2:58997/1‑93 (MQ=255)
                                                                         ttCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCAt               <  1:2730473/93‑1 (MQ=255)
                                                                          tCCGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCAtc              <  2:150374/93‑1 (MQ=255)
                                                                            cGCAGACAGCTCCCCCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCAtctc            <  2:1215666/93‑1 (MQ=255)
                                                                            cGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATTCCAtctc            >  2:4120879/1‑93 (MQ=255)
                                                                            cGCAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCAtctc            <  1:102768/93‑1 (MQ=255)
                                                                             gCAGACAGCTCCCCCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATctct           <  1:1253235/93‑1 (MQ=255)
                                                                             gCAGACAGCTCCCACCTCGCGGGTCCTGATGACATATACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATctct           <  2:4372432/93‑1 (MQ=255)
                                                                              cAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTg          <  1:1720447/93‑1 (MQ=255)
                                                                              cAGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTg          >  1:4477702/1‑93 (MQ=255)
                                                                               aGACAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTgg         <  2:4011367/93‑1 (MQ=255)
                                                                                  cAGCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTaa      >  1:281011/1‑93 (MQ=255)
                                                                                   aGCGCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAg     <  1:2223409/93‑1 (MQ=255)
                                                                                    gCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGa    >  2:533217/1‑93 (MQ=255)
                                                                                    gCTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGa    >  2:4235087/1‑93 (MQ=255)
                                                                                     cTCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGAt   >  2:1187244/1‑93 (MQ=255)
                                                                                      tCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGAtt  >  2:2081604/1‑93 (MQ=255)
                                                                                         |                                                                                         
ACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATT  >  NC_000913/3964156‑3964334

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