Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 733,649 A→T 7.9% intergenic (+20/‑27) 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,6490AT 7.9% 249.5 / 5.6 90intergenic (+20/‑27)PP_t09/PP_t10tRNA‑Lys/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base A (45/37);  new base T (4/4);  total (49/41)
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

TGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACC  >  NC_002947/733500‑733797
                                                                                                                                                     |                                                                                                                                                     
tGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCaa                                                                                                                                                       <  2:4230152/150‑1 (MQ=255)
 gCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCaaa                                                                                                                                                      <  1:724020/150‑1 (MQ=255)
 gCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCaaa                                                                                                                                                      <  2:2295900/150‑1 (MQ=255)
 gCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCaaa                                                                                                                                                      <  2:2965227/150‑1 (MQ=255)
 gCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCaaa                                                                                                                                                      <  1:2713241/150‑1 (MQ=255)
 gCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCaaa                                                                                                                                                      <  1:724372/150‑1 (MQ=255)
 gCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCaaa                                                                                                                                                      <  1:1252379/150‑1 (MQ=255)
   tttGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAaga                                                                                                                                                    <  2:1007952/150‑1 (MQ=255)
      gggAGGAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGAt                                                                                                                                                 <  2:3563146/150‑1 (MQ=255)
      gggAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGAt                                                                                                                                                 <  1:2012762/150‑1 (MQ=255)
          gCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGTTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTTACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGccc                                                                                                                                             <  2:624151/150‑1 (MQ=255)
          gCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGccc                                                                                                                                             <  2:4134291/150‑1 (MQ=255)
          gCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTAGTAGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGccc                                                                                                                                             <  1:3487821/150‑1 (MQ=255)
           cAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTGAGCAATTGGTCGAACGTGCGAATACAAGACGACGCTCAATTTAGATTAGGGCGTCTCATGGAGTTGCCCt                                                                                                                                            >  1:1982519/1‑150 (MQ=14)
            aGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCtt                                                                                                                                           <  1:1306803/150‑1 (MQ=255)
            aGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCtt                                                                                                                                           <  1:1306787/150‑1 (MQ=255)
             ggATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCtt‑t                                                                                                                                         <  2:1677296/150‑1 (MQ=255)
              gATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCACAAACAAATGCCCtt‑tt                                                                                                                                        >  1:4863398/1‑150 (MQ=255)
              gATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGGCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGAGGCCCtt‑tt                                                                                                                                        >  1:2039512/1‑150 (MQ=255)
                 gTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCtt‑ttttt                                                                                                                                     <  1:1762840/150‑1 (MQ=255)
                    aGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAAAAAAGGGCATCTCTTTGAGATGCCCtt‑ttctacat                                                                                                                                  <  1:4845094/150‑7 (MQ=14)
                    aGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCGCAGTTGGTGGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTGCGAAACCTACACGACCCACCATGTAGAAATGGGCATCTCGAAGAGATGCCCTT‑TTTTTTGa                                                                                                                                  >  1:356149/1‑150 (MQ=17)
                     ggAGTTCGTATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGaa                                                                                                                                 <  2:3876035/150‑1 (MQ=255)
                     ggAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGaa                                                                                                                                 <  2:3076886/150‑1 (MQ=255)
                     ggAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGaa                                                                                                                                 <  1:5122013/150‑1 (MQ=255)
                     ggAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGaa                                                                                                                                 <  2:3523778/150‑1 (MQ=255)
                         ttCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAacac                                                                                                                             >  2:549985/1‑150 (MQ=255)
                         ttCGAATCCACTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTAAACGACCCACCATGTAGAAAAGGGAATCTCAAACAGATCCCCTTATTTTTTGAAaca                                                                                                                              >  1:3481278/1‑150 (MQ=18)
                          tCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGTTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAATGGGAATCTCAAATTGATGCACTT‑TTTTTTGAAACACt                                                                                                                            >  2:4190400/1‑150 (MQ=18)
                            gAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGGCCTT‑TTTTTTGAAACACTCg                                                                                                                          >  2:4945220/1‑150 (MQ=255)
                            gAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGTAACACTCg                                                                                                                          >  1:4672957/1‑150 (MQ=255)
                            gAATCCCCTCACCCCGACCATTTTCCGGGGCGTAAGCTCAGTTGGTAGAGCAGTTGGATTTAAACCAATGGGTCGTAGGTTCGAATCCTACACGTGCCCGCAGGTAGAAATGGGCAGCTCGTGGAGGTGCCCTT‑TGTTTTGAATCACTCg                                                                                                                          >  2:4740704/1‑150 (MQ=17)
                             aaTCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGTAAAGGGCATCTCAAAGAGATGCCCTT‑TTGTGTGAAACACTCgg                                                                                                                         >  1:4520351/1‑150 (MQ=255)
                              aTCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCggg                                                                                                                        <  2:852398/150‑1 (MQ=255)
                               tCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCAAGATGTAGAAAAGGGAATATCAAAGAGATGCCCTT‑TTTTTTGAAAAACTCgggg                                                                                                                       >  2:2004814/1‑150 (MQ=18)
                                ccccTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGgt                                                                                                                      <  1:4224224/150‑1 (MQ=255)
                                ccccTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGgt                                                                                                                      <  2:1526079/150‑1 (MQ=255)
                                 cccTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAACAGAGGCCCTT‑TATGATGAAACACTCGGGgtg                                                                                                                     >  1:2371105/1‑150 (MQ=18)
                                   cTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTa                                                                                                                   <  1:2117958/150‑1 (MQ=255)
                                      aCCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCAGGTAGAAAAGGGCATCACAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAgag                                                                                                                >  1:647825/1‑148 (MQ=255)
                                       ccccGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAgcgc                                                                                                               >  1:2185905/1‑150 (MQ=255)
                                        cccGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGTCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACACGGGGTGTAGCGCa                                                                                                              >  2:1303441/1‑150 (MQ=255)
                                        cccGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGATATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCa                                                                                                              <  1:622478/150‑1 (MQ=255)
                                        cccGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTGAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCa                                                                                                              <  2:1980979/150‑1 (MQ=255)
                                          cGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGt                                                                                                            >  1:640966/1‑150 (MQ=255)
                                            aCCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGTAGTGGGCAGCAAAAAGAGATCCCCAT‑TGTGTTGTAACGCTCGGGGTGTCgcgcggggc                                                                                                          >  1:1736213/1‑145 (MQ=17)
                                                          cGTCAGCCCCGTTGGGAGCGCAGTGGGCTGTGACCCGTTGGGCGGTAGGTTCGACTGCTCCCCGACCCACCCTGTAGAAATGGGCACCTGCAAGCGCTGCCGTT‑CTTTTTGACACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGc                                                                                            <  2:2517931/150‑1 (MQ=255)
                                                            ttAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGAGCGCCTGCtt                                                                                          >  2:4408255/1‑150 (MQ=255)
                                                              aGCTCAGTTGGTCGAGCAGTTGTCTTTTTACCAATTTGTCGTAGTTTCGTATCCTACACGCCCCTCCTTTTAGAATAGGGTCTCTCCTAGTGATGCACTT‑TTTTTTGAAACTCTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTg                                                                                        <  2:4665452/150‑1 (MQ=25)
                                                               tcTCAGCTGTTAGCGCATTTGGCTTTTAACCCATTGGTGGTAGGTCCGAATCCTCCTCTACCGACCATGTATAAAAGGGCATCTCCTAGCTATGCCGTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTgg                                                                                       <  1:333139/149‑1 (MQ=25)
                                                               gCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTgg                                                                                       <  2:14867/150‑1 (MQ=255)
                                                               gCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTgg                                                                                       <  1:4205794/150‑1 (MQ=255)
                                                               gCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTgg                                                                                       <  1:1495093/150‑1 (MQ=255)
                                                               gCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTgg                                                                                       <  1:4446323/150‑1 (MQ=255)
                                                                cTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCTTGCTTTggg                                                                                      <  1:1462010/150‑1 (MQ=255)
                                                                 tCAGTTGGTAGATCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGTATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGa                                                                                     <  2:3818603/150‑1 (MQ=255)
                                                                 tCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCAGGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGGAGCGCGCCTGCTTTGGGa                                                                                     >  2:2788692/1‑150 (MQ=255)
                                                                    cttGGTAGAGCAGTTGGATTTTAACCAATTGGTCGTATGTTCGAATCCTACACGACCCACCATGTATAAAAGGGCATCTCAAAGAGATGCACTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCa                                                                                  <  1:1282204/149‑1 (MQ=255)
                                                                         tAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGACACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATg                                                                             <  2:1324556/150‑1 (MQ=255)
                                                                           gaggagTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCCTGTAGAAAAGGGGATCTCAAAGATATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTc                                                                           <  2:4136256/150‑1 (MQ=255)
                                                                               aGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGGAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGa                                                                       <  2:4053805/150‑1 (MQ=255)
                                                                                   ggCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTc                                                                   >  2:2193043/1‑150 (MQ=255)
                                                                                      ttttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGaa                                                                >  1:505064/1‑150 (MQ=255)
                                                                                      ttttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGaa                                                                >  1:636677/1‑150 (MQ=255)
                                                                                         tAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCTGGGTGTAGCGCAGTCCGGTAGAGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATcc                                                             >  2:1491275/1‑150 (MQ=255)
                                                                                           aCCAATTGGTCGTAGGTTCGAATGCTACTCGACCCTCCATGTAGAAAAGGGGATCTCAAAGTGTTGCAGTT‑ATTTTTGAAACTCTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATcccc                                                           <  2:3652696/150‑1 (MQ=21)
                                                                                             cagttGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGTAAAGGGCATCTCAATGAGATGCCCTT‑TTGTTTGAAACACTCGGGGTGTAGCGCGGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTc                                                         >  1:629077/4‑150 (MQ=37)
                                                                                             cAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTc                                                         >  1:4038594/1‑150 (MQ=255)
                                                                                              aaTTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCACAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCACTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCa                                                        >  2:1187185/1‑150 (MQ=255)
                                                                                                ttGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAATGGGCATCACAAAGAGATGCCCTT‑TTTTTTGAAACACTAGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCATAGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCAcc                                                      >  2:4391816/1‑150 (MQ=37)
                                                                                                        aGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCAt                                              >  1:2022182/1‑150 (MQ=255)
                                                                                                           ttCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑T‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTCGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCAttttt                                          >  2:1337376/1‑150 (MQ=255)
                                                                                                           ttCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑T‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCAttttt                                          >  2:1337364/1‑150 (MQ=255)
                                                                                                             cgtccccTCCACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTg                                         <  2:4205752/145‑1 (MQ=255)
                                                                                                              gAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTgg                                        >  2:1323174/1‑150 (MQ=255)
                                                                                                               aaTCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTggg                                       >  1:494391/1‑150 (MQ=255)
                                                                                                               aaTCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTGTTggg                                       >  2:640247/1‑150 (MQ=255)
                                                                                                                  ccTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCg                                    >  2:4365578/1‑150 (MQ=255)
                                                                                                                   cTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGt                                   >  2:620227/1‑150 (MQ=255)
                                                                                                                     acacGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTgg                                                                                       <  1:2463771/96‑1 (MQ=255)
                                                                                                                     acacGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTgg                                                                                       >  2:2463771/1‑96 (MQ=255)
                                                                                                                        cGACCCACCCTGTAGAACAGGGCATCTCAAAGAGATGCACTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCggggtga                              <  1:2640332/150‑7 (MQ=21)
                                                                                                                            ccaccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCACAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGt                          >  1:3437617/1‑150 (MQ=255)
                                                                                                                             caccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGtt                         >  2:1278065/1‑150 (MQ=255)
                                                                                                                               ccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTgg                       >  2:1131369/1‑150 (MQ=255)
                                                                                                                                ttcaaaaaaaaGGGCATCTCTTTGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGt                      >  1:1461195/8‑150 (MQ=18)
                                                                                                                                   gTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTaga                   >  2:297997/1‑150 (MQ=255)
                                                                                                                                      gAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGcc                >  1:1754645/1‑149 (MQ=255)
                                                                                                                                         aaGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGtt             <  1:549985/150‑1 (MQ=255)
                                                                                                                                            ggCATCTCTTTGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGc          >  2:721295/1‑150 (MQ=255)
                                                                                                                                              cATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCTGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCtt        <  2:505064/150‑1 (MQ=255)
                                                                                                                                              cATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTGCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCtt        >  2:4379821/1‑150 (MQ=255)
                                                                                                                                                tctcAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCtttt      >  1:887585/1‑150 (MQ=255)
                                                                                                                                                 ctcAACGAGATGCCCTT‑ATTTTTGCAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAAGCCCCACACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTa     <  2:4520351/150‑1 (MQ=17)
                                                                                                                                                  tcAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTat    >  1:126078/1‑149 (MQ=255)
                                                                                                                                                  tcAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTaa    >  2:3141219/1‑150 (MQ=255)
                                                                                                                                                  tcAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTaa    >  2:1053937/1‑150 (MQ=255)
                                                                                                                                                    aaaGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAAcc  >  1:2342279/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTT‑TTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACC  >  NC_002947/733500‑733797

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