Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,103,527 T→A 5.5% intergenic (+31/+75) yiiM → / ← cpxA 6‑N‑hydroxylaminopurine resistance protein/sensory histidine kinase in two‑component regulatory system with CpxR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,103,5270TA 5.5% 290.2 / 5.1 91intergenic (+31/+75)yiiM/cpxA6‑N‑hydroxylaminopurine resistance protein/sensory histidine kinase in two‑component regulatory system with CpxR
Reads supporting (aligned to +/- strand):  ref base T (48/38);  new base A (3/2);  total (51/40)
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.91e-01

TCAATATCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCC  >  NC_000913/4103379‑4103675
                                                                                                                                                    |                                                                                                                                                    
tCAATATCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTAt                                                                                                                                                      >  2:585726/1‑149 (MQ=255)
  aaTATCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCa                                                                                                                                                    <  1:744295/149‑1 (MQ=255)
    tatCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGc                                                                                                                                                  <  2:942343/149‑1 (MQ=255)
     atCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCg                                                                                                                                                 >  1:1066800/1‑149 (MQ=255)
      tCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGctt                                                                                                                                                >  1:680021/1‑147 (MQ=255)
       cACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTc                                                                                                                                               <  1:21174/149‑1 (MQ=255)
        aCCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCATAAGCGCCGACTGAGCGGCAAGATCGAAGATTATTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCg                                                                                                                                              <  1:262280/149‑1 (MQ=255)
        aCCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCg                                                                                                                                              >  2:536795/1‑149 (MQ=255)
         ccGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGc                                                                                                                                             <  2:413179/149‑1 (MQ=255)
          cGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTGATAGCGTctta                                                                                                                                            <  2:749004/149‑4 (MQ=255)
            tttACTCTACGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCAGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATc                                                                                                                                          >  2:1134272/1‑149 (MQ=255)
             ttACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCa                                                                                                                                         >  1:57101/1‑149 (MQ=255)
              tACTCTCCGCTTCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAg                                                                                                                                        <  2:808741/149‑1 (MQ=255)
              tACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTATGCCTGATGCGACGCTGATAGCGTCTTATCAg                                                                                                                                        <  2:441644/149‑1 (MQ=255)
                ctctCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGATGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGc                                                                                                                                      >  2:788819/1‑149 (MQ=255)
                 tctcCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCa                                                                                                                                     >  1:1477218/1‑149 (MQ=255)
                   tcCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATc                                                                                                                                   <  1:387975/149‑1 (MQ=255)
                       cTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGc                                                                                                                               >  2:615008/1‑149 (MQ=255)
                         gCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTc                                                                                                                             <  1:585726/149‑1 (MQ=255)
                             ggTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAAt                                                                                                                         >  1:208640/1‑149 (MQ=255)
                                 aTCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCg                                                                                                                     <  2:1494992/149‑1 (MQ=255)
                                  tCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCg                                                                                                                     <  1:60429/148‑1 (MQ=255)
                                  tCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCg                                                                                                                     >  2:60429/1‑148 (MQ=255)
                                    gAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCCGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCAGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGAt                                                                                                                  >  2:931961/1‑149 (MQ=255)
                                     aaaaGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATc                                                                                                                                   >  2:403152/1‑131 (MQ=255)
                                     aaaaGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATc                                                                                                                                   <  1:403152/131‑1 (MQ=255)
                                       aaGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATgcg                                                                                                               <  2:862932/149‑1 (MQ=255)
                                        aGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATgcgg                                                                                                              <  1:1204314/149‑1 (MQ=255)
                                        aGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATgcgg                                                                                                              <  1:1229865/149‑1 (MQ=255)
                                          cTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATgcggcg                                                                                                            <  1:2845/149‑1 (MQ=255)
                                           tgtaagcgGACGATGAAAAAGCGCAGACTGAGCGGCAAGTTCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCTAATGCCGGATGCGGCGt                                                                                                           <  2:1211046/144‑1 (MQ=39)
                                             gACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGCTAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTaa                                                                                                         >  2:1220264/1‑149 (MQ=255)
                                              aCGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTaaa                                                                                                        <  1:1566411/149‑1 (MQ=255)
                                                 cgGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGc                                                                                                     <  1:1354732/149‑1 (MQ=255)
                                                   gACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCt                                                                                                   >  2:570803/1‑149 (MQ=255)
                                                       aTGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATc                                                                                               <  1:874868/149‑1 (MQ=255)
                                                          cAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATcctg                                                                                            <  1:814902/149‑1 (MQ=255)
                                                           aaaaGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATcctgc                                                                                           >  2:1232099/1‑149 (MQ=255)
                                                            aaaGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATcctgcc                                                                                          >  1:48526/1‑149 (MQ=255)
                                                                cgcCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCa                                                                                      <  2:89267/149‑1 (MQ=255)
                                                                 gcCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCaa                                                                                     <  1:866630/149‑1 (MQ=255)
                                                                  ccGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCaaa                                                                                    >  1:317064/1‑149 (MQ=255)
                                                                   cGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAAt                                                                                   >  1:544079/1‑149 (MQ=255)
                                                                   cGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAAt                                                                                   >  2:435404/1‑149 (MQ=255)
                                                                    gACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATg                                                                                  >  1:1575306/1‑149 (MQ=255)
                                                                     aCTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGc                                                                                 >  2:1049020/1‑149 (MQ=255)
                                                                       tGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCGTATCCTGCCTGCAAATGCGa                                                                               >  1:203201/1‑149 (MQ=255)
                                                                        gAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGaa                                                                              <  2:833768/149‑1 (MQ=255)
                                                                        gAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGaa                                                                              <  1:47681/149‑1 (MQ=255)
                                                                         aGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAg                                                                             >  2:111295/1‑149 (MQ=255)
                                                                         aGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAg                                                                             >  2:150860/1‑149 (MQ=255)
                                                                           cGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGtt                                                                           <  2:298602/149‑1 (MQ=255)
                                                                             gCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTa                                                                         >  1:1163516/1‑149 (MQ=255)
                                                                              cAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTaa                                                                        >  1:1384846/1‑149 (MQ=255)
                                                                              cAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTaa                                                                        <  1:592365/149‑1 (MQ=255)
                                                                              cAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTaa                                                                        >  1:1160202/1‑149 (MQ=255)
                                                                                 gATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTc                                                                     >  2:308038/1‑149 (MQ=255)
                                                                                  aTCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTcc                                                                    >  2:396811/1‑149 (MQ=255)
                                                                                   tCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCg                                                                   >  1:635016/1‑149 (MQ=255)
                                                                                      aaGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAACGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCtt                                                                >  2:1438181/1‑149 (MQ=255)
                                                                                             ttCCCGCCTGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGc                                                         <  1:1134272/149‑1 (MQ=255)
                                                                                                  cccgGTTGTTGGGAAAATAACACCCGGATTGTATGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCaa                                                    <  2:255276/148‑1 (MQ=255)
                                                                                                           ggggAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCAc                                           <  2:57101/149‑1 (MQ=255)
                                                                                                             ggAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCa                                         <  1:656573/149‑1 (MQ=255)
                                                                                                              gAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAg                                        >  1:476519/1‑149 (MQ=255)
                                                                                                                   tAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTa                                   >  2:888122/1‑149 (MQ=255)
                                                                                                                     aCCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTaaa                                 <  2:1477218/149‑1 (MQ=255)
                                                                                                                      cccccGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAAc                                >  2:1115936/1‑149 (MQ=255)
                                                                                                                        cccGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCGTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTactccg                              >  1:250777/1‑144 (MQ=255)
                                                                                                                        cccGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGGTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCg                              >  1:250778/1‑149 (MQ=255)
                                                                                                                               tgtAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCg                       >  2:754594/1‑149 (MQ=255)
                                                                                                                                 tAGGCCTGAT‑AGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCGGCCGTAAACCGCCCAGCGGGc                    >  1:1361571/1‑149 (MQ=255)
                                                                                                                                 gatgCCTGATGCGACGCTGATAGCGTCTTATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCggg                     >  2:334011/4‑149 (MQ=38)
                                                                                                                                  atgCCTGATGCGACGCTGATAGCGTCTTATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGc                    >  1:1281382/3‑149 (MQ=38)
                                                                                                                                   tgCCTGATGCGACGCTGATAGCGTCTTATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCt                   >  2:1398240/2‑149 (MQ=38)
                                                                                                                                   ggCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCt                   >  1:877562/1‑149 (MQ=255)
                                                                                                                                    gCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCACCCAAATCACCAGCCGTAATCCGCCCAGCGGGCTg                  >  1:509692/1‑149 (MQ=255)
                                                                                                                                    gCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTg                  >  1:834797/1‑149 (MQ=255)
                                                                                                                                     ccTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGt                 >  2:1388565/1‑149 (MQ=255)
                                                                                                                                      cTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtc                <  1:615008/149‑1 (MQ=255)
                                                                                                                                      cTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtc                >  2:188652/1‑149 (MQ=255)
                                                                                                                                       tGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCTGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtct               <  1:788819/149‑1 (MQ=255)
                                                                                                                                       tGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtct               <  1:527189/149‑1 (MQ=255)
                                                                                                                                       tGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtct               <  1:920352/149‑1 (MQ=255)
                                                                                                                                       tGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtct               >  2:1182401/1‑149 (MQ=255)
                                                                                                                                        gATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtctt              >  1:1454720/1‑149 (MQ=255)
                                                                                                                                         aTAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtcttc             >  2:508233/1‑149 (MQ=255)
                                                                                                                                          tAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTAAACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtcttct            <  2:48526/149‑1 (MQ=255)
                                                                                                                                             gACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGcc         <  1:210714/149‑1 (MQ=255)
                                                                                                                                             gACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGcc         >  1:1318388/1‑149 (MQ=255)
                                                                                                                                                  tATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAACTCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCAc    >  1:795006/1‑149 (MQ=255)
                                                                                                                                                   aTCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCAcc   <  2:1575306/149‑1 (MQ=255)
                                                                                                                                                    tCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAACTCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCAccc  >  1:827350/1‑149 (MQ=255)
                                                                                                                                                    tCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCAccc  >  2:1213496/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCAATATCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCC  >  NC_000913/4103379‑4103675

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