Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 733,650 A→T 5.1% intergenic (+21/‑26) PP_t09 → / → PP_t10 tRNA‑Lys/tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947733,6500AT 5.1% 291.0 / 3.1 94intergenic (+21/‑26)PP_t09/PP_t10tRNA‑Lys/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base A (48/40);  new base T (3/2);  total (52/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACC  >  NC_002947/733510‑733797
                                                                                                                                             |                                                                                                                                                   
gCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGccc                                                                                                                                            <  2:2122113/150‑1 (MQ=255)
gCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGccc                                                                                                                                            >  1:1441619/1‑150 (MQ=255)
gCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGccc                                                                                                                                            <  1:1294111/150‑1 (MQ=255)
 cAGGTTGTCAGTAGTTCGAATCCCCTCACCCCGACCATTTTTCGGGCGGGTACGTAAGTTCGTCGAGCAGGTGACTTGTATCCAATTGGGTGTAGGTCTCG‑GTGCTCGACGACCCATTATGTAAAAAAGGGCATCTCAAAGAGATGCCCt                                                                                                                                           >  2:35079/1‑150 (MQ=9)
    gATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGCCCtttt                                                                                                                                        >  2:330698/1‑150 (MQ=255)
    gATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCtttt                                                                                                                                        >  1:472956/1‑150 (MQ=255)
    gATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCtttt                                                                                                                                        <  1:356047/150‑1 (MQ=255)
     aTGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCttttt                                                                                                                                       <  2:1977895/150‑1 (MQ=255)
      tGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCtttttt                                                                                                                                      >  1:2371400/1‑150 (MQ=255)
        tCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAaga                                                                                                                                                   <  2:2247028/135‑1 (MQ=255)
        tCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCtttttttt                                                                                                                                    <  1:332675/150‑1 (MQ=255)
        tCAGGAGTTCGAATCCCCTCACCCCGACCAGTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAaga                                                                                                                                                   >  1:2247028/1‑135 (MQ=255)
         cAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCAGCTCAAAGAGATGCCCTTTTTTTTg                                                                                                                                   <  1:584125/150‑1 (MQ=255)
            gAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGaaa                                                                                                                                <  1:1210069/150‑1 (MQ=255)
              gTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGTCATCTCAAAGAGGTGCCCTTTTTTTTGAAaca                                                                                                                              >  1:130499/1‑150 (MQ=255)
              gTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAaca                                                                                                                              >  1:1596742/1‑150 (MQ=255)
               ttCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGaaa                                                                                                                                <  2:744432/147‑1 (MQ=255)
               ttCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGaaa                                                                                                                                >  1:744432/1‑147 (MQ=255)
               ttCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAacac                                                                                                                             >  2:752754/1‑150 (MQ=255)
               ttCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAacac                                                                                                                             >  2:2164887/1‑150 (MQ=255)
                  gAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAGAAGGGCATCTCaaa                                                                                                                                                     <  2:1723170/123‑1 (MQ=255)
                  gAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAGAAGGGCATCTCaaa                                                                                                                                                     >  1:1723170/1‑123 (MQ=255)
                  gAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGACCAAGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCg                                                                                                                          >  1:635239/1‑150 (MQ=255)
                   aaTCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCgg                                                                                                                         <  1:2164887/150‑1 (MQ=255)
                    aTCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCggg                                                                                                                        >  1:179133/1‑150 (MQ=255)
                    aTCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCggg                                                                                                                        <  2:1441619/150‑1 (MQ=255)
                      ccccTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGgt                                                                                                                      <  2:179133/150‑1 (MQ=255)
                      ccccTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTGAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGgt                                                                                                                      <  2:2371400/150‑1 (MQ=255)
                       cccTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGgtg                                                                                                                     >  1:2238245/1‑150 (MQ=255)
                         cTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCggg                                                                                                                        <  2:1563371/145‑1 (MQ=255)
                         cTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCggg                                                                                                                        >  1:1563371/1‑145 (MQ=255)
                         cTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTa                                                                                                                   <  2:1860213/150‑1 (MQ=255)
                              cccGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCa                                                                                                              >  1:234956/1‑150 (MQ=255)
                               ccGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCTATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAg                                                                                                             >  1:1241250/1‑150 (MQ=255)
                                  aCCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAacac                                                                                                                             <  2:548998/131‑1 (MQ=255)
                                  aCCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAacac                                                                                                                             >  1:548998/1‑131 (MQ=255)
                                    cATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATg                                                                                                                                               >  1:441392/1‑111 (MQ=255)
                                    cATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATg                                                                                                                                               <  2:441392/111‑1 (MQ=255)
                                            gggTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCAACCATGTAGAAAAGGGCATCTGAAAGAGATGCCCTTATTTTTGAAACACTCGGGGAGTAGCGCAGTCCGGTAgcgcgc                                                                                                >  1:109047/1‑150 (MQ=255)
                                               tCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTg                                                                                             >  2:1607080/1‑150 (MQ=255)
                                                 gTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTTAACCACTCGGGGTTTAGCTCGATCCGGAAGCGCGCCTGCt                                                                                           >  1:1542340/1‑150 (MQ=18)
                                                     gCTCAGTTGGTAGAGCAGTTGGCTTTTAACCTATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTgg                                                                                       <  2:1241250/150‑1 (MQ=255)
                                                     gCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTgg                                                                                       <  1:752754/150‑1 (MQ=255)
                                                      cTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTCGAAAAAGGCATCTCAAAGGGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTggg                                                                                      >  2:1399723/1‑150 (MQ=255)
                                                       tCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGccc                                                                                                                                            <  1:326331/95‑1 (MQ=14)
                                                       tCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGccc                                                                                                                                            >  2:326331/1‑95 (MQ=14)
                                                          gTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCa                                                                                  <  1:1819464/150‑1 (MQ=255)
                                                          gTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCa                                                                                  <  2:2238245/150‑1 (MQ=255)
                                                                 gagCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAATGAGATGCCCTTTTTTTTGAAACCCTCGGGGTCTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTc                                                                           <  2:1842845/150‑1 (MQ=255)
                                                                 gagCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCgg                                                                                                                         >  1:1213134/1‑104 (MQ=255)
                                                                 gagCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCgg                                                                                                                         >  2:1382962/1‑104 (MQ=255)
                                                                 gagCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCgg                                                                                                                         <  1:1382962/104‑1 (MQ=255)
                                                                 gagCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCgg                                                                                                                         <  2:1213134/104‑1 (MQ=255)
                                                                       ttGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCggg                                                                                                                        <  2:1095472/99‑1 (MQ=255)
                                                                       ttGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCggg                                                                                                                        >  1:1095472/1‑99 (MQ=255)
                                                                         ggCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTc                                                                   <  1:1607080/150‑1 (MQ=255)
                                                                            ttttAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGATTTGGGAGCAGGGTGTCAGGAGGTCGaa                                                                >  2:1095147/1‑150 (MQ=255)
                                                                            ttttAACCAATTGGTCGGAGGT‑TCGAATCCTACACGACCAACCATGTAGAAAAGGGCATCTCAAAGACATGCCCTTTTTTTTGAAACACTTGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGaa                                                                >  2:1197760/1‑150 (MQ=255)
                                                                                  ccAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGCAAAGGGCACTTCAAAGAGATGCCCTTTTTCTTGAAACACTGGAAGTGTAGCGTAGACGGGTAGCGCGCCTACATTGGGATCAGTATGGGAGGGGGTCGAATccttt                                                          >  2:535021/1‑147 (MQ=18)
                                                                                  ccAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCt                                                          >  2:471439/1‑150 (MQ=255)
                                                                                   cAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCa                                                                          >  2:6076/1‑133 (MQ=255)
                                                                                   cAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCa                                                                          <  1:6076/133‑1 (MQ=255)
                                                                                   cAATTGGACGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAATACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTc                                                         >  2:1126147/1‑150 (MQ=255)
                                                                                    aaTTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCa                                                        >  1:1072134/1‑150 (MQ=255)
                                                                                         gTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCGGACGGAGATGCCCTTTTTTCTGAAAAACACGGCGTGTGGAGCACTCCGTTAGCGCCCCCGCTTTAAGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCg                                                   >  2:955540/1‑150 (MQ=21)
                                                                                              aGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAg                                                                                 <  1:1066816/115‑1 (MQ=255)
                                                                                              aGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAg                                                                                 >  2:1066816/1‑115 (MQ=255)
                                                                                               ggT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTAAGAGATGCCCTTTTTTTTGTAACACTCGGGGGGTAGCGCAGTCCGGTAGCGCGCCTGCTCTGGGAGCAGCATGGCAGGAGTTCGAATCCTCTCACCCCGACCAtt                                             >  2:1747447/1‑150 (MQ=25)
                                                                                                gT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGATTCCCCTCACCCCGACCAttt                                            >  2:1947232/1‑150 (MQ=255)
                                                                                                     gAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCGGTCCGGTAGCGCGCCCCCTTTGGGAGCAGGATGGCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTgg                                        <  2:1596742/150‑1 (MQ=255)
                                                                                                       aTCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGcc                                                                                               >  1:2137146/1‑93 (MQ=255)
                                                                                                       aTCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGcc                                                                                               <  2:2137146/93‑1 (MQ=255)
                                                                                                        tCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTc                                     <  1:35079/150‑1 (MQ=255)
                                                                                                        tCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTc                                     >  2:1544002/1‑150 (MQ=255)
                                                                                                        tCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAAATGCACTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTGCGAATCCCCTCACCCCCACCATTTTTGGGTc                                     <  1:1095147/150‑1 (MQ=255)
                                                                                                          cTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGt                                   >  2:1671585/1‑150 (MQ=255)
                                                                                                             cacGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTgg                                        >  1:810241/1‑142 (MQ=255)
                                                                                                             cacGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTgg                                        <  2:810241/142‑1 (MQ=255)
                                                                                                              acGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGAt                                                                              >  1:1675641/1‑103 (MQ=255)
                                                                                                              acGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGAt                                                                              <  2:1675641/103‑1 (MQ=255)
                                                                                                                 aCCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCa                            <  1:471439/150‑1 (MQ=255)
                                                                                                                   ccaccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATc                                                              >  1:1239502/1‑114 (MQ=255)
                                                                                                                   ccaccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATc                                                              <  2:1239502/114‑1 (MQ=255)
                                                                                                                     accaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTc                                                                                                           <  2:1458088/67‑1 (MQ=255)
                                                                                                                     accaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTc                                                                                                           >  1:1458088/1‑67 (MQ=255)
                                                                                                                     accaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTg                        >  2:764504/1‑150 (MQ=255)
                                                                                                                      ccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGgt                                                                                                                      <  2:640767/55‑1 (MQ=255)
                                                                                                                      ccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGgt                                                                                                                      >  1:640767/1‑55 (MQ=255)
                                                                                                                      ccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTgg                       <  1:1671585/150‑1 (MQ=255)
                                                                                                                      ccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTgg                       <  1:1126147/150‑1 (MQ=255)
                                                                                                                        aTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTa                     <  1:1197760/150‑1 (MQ=255)
                                                                                                                          aaaaaaaaGGGCATCTCTTTGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTaga                   >  1:626154/5‑150 (MQ=18)
                                                                                                                            aGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGc                 >  2:1692516/1‑150 (MQ=255)
                                                                                                                                aaGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGtt             >  2:372910/1‑150 (MQ=255)
                                                                                                                                aaGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGtt             >  2:2025424/1‑150 (MQ=255)
                                                                                                                                  gggCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTgg           >  1:2227776/1‑150 (MQ=255)
                                                                                                                                  gggCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCAGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTgg           >  1:1972053/1‑150 (MQ=255)
                                                                                                                                     gagctcGAAGAGATGCCCTTTTTTTTGCAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCtt        <  1:158746/147‑1 (MQ=255)
                                                                                                                                           aaaGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAAcc  >  2:2148673/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACC  >  NC_002947/733510‑733797

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