Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,186,578 A→G 89.3% D410G (GAT→GGT)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,186,5780AG89.3% 259.1 / 17.8 84D410G (GAT→GGT) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (6/3);  new base G (38/37);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

ATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCG  >  NC_000913/4186440‑4186714
                                                                                                                                          |                                                                                                                                        
aTCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGat                                                                                                                                         >  1:1361210/1‑140 (MQ=255)
    gTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATcc                                                                                                                                     >  2:970267/1‑140 (MQ=255)
    gTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATcc                                                                                                                                     >  2:970089/1‑140 (MQ=255)
    gTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATcc                                                                                                                                     >  2:1426988/1‑140 (MQ=255)
    gTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCGGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATcc                                                                                                                                     >  1:671896/1‑140 (MQ=255)
        ggTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATTGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGa                                                                                                                                 <  2:132703/140‑1 (MQ=255)
                       aTGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGa                                                                                                                  <  1:1442783/140‑1 (MQ=255)
                       aTGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGa                                                                                                                  <  1:1444349/140‑1 (MQ=255)
                          gCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAaca                                                                                                               <  1:1143515/140‑1 (MQ=255)
                          gCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAaca                                                                                                               <  2:1528867/140‑1 (MQ=255)
                            aCTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACAcc                                                                                                             <  2:567995/140‑1 (MQ=255)
                            aCTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACAcc                                                                                                             <  1:992987/140‑1 (MQ=255)
                              tGGAGCTGTTCAACCCGTTCGTCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCGGACGAAGTTATCCGCGAACACCCgg                                                                                                          <  1:19966/140‑1 (MQ=255)
                               ggcgcagTTCACACCGTACCTCCCCGGCACGCTGGAACAGCGTGGTCTTCCTACCACCATTAAAGCTGCGAAGAAAAGGGTTGAGCGCGAAGAAGCTGTCGATTGGGGTATCCCGGACGAAGTTATCCGCGAACACCCgg                                                                                                          <  2:367328/134‑1 (MQ=255)
                               ggAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCgg                                                                                                          <  1:1305059/140‑1 (MQ=255)
                               ggAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCgg                                                                                                          <  1:1510740/140‑1 (MQ=255)
                                      ttCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTActgct                                                                                                   <  1:166337/140‑1 (MQ=255)
                                          aaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAAc                                                                                               <  1:277720/140‑1 (MQ=255)
                                          aaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAAc                                                                                               >  1:98875/1‑140 (MQ=255)
                                          aaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAAc                                                                                               >  1:1111992/1‑140 (MQ=255)
                                          aaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAAc                                                                                               >  1:24968/1‑140 (MQ=255)
                                           aCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAAcc                                                                                              >  1:1129453/1‑140 (MQ=255)
                                           aCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAAcc                                                                                              >  1:1526740/1‑140 (MQ=255)
                                              gTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTg                                                                                           >  2:28106/1‑140 (MQ=255)
                                                tCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGaaaa                                                                                                               >  1:1265373/1‑116 (MQ=255)
                                                tCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAaca                                                                                                               <  2:1265373/118‑1 (MQ=255)
                                                tCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCa                                                                                         >  1:1627471/1‑140 (MQ=255)
                                                 caacTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCAc                                                                                        <  1:1170953/137‑1 (MQ=255)
                                                 cATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAACGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCAc                                                                                        <  1:966357/140‑1 (MQ=255)
                                                 cATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCAc                                                                                        <  1:1551738/140‑1 (MQ=255)
                                                 cATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCAc                                                                                        <  1:516445/140‑1 (MQ=255)
                                                 cATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCAc                                                                                        <  2:154120/140‑1 (MQ=255)
                                                 cATCTACCGCAAGCTGGAACTGCGTTGTCTAGCTACCACCAATAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCAGGACTAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCAc                                                                                        <  2:1581719/140‑1 (MQ=255)
                                                  aTCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCAcc                                                                                       <  1:1608166/140‑1 (MQ=255)
                                                     tACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGAc                                                                                    <  1:379781/140‑1 (MQ=255)
                                                       cGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGActc                                                                                  >  1:281182/1‑140 (MQ=255)
                                                       cGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGActc                                                                                  >  1:281195/1‑140 (MQ=255)
                                                         gCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTg                                                                                >  2:942197/1‑140 (MQ=255)
                                                         gCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTg                                                                                >  2:942263/1‑140 (MQ=255)
                                                            aGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCAc                                                                             <  1:494563/140‑1 (MQ=255)
                                                             gCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCAcc                                                                            >  1:1363807/1‑140 (MQ=255)
                                                             gCTGGAACTGAGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCAcc                                                                            >  2:1312361/1‑140 (MQ=255)
                                                              cTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCg                                                                           <  1:629691/140‑1 (MQ=255)
                                                                 gAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCt                                                                        <  1:271011/140‑1 (MQ=255)
                                                                       cGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTTGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTAt                                                                  >  1:6162/1‑140 (MQ=255)
                                                                       cGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTAt                                                                  >  1:1591263/1‑140 (MQ=255)
                                                                       cGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTAt                                                                  >  1:598264/1‑140 (MQ=255)
                                                                       cGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTAt                                                                  >  1:1591758/1‑140 (MQ=255)
                                                                           gTCTTGCTACCACCATTAAAGCTGCGAAGAAACTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAg                                                              <  2:1316300/140‑1 (MQ=255)
                                                                                gCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAtt                                                         <  1:494315/140‑1 (MQ=255)
                                                                                 cTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAttt                                                        >  1:810397/1‑140 (MQ=255)
                                                                                    ccaccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGaa                                                     <  1:1378715/140‑1 (MQ=255)
                                                                                     caccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAAc                                                    <  2:658133/140‑1 (MQ=255)
                                                                                     caccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAAc                                                    >  1:70566/1‑140 (MQ=255)
                                                                                      accaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAAcc                                                   >  2:517690/1‑140 (MQ=255)
                                                                                          ttAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTa                                               <  2:661079/140‑1 (MQ=255)
                                                                                          tgaaaGCTGCGAAGAAAATGTGTGAGCGCGAAGAAGCTGTCGTAAGGGATATCCTGGACGAAGTAATCCGCGAACACCCGGACCTCCTCAACCGGGCACCGCCTCTGCACCGTCTGGGTATCCCGGCCTGGGAACCgggc                                               >  2:1459090/3‑138 (MQ=255)
                                                                                           tAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTAc                                              <  2:1479411/140‑1 (MQ=255)
                                                                                           tAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTAc                                              <  2:1479400/140‑1 (MQ=255)
                                                                                                cTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATc                                         <  1:629922/140‑1 (MQ=255)
                                                                                                  gCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTGGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGa                                       >  2:10258/1‑140 (MQ=255)
                                                                                                  gCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTTCACCGTCTGGGTATCCAGGCAATTTAACCGGTACTGATAGa                                       >  1:294937/1‑140 (MQ=255)
                                                                                                  gCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGa                                       >  2:232033/1‑140 (MQ=255)
                                                                                                  gCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGa                                       >  2:232169/1‑140 (MQ=255)
                                                                                                  gCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGa                                       >  2:1031545/1‑140 (MQ=255)
                                                                                                  gCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGa                                       >  2:1225910/1‑140 (MQ=255)
                                                                                                  gCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGa                                       >  1:489854/1‑140 (MQ=255)
                                                                                                  gCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGa                                       >  1:518837/1‑140 (MQ=255)
                                                                                                  gCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGa                                       >  1:745167/1‑140 (MQ=255)
                                                                                                  gCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGa                                       >  2:1232209/1‑140 (MQ=255)
                                                                                                  gCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGa                                       >  1:273845/1‑140 (MQ=255)
                                                                                                  gCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTGTGGGGTATCCTGGACGAAGTTATCCGCGAACCCCCGGTACAGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGa                                       >  1:978245/1‑140 (MQ=255)
                                                                                                   cGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGaa                                      >  1:1733459/1‑140 (MQ=255)
                                                                                                     aagaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAgg                                    >  1:1514402/1‑140 (MQ=255)
                                                                                                     aagaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAgg                                    <  2:1624568/140‑1 (MQ=255)
                                                                                                     aagaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCAGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGAGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACAGATAGAAgg                                    >  1:1608687/1‑140 (MQ=255)
                                                                                                        aaaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTaa                                 <  2:199950/140‑1 (MQ=255)
                                                                                                               ttGAGCGCGAAGACGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATc                          <  1:1668472/140‑1 (MQ=255)
                                                                                                                  aGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAg                       >  1:871060/1‑140 (MQ=255)
                                                                                                                    cgcgAAGAACATGTCGATTGGGATAACCTGGACGAAGTTATCCTCGAACACCCTCTACTGCTGAACCGTGCACCGACTCTGCACCGTCTTGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCt                     <  2:1442758/140‑1 (MQ=255)
                                                                                                                           aagCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCg              <  2:1672335/140‑1 (MQ=255)
                                                                                                                                gTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGt         <  1:231615/140‑1 (MQ=255)
                                                                                                                                 tCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGtt        >  2:287583/1‑140 (MQ=255)
                                                                                                                                 tCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGtt        >  2:389502/1‑140 (MQ=255)
                                                                                                                                     ttGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTtgtg    <  1:1453546/140‑1 (MQ=255)
                                                                                                                                       ggggTATCCTGGACGCAGTAATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCg  <  2:1204035/140‑1 (MQ=255)
                                                                                                                                       ggggTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCg  <  2:1028697/140‑1 (MQ=255)
                                                                                                                                       ggggTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCg  <  2:1204067/140‑1 (MQ=255)
                                                                                                                                       ggggTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCg  <  1:1188606/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                        
ATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCG  >  NC_000913/4186440‑4186714

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