Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,189,605 T→C 5.1% intergenic (+32/‑181) 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,6050TC 5.1% 268.0 / 4.7 99intergenic (+32/‑181)rpoC/yjaZRNA polymerase, beta prime subunit/stationary phase growth adaptation protein
Reads supporting (aligned to +/- strand):  ref base T (51/41);  new base C (2/3);  total (54/45)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

GCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGATTCG  >  NC_000913/4189459‑4189750
                                                                                                                                                  |                                                                                                                                                 
gCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGAAGGTCTGGGCGGTTCTGATAAAGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTc                                                                                                                                                 >  2:1083738/1‑149 (MQ=255)
   tCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTcggc                                                                                                                                              <  2:109128/149‑1 (MQ=255)
   tCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTcggc                                                                                                                                              <  2:507761/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                                                                                                                                <  1:972422/149‑1 (MQ=255)
                 gTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATgg                                                                                                                                <  2:616550/149‑1 (MQ=255)
                   gAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATgggg                                                                                                                              <  1:166131/149‑1 (MQ=255)
                   gAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATgggg                                                                                                                              >  2:18640/1‑149 (MQ=255)
                   gAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATgggg                                                                                                                              >  1:347255/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)
                                                               cAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTAT‑GGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCt                                                                                 >  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:363865/149‑1 (MQ=255)
                                                                               ctgAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCa                                                                  >  1:343097/1‑149 (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:525945/149‑1 (MQ=255)
                                                                                    cGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTa                                                             <  2:198004/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                                         >  2:686241/1‑149 (MQ=255)
                                                                                                        aTAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGaa                                         <  1:909295/149‑1 (MQ=255)
                                                                                                        aTAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGaa                                         >  2:665374/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)
                                                                                                                                                  |                                                                                                                                                 
GCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGATTCG  >  NC_000913/4189459‑4189750

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: ‑