Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,189,608 G→A 5.3% intergenic (+35/‑178) rpoC → / → yjaZ RNA polymerase, beta prime subunit/stationary phase growth adaptation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,189,6080GA 5.3% 256.9 / 4.0 94intergenic (+35/‑178)rpoC/yjaZRNA polymerase, beta prime subunit/stationary phase growth adaptation protein
Reads supporting (aligned to +/- strand):  ref base G (50/38);  new base A (2/3);  total (52/42)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

TCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGATTCGATACTA  >  NC_000913/4189462‑4189756
                                                                                                                                                  |                                                                                                                                                    
tCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTcggc                                                                                                                                                    <  2:507761/149‑1 (MQ=255)
tCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTcggc                                                                                                                                                    <  2:109128/149‑1 (MQ=255)
 cGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGCCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATACCGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTcggcg                                                                                                                                                   >  2:1043026/1‑149 (MQ=255)
  gCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTcggcgg                                                                                                                                                  <  1:470039/149‑1 (MQ=255)
    cGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGt                                                                                                                                                <  1:1132643/149‑1 (MQ=255)
     gTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGtt                                                                                                                                               <  2:449878/149‑1 (MQ=255)
      ggctgcGGGGGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTGCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAACACCCGCCGAAGCGGGttt                                                                                                                                              <  2:700001/148‑1 (MQ=255)
       gctgcGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGCACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGtttt                                                                                                                                             >  2:293010/1‑149 (MQ=255)
        ctgcGGGTGAAGCTCCGGCTGCACCGCAGGTGACTTCAGAATACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGttttt                                                                                                                                            <  1:1009137/149‑1 (MQ=255)
        ctgcGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGttttt                                                                                                                                            <  2:826691/149‑1 (MQ=255)
         tgcGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGtttttt                                                                                                                                           >  1:433568/1‑149 (MQ=255)
         ggcGGGGGAAGCGCCGGCTTCACCGCAGTTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTACTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCCGAAGCGGGttttta                                                                                                                                           <  1:760496/148‑2 (MQ=255)
          gcGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCCGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCCGAAGCGGGtttttac                                                                                                                                          <  1:183019/149‑3 (MQ=255)
            gggTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTAt                                                                                                                                        >  1:1175873/1‑149 (MQ=255)
             ggTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGGCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATg                                                                                                                                       <  1:1137174/149‑1 (MQ=255)
              gTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATgg                                                                                                                                      >  2:760461/1‑149 (MQ=255)
              gTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATgg                                                                                                                                      <  2:616550/149‑1 (MQ=255)
              gTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATgg                                                                                                                                      <  1:972422/149‑1 (MQ=255)
                gAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATgggg                                                                                                                                    >  1:347255/1‑149 (MQ=255)
                gAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATgggg                                                                                                                                    <  1:166131/149‑1 (MQ=255)
                gAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATgggg                                                                                                                                    >  2:18640/1‑149 (MQ=255)
                gAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGGCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATgggg                                                                                                                                    <  2:373473/149‑1 (MQ=255)
                   gCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGa                                                                                                                                 >  2:244128/1‑149 (MQ=255)
                        ggCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTa                                                                                                                            <  1:382121/149‑1 (MQ=255)
                           tgcaaagCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATACGCAAATAAAGTAAAAACCCGCTTCGGCGGCTTTATTTATGGGGGGAGTTTAggg                                                                                                                         >  2:719099/7‑149 (MQ=255)
                                cGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCAGCATATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAtga                                                                                                                    >  2:307464/1‑149 (MQ=255)
                                 gCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAagag                                                                                                                   >  1:997286/1‑149 (MQ=255)
                                 gCAGGTGACAGGATTAGCCGCATCTGCCAGCCAGGCAGTACTGCTGCCGGCAGGTCTGGGCGGGGGAGATAACGAGTTATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCTGCGGGTTTTTTTATGGGGGGAGTTTACGGATagag                                                                                                                   <  2:1096289/149‑1 (MQ=255)
                                  cAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGc                                                                                                                  >  2:989070/1‑149 (MQ=255)
                                  cAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGc                                                                                                                  <  2:501785/149‑1 (MQ=255)
                                          tGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTGTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCa                                                                                                          >  2:1079001/1‑149 (MQ=255)
                                            cAGAAGACGCATCTGCCATCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATAGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTGTAGGGAAAGAGCATTTGTCAGa                                                                                                        <  1:527782/149‑1 (MQ=255)
                                                  aCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAATTAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATAttt                                                                                                  <  1:67126/149‑1 (MQ=255)
                                                    gCATCTGCCCGCCTGGCAGAACTTCTGAACGCAGGTCGGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTaa                                                                                                <  2:273907/149‑1 (MQ=255)
                                                    gCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTaa                                                                                                <  1:441187/149‑1 (MQ=255)
                                                    gCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTaa                                                                                                >  1:201029/1‑149 (MQ=255)
                                                    gCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGCTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCGTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTaa                                                                                                >  1:643733/1‑149 (MQ=255)
                                                          gCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAAtt                                                                                          >  2:1046322/1‑149 (MQ=255)
                                                            cAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCt                                                                                       >  2:246769/1‑149 (MQ=255)
                                                                cTGGCAGAACTGCTGAACGAAGGTCTGGGCGGTTCTGATAACGAGTAATAGTTAATCCGCAAATAACTTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCACTTGTCAGAATATTTAAGGAATTTCTGaa                                                                                    <  1:1159214/149‑1 (MQ=255)
                                                                 tGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTTATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAAt                                                                                   >  2:653356/1‑149 (MQ=255)
                                                                  ggCATAACTGGTGAACGCAGGTCTGTGCGGTTCTGATAACTAGTAATCGTTAATCCGCAAATAACGTAAAATCCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTATGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATa                                                                                  <  2:1043328/149‑1 (MQ=255)
                                                                   gCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATAc                                                                                 <  2:1098134/149‑1 (MQ=255)
                                                                            ctgAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCa                                                                        >  1:343097/1‑149 (MQ=255)
                                                                            ctgAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCa                                                                        <  1:363865/149‑1 (MQ=255)
                                                                               aaCGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATg                                                                     >  2:592662/1‑149 (MQ=255)
                                                                               aaCGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATg                                                                     >  1:1012015/1‑149 (MQ=255)
                                                                                 cgaagttCTTGGAGGTTCTGATAACGAGTAATCGTTAATCCTAAAATAACGTAAAAACCCGATTCGGCGGGTTTTTTTATGGGGGGAGTTTATGGAAATAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTa                                                                   <  1:18640/143‑1 (MQ=255)
                                                                                 cGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTa                                                                   <  2:198004/149‑1 (MQ=255)
                                                                                 cGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTa                                                                   <  2:525945/149‑1 (MQ=255)
                                                                                    aGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGTGTGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAAGTTCTGAATACTCATCATCAAAGGagag                                                                >  2:1089740/1‑149 (MQ=255)
                                                                                      gTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGCGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGAt                                                              >  1:133768/1‑149 (MQ=255)
                                                                                       tCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATg                                                             <  2:209578/149‑1 (MQ=255)
                                                                                        cTGTGCGGTTCTTATATCGAGTTATAGTTAATCCGCAAATTACGTAAAAACCCGCTTCGGCGTGTTTTTTTATGGGGGGAGTTTAGTGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGa                                                            <  1:773399/149‑1 (MQ=255)
                                                                                         tGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGAc                                                           <  1:1083738/149‑1 (MQ=255)
                                                                                         tGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTATTTATGGGGGCAGTTTAGGGAAACAGCATTTGTCAGAATATTAAAGGAATTTCAGAATACTCATAATCAATGTAGAGATGAc                                                           >  2:1109850/1‑149 (MQ=255)
                                                                                           ggCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTa                                                         >  1:1059549/1‑149 (MQ=255)
                                                                                           ggCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTa                                                         >  2:107506/1‑149 (MQ=255)
                                                                                              ggTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAata                                                      <  1:713784/149‑1 (MQ=255)
                                                                                                 tCTGATAACGAGTAATCGGTAATCCGCAAATAAAGTAAAAACCCGCTTAGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATcc                                                   <  1:645972/149‑1 (MQ=255)
                                                                                                   tGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATAAATGTAGAGATTACTAATATCCTg                                                 >  1:224937/1‑149 (MQ=255)
                                                                                                   gattAAAGATTACTCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGCGCATTTGTCAGAATATTAAAGGACATTCTGAATACTCATAATCAATCTAGAGATGACTAATATCCTg                                                 >  2:643711/4‑149 (MQ=255)
                                                                                                    gATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGa                                                >  1:261062/1‑149 (MQ=255)
                                                                                                     aTAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTATAGATGACTAATATCCTGaa                                               >  2:230952/1‑149 (MQ=255)
                                                                                                     aTAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGaa                                               <  1:909295/149‑1 (MQ=255)
                                                                                                     aTAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGaa                                               >  2:665374/1‑149 (MQ=255)
                                                                                                     aTAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGaa                                               >  2:686241/1‑149 (MQ=255)
                                                                                                     aTAACGAGTAATCGTTAAACCGCAAATAACGTAAAAACCAGATTCGGCTGGTTTTTTTATGGGGTGAGTTTAGTGAAAGAGCATTTGTCATAATATTTTAGGAATTTATGAATACTCAAAATCAATGTAGAGATGAATAATATcaagaa                                               >  2:642616/1‑144 (MQ=255)
                                                                                                      tAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGaaa                                              <  1:636984/149‑1 (MQ=255)
                                                                                                       aaCGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACGAATATCCTGAAac                                             >  1:143124/1‑149 (MQ=255)
                                                                                                        aCGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAact                                            >  2:875062/1‑149 (MQ=255)
                                                                                                         cGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAactg                                           >  2:1080609/1‑149 (MQ=255)
                                                                                                          gAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAactga                                          >  2:422728/1‑149 (MQ=255)
                                                                                                            gTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAactgact                                        <  1:244128/149‑1 (MQ=255)
                                                                                                            gTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAactgact                                        >  1:553227/1‑149 (MQ=255)
                                                                                                              aaTCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAGTACTCATAATCAATGTAGAGATGACTAATATCCTGAAActgactga                                      >  1:307553/1‑149 (MQ=255)
                                                                                                               aTCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGaa                                     >  1:174053/1‑149 (MQ=255)
                                                                                                               aTCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGaa                                     <  2:18561/149‑1 (MQ=255)
                                                                                                               aTCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTAAGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGTAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGaa                                     <  2:875016/149‑1 (MQ=255)
                                                                                                                tCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAAc                                    >  1:769532/1‑149 (MQ=255)
                                                                                                                tCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAAc                                    >  2:1110496/1‑149 (MQ=255)
                                                                                                                   ttAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGGACTaa                                 >  2:1112215/1‑149 (MQ=255)
                                                                                                                     aaTCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAAtt                               <  2:997286/149‑1 (MQ=255)
                                                                                                                          gCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTc                          <  2:709328/149‑1 (MQ=255)
                                                                                                                          gCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACAGACTGAACTAATTGAGTc                          >  2:877236/1‑149 (MQ=255)
                                                                                                                              aTAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAAc                      >  1:168696/1‑149 (MQ=255)
                                                                                                                               tAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACt                     <  2:201029/149‑1 (MQ=255)
                                                                                                                                 ataaaaaaaCCCGCCGAAGCGGGTTTTATTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAATGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCg                   >  2:681713/5‑149 (MQ=255)
                                                                                                                                 aCGTAAAAACCCGCTTCGGCGTGTTTTTATATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCg                   <  1:941427/149‑1 (MQ=255)
                                                                                                                                  cGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGATTCAAACTCgg                  >  1:293315/1‑149 (MQ=255)
                                                                                                                                  cGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCgg                  <  2:736254/149‑1 (MQ=255)
                                                                                                                                   gTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCTATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGc                 >  2:983589/1‑149 (MQ=255)
                                                                                                                                     aaaaaCCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCaa               >  1:157121/1‑149 (MQ=255)
                                                                                                                                      aaaaCCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGCATACGCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCCAACTCGGCAAg              >  1:1056447/1‑149 (MQ=255)
                                                                                                                                      aaaaCCCGCCGAAGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAg              >  1:595153/1‑149 (MQ=255)
                                                                                                                                      aaaaCCCGCATCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAAATTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAg              <  2:1069988/149‑1 (MQ=255)
                                                                                                                                       aaaCCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAACCTCGGCAAgg             >  2:154148/1‑149 (MQ=255)
                                                                                                                                            cGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATAATCATAATCAATGTAGAGATGACTAATATACTGAAACTGACTGAACTAATTGAGTCAAACACGGCAAGGATTct        >  1:780059/1‑148 (MQ=255)
                                                                                                                                                  ggcggGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATACTGAAACTGACTGAAATAATTGAGTCAAACTCGGCAAGGATTCGATACTa  >  1:1170160/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGATTCGATACTA  >  NC_000913/4189462‑4189756

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