Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,966,943 T→G 6.8% I176S (ATT→AGT)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,966,9430TG 6.8% 650.3 / 32.5 193I176S (ATT→AGT) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base T (87/93);  new base G (8/5);  total (95/98)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-01

Only 100 of 194 total aligned reads displayed.

CGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGT  >  NC_000913/3966794‑3967091
                                                                                                                                                     |                                                                                                                                                    
cGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGAt                                                                                                                                                      <  1:2319621/150‑1 (MQ=255)
     gCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTgg                                                                                                                                                 <  2:1853411/150‑1 (MQ=255)
     gCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTgg                                                                                                                                                 <  1:2314234/150‑1 (MQ=255)
        acaaAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAc                                                                                                                                              <  2:2243986/150‑1 (MQ=255)
                ctctTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGaaa                                                                                                                                      >  2:618880/1‑150 (MQ=255)
                   tttGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGTATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGcc                                                                                                                                   >  2:1360918/1‑150 (MQ=255)
                   tttGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGcc                                                                                                                                   >  2:807974/1‑150 (MQ=255)
                   tttGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGcc                                                                                                                                   >  2:319556/1‑150 (MQ=255)
                   tttGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGcc                                                                                                                                   >  2:1109104/1‑150 (MQ=255)
                   tttGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGcc                                                                                                                                   <  1:2700261/150‑1 (MQ=255)
                    ttGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCg                                                                                                                                  <  1:2875078/150‑1 (MQ=255)
                      gagaACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGt                                                                                                                                >  2:1693827/1‑150 (MQ=255)
                          aCTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTaaaa                                                                                                                            >  1:1039709/1‑150 (MQ=255)
                            ttAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAAcc                                                                                                                          >  1:755525/1‑150 (MQ=255)
                            ttAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAAcc                                                                                                                          <  1:2266893/150‑1 (MQ=255)
                                      tGCTCGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGTAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAtgctgctgc                                                                                                                <  2:2398465/150‑1 (MQ=255)
                                       gCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCa                                                                                                               <  2:1829190/150‑1 (MQ=255)
                                           gCAAACTCTCGTCTGCGTATGGAACGTGCTAACGGTTCTACTGAAGATTTAACTGCTCCCGAACAGCAACTCGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAAc                                                                                                           >  2:2566102/1‑150 (MQ=255)
                                            cAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACa                                                                                                          <  1:2038335/150‑1 (MQ=255)
                                            cAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACa                                                                                                          <  2:641567/150‑1 (MQ=255)
                                              aaCTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACAtt                                                                                                        >  1:759454/1‑150 (MQ=255)
                                              aaCTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACAtt                                                                                                        >  1:2400756/1‑150 (MQ=255)
                                                 tctcGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCt                                                                                                     >  1:409148/1‑150 (MQ=255)
                                                   tcGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCa                                                                                                   >  1:326973/1‑150 (MQ=255)
                                                   tcGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCa                                                                                                   >  2:1577531/1‑150 (MQ=255)
                                                      tCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCagag                                                                                                <  1:2417057/150‑1 (MQ=255)
                                                       cTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGc                                                                                               >  2:2241316/1‑150 (MQ=255)
                                                          cGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCAtt                                                                                            >  2:1555997/1‑150 (MQ=255)
                                                          cGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCAtt                                                                                            <  2:1077813/150‑1 (MQ=255)
                                                             aTGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCt                                                                                         <  1:1636698/150‑1 (MQ=255)
                                                              tGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCtt                                                                                        >  1:1679683/1‑150 (MQ=255)
                                                              tGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGAGTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCtt                                                                                        >  2:538364/1‑150 (MQ=255)
                                                              tGGAACGTGGTAACGGTTCTACTGAAGATATAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGAGTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCtt                                                                                        >  1:722184/1‑150 (MQ=255)
                                                                gAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTAc                                                                                      >  2:2370617/1‑150 (MQ=255)
                                                                gAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGACATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTAc                                                                                      >  2:744275/1‑150 (MQ=255)
                                                                    gTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAacc                                                                                  <  1:732389/150‑1 (MQ=255)
                                                                        tAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCAccc                                                                              >  2:1608636/1‑150 (MQ=255)
                                                                            ggTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGAt                                                                          <  1:885279/150‑1 (MQ=255)
                                                                               tCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGAGTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATtgt                                                                       <  2:780956/150‑1 (MQ=255)
                                                                                cttcTGATGATTTATCTGCTCGGGTAGTGGTTCTGGCTTCTCCTTTCGGTCGTGGTCTGCGTGGTCTGATTGTGGCACCTCCGTATGCCGGTTTTTCCTTGGTGCTGCAGATCATTGCTCAGAGCATTGCTTACATCCAGCGGGTTtgtg                                                                      <  2:1783268/147‑1 (MQ=255)
                                                                                cTACTGAATATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCTAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCTTTGCTAACAACGACCCGGATtgtg                                                                      <  1:322286/150‑1 (MQ=255)
                                                                                   cTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGc                                                                   <  2:322009/150‑1 (MQ=255)
                                                                                   cTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGc                                                                   <  1:2723417/150‑1 (MQ=255)
                                                                                      aaGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGa                                                                <  2:1737743/150‑1 (MQ=255)
                                                                                       agaataaaCTGCTCGCGTACTGGATCTGCCATCACCTATCGGTCGTGATCAGCGTGGACTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGAt                                                               <  2:2261831/144‑1 (MQ=255)
                                                                                       aGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGAt                                                               <  2:2509883/150‑1 (MQ=255)
                                                                                       aGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGAt                                                               >  1:2526058/1‑150 (MQ=255)
                                                                                             aaCTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTct                                                         <  1:1338938/150‑1 (MQ=255)
                                                                                                tGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTctgct                                                      >  2:1325573/1‑150 (MQ=255)
                                                                                                   tCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGAt                                                   <  1:1742402/150‑1 (MQ=255)
                                                                                                   tCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGAt                                                   <  1:2670293/150‑1 (MQ=255)
                                                                                                   tCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGAt                                                   <  1:440786/150‑1 (MQ=255)
                                                                                                   tCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGAt                                                   <  2:2757054/150‑1 (MQ=255)
                                                                                                     gcgTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATcg                                                 <  2:2681989/150‑1 (MQ=255)
                                                                                                     gcgTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATcg                                                 <  1:190935/150‑1 (MQ=255)
                                                                                                     gcgTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGAGTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATcg                                                 >  1:985842/1‑150 (MQ=255)
                                                                                                     gcgTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGAGTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATcg                                                 >  1:489312/1‑150 (MQ=255)
                                                                                                      cgTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATcga                                                <  1:215624/150‑1 (MQ=255)
                                                                                                      cgTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATcga                                                <  1:807974/150‑1 (MQ=255)
                                                                                                         aCTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATcgacga                                             <  2:2358810/150‑1 (MQ=255)
                                                                                                           tGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACGACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAAc                                           >  1:1540866/1‑150 (MQ=255)
                                                                                                           tGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAAc                                           <  1:114497/150‑1 (MQ=255)
                                                                                                           tGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAAc                                           <  2:2616134/150‑1 (MQ=255)
                                                                                                           tGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGAGTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAAc                                           >  1:1943428/1‑150 (MQ=255)
                                                                                                               tCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGAGTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTcc                                       >  2:1024133/1‑150 (MQ=255)
                                                                                                                 tGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCgg                                     >  2:691593/1‑150 (MQ=255)
                                                                                                                   gCATCACCTATCGGTCGTGGTCTGGGTGGTGTGTTTGTGGCACCGCCGAAAGCCGGTAAAACCTTGCTGCTGAATAAATTTGCTCTGAGCTTTGATTACAACGACCCGGTTTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGgaa                                   <  1:1235139/150‑1 (MQ=255)
                                                                                                                   gCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGTTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGgaa                                   <  1:629596/150‑1 (MQ=255)
                                                                                                                   gCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGgaa                                   <  1:52744/150‑1 (MQ=255)
                                                                                                                      tCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATTGTTCTGCTGATCGACGAACGTCCGgaagaa                                <  1:2370617/150‑1 (MQ=255)
                                                                                                                          cTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGTTAAAACCTTGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTaa                            <  1:2474479/150‑1 (MQ=255)
                                                                                                                          cTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTaa                            >  1:1262937/1‑150 (MQ=255)
                                                                                                                           tATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAAc                           <  1:1058068/150‑1 (MQ=255)
                                                                                                                           tATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAAc                           >  2:920482/1‑150 (MQ=255)
                                                                                                                           tATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAAc                           >  1:1860426/1‑150 (MQ=255)
                                                                                                                             tCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCg                         >  1:575446/1‑150 (MQ=255)
                                                                                                                              cGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCga                        <  1:2727702/150‑1 (MQ=255)
                                                                                                                              cGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCga                        <  2:1454127/150‑1 (MQ=255)
                                                                                                                                gTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCgaga                      <  1:2127768/150‑1 (MQ=255)
                                                                                                                                gTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCgaga                      >  2:1282899/1‑150 (MQ=255)
                                                                                                                                   gTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGc                   <  2:1674982/150‑1 (MQ=255)
                                                                                                                                   gTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGc                   >  1:1461328/1‑150 (MQ=255)
                                                                                                                                      gTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGc                >  2:2818811/1‑150 (MQ=255)
                                                                                                                                       tCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCg               >  1:1598688/1‑150 (MQ=255)
                                                                                                                                       tCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCg               >  2:589026/1‑150 (MQ=255)
                                                                                                                                         aGCGTGGTCTGAGTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTc             <  1:2642085/150‑1 (MQ=255)
                                                                                                                                           cGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTg           >  2:28361/1‑150 (MQ=255)
                                                                                                                                            gTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTgg          >  1:837295/1‑150 (MQ=255)
                                                                                                                                            gTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTgg          >  1:2616159/1‑150 (MQ=255)
                                                                                                                                             tGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGt         >  1:588767/1‑150 (MQ=255)
                                                                                                                                             tGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGt         <  1:1109104/150‑1 (MQ=255)
                                                                                                                                              ggTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACGACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTa        <  2:227634/150‑1 (MQ=255)
                                                                                                                                                tCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTaaa      >  1:2758139/1‑150 (MQ=255)
                                                                                                                                                 cTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTaaaa     <  1:1005580/150‑1 (MQ=255)
                                                                                                                                                 cTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTaaaa     <  1:2682858/150‑1 (MQ=255)
                                                                                                                                                  tGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAg    <  2:1970839/150‑1 (MQ=255)
                                                                                                                                                  tGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAg    >  2:1955524/1‑150 (MQ=255)
                                                                                                                                                  tGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGAGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAg    >  1:1637953/1‑150 (MQ=255)
                                                                                                                                                   gATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAgg   >  1:1376155/1‑150 (MQ=255)
                                                                                                                                                    aTTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGt  >  2:1816710/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGT  >  NC_000913/3966794‑3967091

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