Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,638,818 G→T 6.8% A393S (GCG→TCG)  pitA → metal phosphate:H(+) symporter PitA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,638,8180GT 6.8% 412.4 / 4.5 128A393S (GCG→TCG) pitAmetal phosphate:H(+) symporter PitA
Reads supporting (aligned to +/- strand):  ref base G (50/67);  new base T (9/2);  total (59/69)
Fisher's exact test for biased strand distribution p-value = 2.30e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

ATGCTGTGCGTTTCTGACACTATCGACAAAGTGGTGAAGATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGA  >  NC_000913/3638671‑3638965
                                                                                                                                                   |                                                                                                                                                   
aTGCTGTGCGTTTCTGACACTATCGACAAAGTGGTGAAGATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCg                                                                                                                                                   <  2:580947/150‑1 (MQ=255)
 tGCTGTCCGTTTCTGACACTATCGACAAAGTGGTGAAGATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCTCTGGCGt                                                                                                                                                  >  2:1343408/1‑150 (MQ=255)
  gCTGTGCGTTTCTGACACTATCGACAAAGTGGTGAAGATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGtt                                                                                                                                                 <  1:518661/150‑1 (MQ=255)
   cTGTGCGTTTCTGACACTATCGACAAAGTGGTGAAGATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTa                                                                                                                                                >  2:1628486/1‑150 (MQ=255)
     gtgCGTTTCTGACACTATCGACAAAGTGGTGAAGATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAgg                                                                                                                                              <  1:101375/150‑1 (MQ=255)
       gCGTTTCTGACACTATCGACAAAGTGGTGAAGATGCCTGGCGTGAGTGCTGACGATCAGCGCCGGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTa                                                                                                                                            <  1:1652542/150‑1 (MQ=255)
         gTTTCTGACACTATCGACAAAGTGGTGAAGATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATc                                                                                                                                          <  2:1477000/150‑1 (MQ=255)
           ttCTGACACTATCGACAAAGTGGTGAAGATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCgg                                                                                                                                        <  1:1361480/150‑1 (MQ=255)
           ttCTGACACTATCGACAAAGTGGTGAAGATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCgg                                                                                                                                        <  2:29630/150‑1 (MQ=255)
           ttCTGACACTATCGACAAAGTGGTGAAGATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCgg                                                                                                                                        <  1:1932040/150‑1 (MQ=255)
               gACACTATCGACAAAGTGGTGAAGATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACg                                                                                                                                    >  1:493068/1‑150 (MQ=255)
               gACACTATCGACAAAGTGGTGAAGATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACg                                                                                                                                    <  1:1674227/150‑1 (MQ=255)
                   cTATCGACAAAGTGGTGAAGATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTAGATCAACATGGAGGACGAACTCTAGTTAGGTATAGGTACgatga                                                                                                                                >  2:2420218/1‑150 (MQ=255)
                    tATCGACAAAGTGGTGAAGATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACgatgat                                                                                                                               >  1:1945731/1‑150 (MQ=255)
                      tCGTCATAGTGGTGAAGATGCCTGGCGTGAGTGCTGACTTTCAGCGCCTGTTGAAGAAACTTAAGTCCGACATGCTTAGCACCATCGTGTATGCACCGGTGTGTATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGTTtg                                                                                                                             <  1:882438/150‑1 (MQ=255)
                       cGACAAAGTGGTGAAGATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATtgg                                                                                                                            >  1:2429804/1‑150 (MQ=255)
                        gACAAAGTGGTGAAGATGCCTGGCGTCAGTGCTGACGATCAGCGGCTGTTGAAGCAACTGAAGTCCGACATGCTTCGCCCCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGGTtggg                                                                                                                           >  2:519635/1‑149 (MQ=255)
                          cAAAGTGGTGAAGATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATtggctg                                                                                                                         <  1:1687643/150‑1 (MQ=255)
                           aaaGTGGTGAAGATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATtggctgg                                                                                                                        <  1:614217/150‑1 (MQ=255)
                                ggtgAAGATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCg                                                                                                                   >  2:2066238/1‑150 (MQ=255)
                                 gtgAAGATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCgt                                                                                                                  <  2:117417/150‑1 (MQ=255)
                                     aGATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCTGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTgg                                                                                                              <  1:2101236/150‑1 (MQ=255)
                                     aGATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTgg                                                                                                              <  2:128290/150‑1 (MQ=255)
                                      gATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGc                                                                                                             <  1:820321/150‑1 (MQ=255)
                                      gATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGc                                                                                                             <  1:14139/150‑1 (MQ=255)
                                        tGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCaa                                                                                                           >  2:1391526/1‑150 (MQ=255)
                                         gCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAAc                                                                                                          <  2:93198/150‑1 (MQ=255)
                                          ccTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACg                                                                                                         <  2:1793600/150‑1 (MQ=255)
                                          ccTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACg                                                                                                         <  1:1879791/150‑1 (MQ=255)
                                            tGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGAc                                                                                                       >  2:221804/1‑150 (MQ=255)
                                            tGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGAc                                                                                                       <  1:311262/150‑1 (MQ=255)
                                                   aGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACg                                                                                                         <  1:109647/141‑1 (MQ=255)
                                                    gTGCTGAGGATCAGCGCCTGTTGAAGCAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATCGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTg                                                                                               >  1:325678/1‑150 (MQ=255)
                                                    gTGCTGACGATCAGCGCCTGTTGATGAAACTGAAGTCCGACATGCTTAGCACGATCGTGTATGCACCGGTGTGGATCATCATTGCGGTCGCGCTGTCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGAGTATCGGTg                                                                                               <  2:2182432/150‑1 (MQ=255)
                                                    gTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTg                                                                                               <  2:2234724/150‑1 (MQ=255)
                                                      gCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTgag                                                                                             >  2:738784/1‑150 (MQ=255)
                                                       cTGCCGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTgaga                                                                                            >  1:670082/1‑150 (MQ=255)
                                                             aTCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCg                                                                                      <  1:2404381/150‑1 (MQ=255)
                                                              tCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCgg                                                                                     <  2:2229096/150‑1 (MQ=255)
                                                               cAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGt                                                                                    <  1:2421509/150‑1 (MQ=255)
                                                               cAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGt                                                                                    <  1:2457100/150‑1 (MQ=255)
                                                               cAGCGCCTGTTGAAGAAACTGAAGTCAGACATGCTTAGCACCATCGAGTATGCACCGGTGTAGATCAACATGGACGTCGAGATGTAAAAAAATAAAGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGt                                                                                    >  2:763213/1‑150 (MQ=255)
                                                                aGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTa                                                                                   >  1:1655610/1‑150 (MQ=255)
                                                                aGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTa                                                                                   <  1:191511/150‑1 (MQ=255)
                                                                aGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTAGATCAACATGGACGACGAAAAATAAATAAATAACGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTa                                                                                   >  2:1246002/1‑150 (MQ=255)
                                                                  cgcCTGTTGAAGAAACTGAACTCCGACATGCTTAGCACCATCGAGGATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGTCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAg                                                                                 <  2:1862932/150‑1 (MQ=255)
                                                                      tGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTAGATCAACATGGACGACGAAAAATAAAAAAAAAAAGGTACGATGATTGGCTGGCGCCGTGTGGCAAAGACTATCGGTGAGAAAATCGGTAAGAAAg                                                                             >  2:1825980/1‑150 (MQ=255)
                                                                      tGTTCAACAAAATGAAATCCAACATACTTAGCACAATCTAGAATTCACCGGAGTATATCTAAATAAAAATCGAACAAGAGTTAAGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGGGAGAAAATCGGTAAGAAAg                                                                             >  2:1626400/1‑150 (MQ=255)
                                                                       gTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTTTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAgg                                                                            >  1:1093141/1‑150 (MQ=255)
                                                                         tGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCa                                                                          >  1:2008453/1‑150 (MQ=255)
                                                                          gaagaaACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCgg                                                                                     >  1:492725/1‑138 (MQ=255)
                                                                          gaagaaACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCgg                                                                                     <  2:492725/138‑1 (MQ=255)
                                                                            agaaACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGa                                                                       <  1:2434820/150‑1 (MQ=255)
                                                                             gaaACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGAc                                                                      <  2:2429804/150‑1 (MQ=255)
                                                                              aaaCTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGAcc                                                                     <  2:687553/150‑1 (MQ=255)
                                                                               aaCTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGAGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCt                                                                    >  2:117749/1‑150 (MQ=255)
                                                                                aCTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTa                                                                   <  2:836000/150‑1 (MQ=255)
                                                                                  tGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACg                                                                 <  1:1391526/150‑1 (MQ=255)
                                                                                     aGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTc                                                              <  1:695863/150‑1 (MQ=255)
                                                                                        ccGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGCCTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAgg                                                           >  1:1820266/1‑150 (MQ=255)
                                                                                        ccGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAgg                                                           >  1:1263792/1‑150 (MQ=255)
                                                                                        ccGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAgg                                                           >  2:707294/1‑150 (MQ=255)
                                                                                         cGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAggg                                                          >  1:1158390/1‑150 (MQ=255)
                                                                                         cGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAggg                                                          <  1:221804/150‑1 (MQ=255)
                                                                                         cGACATGCTTAGCACCATCGAGTATGCACCGGTGTAGATCAACATGGACGACGAAAAATAAAAAAAAAAAAGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTATGAAAGCCATGACATACGGGCAggg                                                          >  2:1750790/1‑150 (MQ=255)
                                                                                          gACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCACCGGCTATCGGTGAGAAAATCGGTAAGAAAGGCATGCCCTACGCTCAgggg                                                         >  1:2282488/1‑150 (MQ=255)
                                                                                               gCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTc                                                    <  1:1708372/150‑1 (MQ=255)
                                                                                               gCTTAGCACCATCGAGTATGCACCGGTGTAGATCAACATGGACGACGAAAAATAAAAAAAAAAAGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGTATGACCTACGCTCAGGGGATGTc                                                    >  2:1721985/1‑150 (MQ=255)
                                                                                               gCTTAGCACCATCGAGTATGCACCGGTGTAGATCAACATGGACGACGAAAAATAAAAAAAAAAAAGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTc                                                    >  2:360174/1‑150 (MQ=255)
                                                                                                  tAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGc                                                 <  1:2187240/150‑1 (MQ=255)
                                                                                                     cACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCa                                              <  2:1659908/150‑1 (MQ=255)
                                                                                                     cACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCa                                              <  2:1480045/150‑1 (MQ=255)
                                                                                                     cACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCa                                              <  2:493068/150‑1 (MQ=255)
                                                                                                      aaaaTCGAGTATGCACCGGTGTGGATCATCATTGAAGTAGATCAAGCTATAAGTATATGTACTATGATTGTCTTGCTCCTTGTGGCAAAGACTATAGGTTAGAAAATCGGTAAGAAAGTAATGAACTACTATCAGGTGATTTCTGCACAg                                             >  2:1846178/4‑150 (MQ=255)
                                                                                                       ccATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGa                                            >  2:2148601/1‑150 (MQ=255)
                                                                                                         aTCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATg                                          <  2:2010302/150‑1 (MQ=255)
                                                                                                         aTCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATg                                          >  1:858406/1‑150 (MQ=255)
                                                                                                             aGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGAcgg                                      >  1:446746/1‑150 (MQ=255)
                                                                                                               tATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTAAGCTCAGGGGATGTCTGCCCAGATGAcggcg                                    >  1:1311193/1‑150 (MQ=255)
                                                                                                                 tGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACggcggc                                  <  1:71911/150‑1 (MQ=255)
                                                                                                                 tGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACggcggc                                  >  2:632965/1‑150 (MQ=255)
                                                                                                                 tGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACggcggc                                  <  2:577375/150‑1 (MQ=255)
                                                                                                                 tGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACggcggc                                  >  2:1558184/1‑150 (MQ=255)
                                                                                                                     ccGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCTGGGGATGTCTGCCCAGATGACGGCGGCAgtg                              <  1:1293849/150‑1 (MQ=255)
                                                                                                                     ccGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGGATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAgtg                              <  2:383851/150‑1 (MQ=255)
                                                                                                                       ggTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTc                            >  1:1892090/1‑150 (MQ=255)
                                                                                                                            ggATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCg                       >  1:1584428/1‑150 (MQ=255)
                                                                                                                            ggATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCg                       >  1:240119/1‑150 (MQ=255)
                                                                                                                             gATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCgg                      >  2:1333305/1‑150 (MQ=255)
                                                                                                                               tcatcaTGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGcc                    <  2:2151263/150‑1 (MQ=255)
                                                                                                                                    atGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTTGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCg               >  2:1600744/1‑150 (MQ=255)
                                                                                                                                    atGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCg               >  1:331703/1‑150 (MQ=255)
                                                                                                                                     tGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGa              >  1:2147160/1‑150 (MQ=255)
                                                                                                                                       gCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACTCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGt            <  2:1945731/150‑1 (MQ=255)
                                                                                                                                       gCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGt            <  2:1259194/150‑1 (MQ=255)
                                                                                                                                       gCGGTCGCGCTGGCGTTAGGCATCGGTACGATGACTGGCTGGCGCCGTGTGGCAACGACTATCGGGGAGAAAATCGGTAAGAAAGGCATGACCTACCCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGt            <  2:2295883/150‑1 (MQ=255)
                                                                                                                                         ggTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGAGTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTta          <  1:1909563/150‑1 (MQ=255)
                                                                                                                                                cTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCggg   <  1:1302470/150‑1 (MQ=255)
                                                                                                                                                 tGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGa  <  2:1919024/150‑1 (MQ=255)
                                                                                                                                                 tGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGa  <  2:1158390/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATGCTGTGCGTTTCTGACACTATCGACAAAGTGGTGAAGATGCCTGGCGTGAGTGCTGACGATCAGCGCCTGTTGAAGAAACTGAAGTCCGACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGA  >  NC_000913/3638671‑3638965

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