Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 587,023 A→C intergenic (‑115/+68) envY ← / ← ybcH porin thermoregulatory transcriptional activator/PRK09936 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913587,0230AC98.8% 340.2 / ‑4.0 85intergenic (‑115/+68)envY/ybcHporin thermoregulatory transcriptional activator/PRK09936 family protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (45/39);  total (46/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

TCACCACGCAAGGTTCACTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTAATCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAA  >  NC_000913/586875‑587166
                                                                                                                                                    |                                                                                                                                               
tCACCACGCAAGGTTCACTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTc                                                                                                                                               >  2:344259/1‑151 (MQ=255)
  accacGCAAGGTTCACTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTActct                                                                                                                                              >  1:7551/1‑150 (MQ=255)
   ccacGCAAGGTTCACTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCtt                                                                                                                                             <  1:629748/150‑1 (MQ=255)
    cacGCAAGGTTCACTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTg                                                                                                                                            >  2:90848/1‑150 (MQ=255)
       gCAAGGTTCACTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATgg                                                                                                                                        >  2:1171374/1‑151 (MQ=255)
          aGGTTCACTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAAt                                                                                                                                     <  2:542163/151‑1 (MQ=255)
               cACTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCtga                                                                                                                                <  2:363489/151‑1 (MQ=255)
                aCTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCtgat                                                                                                                               <  1:802465/151‑1 (MQ=255)
                aCTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCtgat                                                                                                                               >  1:99309/1‑151 (MQ=255)
                  tgctgcTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCtgatg                                                                                                                              >  2:1209392/1‑150 (MQ=255)
                  tgctgcTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCtgatg                                                                                                                              >  1:828651/1‑150 (MQ=255)
                    ctgctCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCtgatg                                                                                                                              >  2:695469/1‑148 (MQ=255)
                     tgctCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATgag                                                                                                                            >  2:1203802/1‑149 (MQ=255)
                     tgctCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATgag                                                                                                                            >  1:477435/1‑149 (MQ=255)
                     tgctCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATgag                                                                                                                            >  2:63133/1‑149 (MQ=255)
                      gctCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGagag                                                                                                                          <  1:275295/150‑1 (MQ=255)
                         cAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGa                                                                                                                       >  1:900504/1‑150 (MQ=255)
                          aaTTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGaa                                                                                                                      >  2:54106/1‑150 (MQ=255)
                                 tttCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTa                                                                                                               >  2:606917/1‑150 (MQ=255)
                                 tttCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTa                                                                                                               >  2:114143/1‑150 (MQ=255)
                                  ttCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTaa                                                                                                              <  2:627145/150‑1 (MQ=255)
                                  ttCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTa                                                                                                               >  2:1112665/1‑149 (MQ=255)
                                      cACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAg                                                                                                         <  2:407971/151‑1 (MQ=255)
                                        ctcctcAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGTTGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAgg                                                                                                        >  1:373721/1‑150 (MQ=255)
                                          cctcAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCg                                                                                                     >  1:380047/1‑151 (MQ=255)
                                          cctcAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCg                                                                                                     >  1:891040/1‑151 (MQ=255)
                                          cctcAGATATCAGAAACTCCACTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCg                                                                                                     >  1:634525/1‑151 (MQ=255)
                                              aGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTAcc                                                                                                 >  1:663113/1‑151 (MQ=255)
                                               gATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTAccc                                                                                                <  2:834092/151‑1 (MQ=255)
                                                atatCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTAccc                                                                                                <  2:1004194/150‑1 (MQ=255)
                                                atatCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCg                                                                                               <  1:654430/151‑1 (MQ=255)
                                                 tatCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGACGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCg                                                                                               >  2:921009/1‑150 (MQ=255)
                                                  atCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGAc                                                                                             <  2:1110681/151‑1 (MQ=255)
                                                     aGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACtt                                                                                           >  1:403302/1‑150 (MQ=255)
                                                     aGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTa                                                                                          <  2:132021/151‑1 (MQ=255)
                                                     aGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTa                                                                                          <  2:840292/151‑1 (MQ=255)
                                                      gAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCACGTTCGCCGATAACGATCTTCTTTCTTGAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTAc                                                                                         >  2:260193/1‑151 (MQ=255)
                                                       aaaCTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTAc                                                                                         <  1:995185/150‑1 (MQ=255)
                                                        aaCTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCt                                                                                       >  2:458849/1‑151 (MQ=255)
                                                          cTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTg                                                                                      <  2:1102157/150‑1 (MQ=255)
                                                           tCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTgg                                                                                     >  2:543557/1‑150 (MQ=255)
                                                               cTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTggagg                                                                                  >  1:639146/1‑149 (MQ=255)
                                                                tCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTggagg                                                                                  >  2:372177/1‑148 (MQ=255)
                                                                 cAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTggagg                                                                                  >  2:487295/1‑147 (MQ=255)
                                                                   aaGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGAttt                                                                            <  2:336243/151‑1 (MQ=255)
                                                                          tATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTc                                                                     >  1:890346/1‑151 (MQ=255)
                                                                              cTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCg                                                                    >  2:505218/1‑148 (MQ=255)
                                                                               ttCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCgaga                                                                 >  2:212984/1‑150 (MQ=255)
                                                                                tCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAAt                                                               >  1:449568/1‑151 (MQ=255)
                                                                                 ccTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCTACTTACCTGGAGGAGATTTAATACTCGAGAATg                                                              <  2:504196/151‑1 (MQ=255)
                                                                                 ccTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAAt                                                               <  2:1080948/150‑1 (MQ=255)
                                                                                 ccTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATg                                                              <  1:273159/151‑1 (MQ=255)
                                                                                   tGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGcc                                                            >  1:96431/1‑151 (MQ=255)
                                                                                    gCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCg                                                           <  2:18303/151‑1 (MQ=255)
                                                                                     cATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGt                                                          <  1:66706/151‑1 (MQ=255)
                                                                                        gAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTgcg                                                       <  2:20891/151‑1 (MQ=255)
                                                                                               cGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACggg                                                >  2:290984/1‑151 (MQ=255)
                                                                                               cGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACggg                                                <  1:1074650/151‑1 (MQ=255)
                                                                                                ggCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGc                                               >  1:1134123/1‑151 (MQ=255)
                                                                                                ggCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGc                                               <  1:1209470/151‑1 (MQ=255)
                                                                                                       tCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAg                                         <  2:783196/150‑1 (MQ=255)
                                                                                                       tCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAg                                         <  2:573741/150‑1 (MQ=255)
                                                                                                              tAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCag                                    >  1:922020/1‑148 (MQ=255)
                                                                                                               aaCGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCagaga                                 <  1:786117/150‑1 (MQ=255)
                                                                                                                  gATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGaa                              >  1:904716/1‑150 (MQ=255)
                                                                                                                    tcttctTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGGAAt                            >  1:422324/1‑150 (MQ=255)
                                                                                                                     cttcttTCTTTAGCACGCTTTTTAGCAATTAATCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAAt                            >  2:1193165/1‑149 (MQ=255)
                                                                                                                       tctttcttTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAAt                            <  2:1132347/147‑1 (MQ=255)
                                                                                                                       tctttcttTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATag                          <  1:968201/149‑1 (MQ=255)
                                                                                                                        ctttctttAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATag                          >  1:696659/1‑148 (MQ=255)
                                                                                                                          ttctttAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATagaga                       <  2:278869/149‑1 (MQ=255)
                                                                                                                          ttctttAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAgt                     >  1:55858/1‑151 (MQ=255)
                                                                                                                          ttctttAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGGAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATagagag                      <  1:125590/150‑1 (MQ=255)
                                                                                                                               tAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCt                 <  2:363478/150‑1 (MQ=255)
                                                                                                                                aGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCttt               <  1:597971/151‑1 (MQ=255)
                                                                                                                                 gCATGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTa              <  2:215469/151‑1 (MQ=255)
                                                                                                                                  cACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTAc             >  2:777204/1‑151 (MQ=255)
                                                                                                                                     gCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGc          >  2:826512/1‑151 (MQ=255)
                                                                                                                                        ttttAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCa         <  2:477412/149‑1 (MQ=255)
                                                                                                                                        ttttAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCa         >  2:1137663/1‑149 (MQ=255)
                                                                                                                                         tttAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAag       <  1:838251/150‑1 (MQ=255)
                                                                                                                                           tAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAgg    >  2:834880/1‑151 (MQ=255)
                                                                                                                                           tAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAgg    >  2:873882/1‑151 (MQ=255)
                                                                                                                                            aGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAgg    <  1:344271/150‑1 (MQ=255)
                                                                                                                                            aGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAgg    <  1:114150/150‑1 (MQ=255)
                                                                                                                                             gCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGaa  <  1:469871/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TCACCACGCAAGGTTCACTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTAATCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAA  >  NC_000913/586875‑587166

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