Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,966,571 A→G 5.2% D52G (GAT→GGT)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,966,5710AG 5.2% 244.6 / 2.7 78D52G (GAT→GGT) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base A (32/41);  new base G (2/2);  total (35/43)
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.91e-01

CTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGAT  >  NC_000913/3966423‑3966718
                                                                                                                                                    |                                                                                                                                                   
cTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGa                                                                                                                                                     <  3:396800/149‑1 (MQ=255)
   aCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATg                                                                                                                                                   <  5:26615/148‑1 (MQ=255)
    ccGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAATTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGc                                                                                                                                                 <  6:378527/149‑1 (MQ=255)
    ccGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGCAGATATCTTTGGTGATGGc                                                                                                                                                 <  3:272581/149‑1 (MQ=255)
    ccGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATgg                                                                                                                                                  <  5:383289/148‑1 (MQ=255)
       aaTTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTa                                                                                                                                              >  1:422515/1‑149 (MQ=255)
         ttAAAGAATACGCCGGGTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACt                                                                                                                                            <  2:89111/149‑1 (MQ=255)
           aaaGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTg                                                                                                                                           >  2:310065/1‑148 (MQ=255)
            aaGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGga                                                                                                                                         <  1:350956/149‑1 (MQ=255)
            aaGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGga                                                                                                                                         <  8:88871/149‑1 (MQ=255)
            aaGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGga                                                                                                                                         <  1:239993/149‑1 (MQ=255)
              gAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGgaga                                                                                                                                       <  6:132182/149‑1 (MQ=255)
               aaTACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGGACTGgaga                                                                                                                                       >  6:341256/1‑148 (MQ=255)
                 tACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGatat                                                                                                                                    >  3:341228/1‑149 (MQ=255)
                   cGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGgag                                                                                                                                        <  5:360571/143‑1 (MQ=255)
                   cGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGgag                                                                                                                                        >  6:360571/1‑143 (MQ=255)
                    gCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTg                                                                                                                                  >  6:379170/1‑148 (MQ=255)
                     ccGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCa                                                                                                                                <  8:35842/149‑1 (MQ=255)
                      cgttttCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCGCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAg                                                                                                                               <  5:53106/146‑1 (MQ=255)
                      cGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGGTGGCGTACTGGAGATATTGCat                                                                                                                               <  3:63232/149‑2 (MQ=255)
                      cGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGGTGGCGTACTGGAGATATTGCat                                                                                                                               >  4:63232/1‑148 (MQ=255)
                             tGAGCTGATCACTCTCGGCGAAAATTTGGGGCTGGAAAAACTGGCTCGTATGCGTAAGCAGGACAGTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGATGATATCGTAGGTGCTGGCGTACTGGCGATATTGCAgggtggct                                                                                                                        >  6:310235/1‑147 (MQ=255)
                               aGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggatgga                                                                                                                         <  5:68022/146‑1 (MQ=255)
                                gCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCTTAACCAGGACATTATTTTTGCCATCATGAAGCAGCACGCAAAGAGTGGCGAAGATATGTTTGGGGATGGGGTAAGGGTGATATTGCTGGATGGATTTg                                                                                                                     >  6:104013/1‑149 (MQ=255)
                                 cTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTgg                                                                                                                    <  3:243710/149‑1 (MQ=255)
                                    ctCACTCTCGTCGAAAATATGGGGCTGGAAAAGCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGt                                                                                                                   <  5:177056/146‑1 (MQ=255)
                                    aTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGtt                                                                                                                  <  1:26480/148‑1 (MQ=255)
                                     tCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTtc                                                                                                                <  5:104013/149‑1 (MQ=255)
                                       aCTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTtcct                                                                                                              <  1:310065/149‑1 (MQ=255)
                                         tctcGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTtcctcc                                                                                                            >  1:286895/1‑149 (MQ=255)
                                           tcGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGt                                                                                                          >  1:420618/1‑149 (MQ=255)
                                               cGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGAt                                                                                                                                                    <  7:315845/103‑1 (MQ=255)
                                               cGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGAt                                                                                                                                                    >  8:315845/1‑103 (MQ=255)
                                               cGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTtcctccagat                                                                                                        >  3:220233/1‑143 (MQ=255)
                                               cGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTtcctc                                                                                                             <  4:220233/142‑1 (MQ=255)
                                                gAAACTATGGGGCTGGAAAACCTGGCTCGTATGCGGAAGCAGGACATTATTTTTGCCAGCCTGAAGCAGCTCGCAAAGAGTGGCGAAGATATCTTTGGTGATGTCGTACTGGAGATATTGCAGGATGTATTTGGTTTACTCCGTTCCGc                                                                                                     <  6:105318/149‑1 (MQ=255)
                                                 aaaaTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCa                                                                                                    <  2:422515/149‑1 (MQ=255)
                                                   aaTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGa                                                                                                  <  2:286895/149‑1 (MQ=255)
                                                    atatGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGAc                                                                                                 <  2:420618/149‑1 (MQ=255)
                                                               gAAAACCTGGCTCGTAAGAGTAAGCAGGACATTATTTTTGCCATACTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGAGATGGAGTACAGGAGATAATGCAGGATTGATATGGTTTCCTCCGTTACGCAGACAGCTCATACCt                                                                                      >  8:298836/1‑149 (MQ=255)
                                                                aaaaCCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGCAGATATCTTTGGTGATGGCGTACTGGAGCTATTGCAGGATGGATTTGGGTTCCTCCGTTCCGCAGCCAGCTCCTACCTc                                                                                     >  6:216128/1‑149 (MQ=255)
                                                                aaaaCCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTc                                                                                     >  8:125718/1‑149 (MQ=255)
                                                                  aaCCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGa                                                                                                  <  2:398256/134‑1 (MQ=255)
                                                                  aaCCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGa                                                                                                  >  1:398256/1‑134 (MQ=255)
                                                                     cTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCgg                                                                                >  2:257976/1‑149 (MQ=255)
                                                                        gCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTc                                                                              <  4:312865/148‑1 (MQ=255)
                                                                        gCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTc                                                                              <  7:298836/148‑1 (MQ=255)
                                                                            gTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgat                                                                         <  6:7727/149‑1 (MQ=255)
                                                                              aTGCGTAAGCAGGACATTATTTTTGCCATCCTGGAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgatga                                                                       >  5:130347/1‑149 (MQ=255)
                                                                               tGCGTAAGCAGGACATTATTTTTGCCATCCTGGAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgatga                                                                       <  6:130347/148‑1 (MQ=255)
                                                                                  gTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTtcctcc                                                                                                            <  1:78495/108‑1 (MQ=255)
                                                                                  gTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTtcctcc                                                                                                            >  2:78495/1‑108 (MQ=255)
                                                                                      gCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACg                                                               >  7:388345/1‑149 (MQ=255)
                                                                                       cAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATTGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGt                                                              >  5:46764/1‑149 (MQ=255)
                                                                                         ggACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGATATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGttt                                                            <  1:257976/149‑1 (MQ=255)
                                                                                           aCATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTcc                                                          <  3:91772/149‑1 (MQ=255)
                                                                                            cATTATTTTTGCCATGCTGAAGCTGCAGGGAAATAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCTCCGGTCCTGATGACATCTACGTTTccc                                                         <  1:138338/149‑1 (MQ=255)
                                                                                             attattTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTcccc                                                        >  4:368009/1‑149 (MQ=255)
                                                                                             attattTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTcccc                                                        >  8:249392/1‑149 (MQ=255)
                                                                                                   tttGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGcc                                                   <  4:341228/148‑1 (MQ=255)
                                                                                                   tttGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGcc                                                   <  7:249392/148‑1 (MQ=255)
                                                                                                            cTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccgccg                                         <  5:379170/149‑1 (MQ=255)
                                                                                                            cTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTCCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTACCCCTAGCCAAAAccgccg                                         >  1:406126/1‑149 (MQ=255)
                                                                                                                     caCGCAAAGAGTGGCGAAGATATCTTTGGTTATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgat                                                                         <  8:318502/108‑1 (MQ=255)
                                                                                                                     caCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCt                                >  2:427466/1‑149 (MQ=255)
                                                                                                                     caCGCAAAGAGTGGCGAAGATATCTTTGGTAATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgat                                                                         >  7:318502/1‑108 (MQ=255)
                                                                                                                      aCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATACGCCGTTTCAACCTc                               >  7:389942/1‑149 (MQ=255)
                                                                                                                       cGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCTAAATCCGCCGTTTCAACCTcc                              <  3:368009/149‑1 (MQ=255)
                                                                                                                        gCAAAGAGTGGCGAAGATATCTTTGGTGGTGGCGTACTGGAGATATTGCAGGATGGATTTGGtt                                                                                                                  <  6:117158/64‑1 (MQ=255)
                                                                                                                        gCAAAGAGGGGCGAAGATATCTTTGGTGGTGGCGTACTGGAGATATTGCAGGATGGATTTGGtt                                                                                                                  >  5:117158/1‑64 (MQ=255)
                                                                                                                          aaaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCa                           <  1:427466/149‑1 (MQ=255)
                                                                                                                          aaaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCa                           >  1:115311/1‑149 (MQ=255)
                                                                                                                          aaaGAGTGGCGAAGATATCTGTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCa                           >  1:63067/1‑149 (MQ=255)
                                                                                                                                  gCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGa                    <  3:102722/148‑1 (MQ=255)
                                                                                                                                         tatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCAtct            >  3:247971/1‑149 (MQ=255)
                                                                                                                                           tCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCATAGCCAAATCCGCCTTTTCAACCTCCGCACTGGTGATACCATctct          >  8:383326/1‑149 (MQ=255)
                                                                                                                                             tttGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTgg        >  3:178328/1‑149 (MQ=255)
                                                                                                                                               tGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATctct          <  2:115311/145‑1 (MQ=255)
                                                                                                                                                   gATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGAt  >  4:398351/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGAT  >  NC_000913/3966423‑3966718

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