Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,992,052 T→A 89.3% D300E (GAT→GAA cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,992,0520TA89.3% 280.1 / 14.8 107D300E (GAT→GAAcyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (6/6);  new base A (55/40);  total (61/46)
Fisher's exact test for biased strand distribution p-value = 7.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

TACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTGGCGTC  >  NC_000913/3991906‑3992194
                                                                                                                                                  |                                                                                                                                              
  cAAGAGTATCGATTCCCCATACAACGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCCCGTCTGCTGGAGAAAGAAATCAAACAGCGTTTGCCCGACGGCGAGATTGTATCGTTTGGTATCGAACCa                                                                                                                                             >  6:78464/1‑148 (MQ=255)
  cAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACACACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCa                                                                                                                                             >  3:440/1‑148 (MQ=255)
   aaGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATa                                                                                                                                           <  3:317342/149‑1 (MQ=255)
     gagTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACt                                                                                                                                         <  4:45242/149‑1 (MQ=255)
        tATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCa                                                                                                                                      >  5:102618/1‑149 (MQ=255)
          tCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCat                                                                                                                                     >  7:35161/1‑148 (MQ=255)
                 cccATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTgg                                                                                                                             >  1:259125/1‑149 (MQ=255)
                 cccATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTg                                                                                                                              >  1:245888/1‑148 (MQ=255)
                 cccATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTg                                                                                                                              >  5:75881/1‑148 (MQ=255)
                  ccATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGa                                                                                                                            <  2:166404/149‑1 (MQ=255)
                    aTACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGc                                                                                                                          <  8:147601/149‑1 (MQ=255)
                      aCAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCAAACTGCATGATGCTGGAGCgt                                                                                                                        >  6:305173/1‑149 (MQ=255)
                       cAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGGCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCgt                                                                                                                        <  5:23387/148‑1 (MQ=255)
                        aaaGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCgtgt                                                                                                                      <  7:61840/149‑1 (MQ=255)
                        aaaGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGAAATCAAACAGCGTTAGCACGACGGCGAGATTGTATCGGTTGGTCTCGATCCAAACTGCATGATGCAGGAGCgtgt                                                                                                                      >  3:29578/1‑149 (MQ=255)
                            cGGTACTGAAAACACTGCTGCTGGACGCCTCTTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGCTTGTATCGTTTCGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACt                                                                                                                  <  6:297605/149‑1 (MQ=255)
                            cGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACt                                                                                                                  <  8:79276/149‑1 (MQ=255)
                              gTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGa                                                                                                                <  6:102618/149‑1 (MQ=255)
                               tACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGaa                                                                                                               <  8:40258/149‑1 (MQ=255)
                                    aaaaCACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATAc                                                                                                            <  3:291939/147‑1 (MQ=255)
                                    aaaaCACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCt                                                                                                          >  2:271503/1‑149 (MQ=255)
                                    aaaaCACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCt                                                                                                          >  6:108762/1‑149 (MQ=255)
                                    aaaaCACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTGGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATAc                                                                                                            <  5:258220/147‑1 (MQ=255)
                                     aaaCACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCt                                                                                                          >  8:235705/1‑148 (MQ=255)
                                         acTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACgg                                                                                                     >  5:235615/1‑149 (MQ=255)
                                             ctgctgGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGAt                                                                                                 >  3:332163/1‑149 (MQ=255)
                                             ctgctgGAAGCCTATTCCTGGGAATACCCGAACCCACGGCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGa                                                                                                  <  8:154539/148‑1 (MQ=255)
                                                ctgGCAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATCTCAAACAGCGTTTCCACGACGGCGAGATGGTCTCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGAtt                                                                                                <  2:55205/147‑1 (MQ=255)
                                                ctgGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGAtt                                                                                                <  1:271503/147‑1 (MQ=255)
                                                  ggAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAg                                                                                            <  4:127311/149‑1 (MQ=255)
                                                  ggAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCAAGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTg                                                                                               >  7:192547/1‑146 (MQ=255)
                                                   gAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTg                                                                                               >  7:370405/1‑145 (MQ=255)
                                                   gAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGa                                                                                           <  1:191444/149‑1 (MQ=255)
                                                    aaGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGAt                                                                                          >  3:181918/1‑149 (MQ=255)
                                                     aGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACTGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATTATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGAtt                                                                                         <  7:360450/149‑1 (MQ=255)
                                                      gCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACAACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTAATGAATACCTGACGGCGATTGAAGAttt                                                                                        >  6:107333/1‑149 (MQ=255)
                                                         tATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTAc                                                                                     >  6:326722/1‑149 (MQ=255)
                                                          aTTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTa                                                                                      <  6:182804/147‑1 (MQ=255)
                                                          aTTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTa                                                                                      >  5:182804/1‑147 (MQ=255)
                                                           ttCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTAccc                                                                                   <  8:35161/149‑1 (MQ=255)
                                                              cTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGCCGGCGAGATTGTATCGTTTGGTCTCGAACCATACGGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTc                                                                                >  2:295268/1‑149 (MQ=255)
                                                              cTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTa                                                                                      <  5:190767/143‑1 (MQ=255)
                                                              cTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTa                                                                                      >  6:190767/1‑143 (MQ=255)
                                                               tGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGtt                                                                                                                     <  1:255169/111‑1 (MQ=255)
                                                               tGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCt                                                                               >  7:206732/1‑149 (MQ=255)
                                                                gggAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACAGc                                                                                                                                       >  8:195170/1‑92 (MQ=255)
                                                                 ggAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCt                                                                               <  2:245888/147‑1 (MQ=255)
                                                                    aTACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGAtt                                                                          <  8:370405/149‑1 (MQ=255)
                                                                        ccGAACCGACGTCTGCTGGCGAAAGATATCAAACAGCGATTCAACGACGGCGACATTGTAACGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGAGATACTGAATACCTGAAGGCGAATGAAGATTTTACCCGTCTGGATTTAg                                                                       >  2:293676/1‑148 (MQ=255)
                                                                        ccGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGAGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTac                                                                       >  2:216709/1‑147 (MQ=255)
                                                                             cccACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTc                                                                 >  5:78566/1‑149 (MQ=255)
                                                                              ccACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCg                                                                >  6:328858/1‑149 (MQ=255)
                                                                              ccACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCg                                                                >  5:378228/1‑149 (MQ=255)
                                                                                aCGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCgct                                                              >  2:257729/1‑149 (MQ=255)
                                                                                aCGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCgct                                                              >  6:85144/1‑149 (MQ=255)
                                                                                     tgctgGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTc                                                         >  7:77085/1‑149 (MQ=255)
                                                                                       ctgGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTa                                                       <  8:77085/149‑1 (MQ=255)
                                                                                         ggCGAAAGATATTAAACAGAGTTTGGACGACGGCGTGTTTGTTTCGTTTGGTCTCGATCCATACTGCTTGATGCTTGAGCGTGTTACTGAATACCTGACGGCGATTGAAGAGTTTACCCGTGTGGATTTAGTTCGTCGCTGCTTCTAtt                                                     <  4:29578/149‑1 (MQ=255)
                                                                                          gCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTAt                                                      <  7:165971/147‑1 (MQ=255)
                                                                                          gCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGATTGGGCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTAACCGTCTGGAGTTAGTCCGTCGCTGCTTCTAtt                                                     >  6:383393/1‑148 (MQ=255)
                                                                                            gAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTa                                                   >  5:381132/1‑148 (MQ=255)
                                                                                            gAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTa                                                   >  7:170116/1‑148 (MQ=255)
                                                                                            gAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTa                                                   >  5:260994/1‑148 (MQ=255)
                                                                                             aaaGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTGACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTaa                                                  >  5:280943/1‑148 (MQ=255)
                                                                                             aaaGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTaa                                                  >  3:270284/1‑148 (MQ=255)
                                                                                              aaGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTaaaa                                                <  4:270284/149‑1 (MQ=255)
                                                                                                    tCAAACAGCGTTTGCACGACGGCGAGAGTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGGGTTACTGAAAACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGc                                          <  8:116779/149‑1 (MQ=255)
                                                                                                        aCAGCGGTGGCACGACGGCGAGATTGTATCGGTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGACTACCTGACGGCGATTGACGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGaaa                                      <  2:259125/149‑1 (MQ=255)
                                                                                                          aGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAg                                    <  3:31917/149‑1 (MQ=255)
                                                                                                            cGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCt                                  <  6:378228/149‑1 (MQ=255)
                                                                                                                 gCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGcc                             >  2:141173/1‑149 (MQ=255)
                                                                                                                   acgtcgGCGAGATTGTATCGTTTGGTCTCGAAGCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGTCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGt                           <  5:233374/149‑1 (MQ=255)
                                                                                                                   acgacgGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGt                           <  4:181918/149‑1 (MQ=255)
                                                                                                                   acgacgGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGt                           >  6:97184/1‑149 (MQ=255)
                                                                                                                      acgGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGaa                        <  6:235615/149‑1 (MQ=255)
                                                                                                                       cgGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAAc                       <  6:75881/149‑1 (MQ=255)
                                                                                                                        ggCGCGATTGTATCGTTTGGTCTCGAACCATACTGCCAGATGCTGGAGCGTGGTACAGAATACCTGCCGGCGATTGAAGCATTTACCCGTCTGGATTTAGTACGGCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAAcg                      <  2:324296/149‑1 (MQ=255)
                                                                                                                        ggCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAAcg                      >  1:92026/1‑149 (MQ=255)
                                                                                                                        ggCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAAcg                      >  7:169433/1‑149 (MQ=255)
                                                                                                                        ggCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAAcg                      >  2:328678/1‑149 (MQ=255)
                                                                                                                           gagaTTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGaaaa                                     >  1:332851/1‑131 (MQ=255)
                                                                                                                           gagaTTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGaaaa                                     <  2:332851/131‑1 (MQ=255)
                                                                                                                           gagaTTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACgcgc                   >  7:12734/1‑149 (MQ=255)
                                                                                                                           gagaTTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACgcgc                   >  5:249397/1‑149 (MQ=255)
                                                                                                                           gagaTTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACgcgc                   <  5:49235/149‑1 (MQ=255)
                                                                                                                             gaTTGTATCGTTTGGTCTTGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCt                 >  8:150745/1‑149 (MQ=255)
                                                                                                                               ttGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGc               >  2:220877/1‑149 (MQ=255)
                                                                                                                                 gTATCGTGTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTTAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCg              <  5:383393/148‑1 (MQ=255)
                                                                                                                                  tATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCg              >  4:943/1‑147 (MQ=255)
                                                                                                                                   aTCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCAGCCTac           >  3:208589/1‑148 (MQ=255)
                                                                                                                                     cGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGc         >  1:396582/1‑149 (MQ=255)
                                                                                                                                       tttGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTg       >  8:120894/1‑149 (MQ=255)
                                                                                                                                         tGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTGGc     >  7:179419/1‑149 (MQ=255)
                                                                                                                                         tGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACATGACGGCGATTGAAGAtc                                                                                         >  4:49235/1‑64 (MQ=255)
                                                                                                                                          ggTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTGGc     <  6:381132/148‑1 (MQ=255)
                                                                                                                                          ggTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTGGc     <  3:943/148‑1 (MQ=255)
                                                                                                                                          ggTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTGGc     <  7:235705/148‑1 (MQ=255)
                                                                                                                                           gTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGGCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTGGCg    <  3:29646/148‑1 (MQ=255)
                                                                                                                                            tctcGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTGGCGTc  >  3:37283/1‑149 (MQ=255)
                                                                                                                                            tctcGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTGGCGTc  >  7:183903/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
TACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTGGCGTC  >  NC_000913/3991906‑3992194

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