Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,529,502 T→C 6.0% T248A (ACC→GCC)  sgcC ← putative PTS enzyme IIC component SgcC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,529,5020TC 6.0% 302.9 / 7.0 99T248A (ACC→GCC) sgcCputative PTS enzyme IIC component SgcC
Reads supporting (aligned to +/- strand):  ref base T (44/49);  new base C (2/4);  total (46/53)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 101 total aligned reads displayed.

CAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGCGTCCCTTCATGAAAG  >  NC_000913/4529353‑4529651
                                                                                                                                                     |                                                                                                                                                     
cAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGt                                                                                                                                                       >  2:271989/1‑150 (MQ=255)
   cGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAat                                                                                                                                                    >  2:414756/1‑150 (MQ=255)
    gTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAatc                                                                                                                                                   >  1:826592/1‑150 (MQ=255)
    gTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAatc                                                                                                                                                   >  1:421806/1‑150 (MQ=255)
      cacaGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAatcat                                                                                                                                                 >  2:631238/1‑150 (MQ=255)
       acaGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAatcatc                                                                                                                                                <  2:686919/150‑1 (MQ=255)
       acaGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAatcatc                                                                                                                                                <  2:508286/150‑1 (MQ=255)
        cagcagTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAAtcatca                                                                                                                                               >  2:673560/1‑150 (MQ=255)
        cagcagTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAAATTTCTTGAACCGTCGGAGATGGGCAGCAGGACCTCAACAATCAGGCGGATCAAGCGCGGGAACAACACCATGATTGCCGAAAAGGTAAtcatca                                                                                                                                               >  2:471198/1‑150 (MQ=255)
         agcagTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAg                                                                                                                                              <  1:232014/150‑1 (MQ=255)
          gcagTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGa                                                                                                                                             >  2:1378065/1‑150 (MQ=255)
          gcagTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTACCGTTGCGGCAatcatggtg                                                                                                                                             <  1:271810/150‑5 (MQ=255)
           cagTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGAc                                                                                                                                            <  2:1130664/150‑1 (MQ=255)
           cagTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTACCGTTGCGGCAatcatggtgt                                                                                                                                            <  1:836610/150‑6 (MQ=255)
             gTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTg                                                                                                                                          >  2:1389202/1‑150 (MQ=255)
             gTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTACCGTTGCGGCAatcatggtgttg                                                                                                                                          <  2:432998/150‑8 (MQ=255)
              tATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTgg                                                                                                                                         >  2:795892/1‑150 (MQ=255)
                tCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGgcg                                                                                                                                       >  2:1145559/1‑150 (MQ=255)
                 ccAGACCAATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGgcgc                                                                                                                                      >  2:1350014/1‑150 (MQ=255)
                         atatACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTa                                                                                                                              <  1:1380358/150‑1 (MQ=255)
                           ataCACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTTAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTaaa                                                                                                                            <  2:1269052/150‑1 (MQ=255)
                           ataCACTTCACGTCCTCTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTaaa                                                                                                                            >  1:390360/1‑150 (MQ=255)
                            taCACTTCACTACATTTGAAATATTTCTGTAAGAACTTTCTTTCACCGTCGGAGATGGGCAGAAATCCCTAAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAg                                                                                                                           >  2:848798/1‑150 (MQ=255)
                            taCACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAg                                                                                                                           >  1:630820/1‑150 (MQ=255)
                            taCACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAg                                                                                                                           >  2:1520795/1‑150 (MQ=255)
                             acacTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGc                                                                                                                          >  1:1469794/1‑150 (MQ=255)
                               acTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCt                                                                                                                        <  2:376793/150‑1 (MQ=255)
                                  tcccGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACTATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCa                                                                                                                     <  2:1124075/147‑1 (MQ=255)
                                   cACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCAc                                                                                                                    <  1:1145559/150‑1 (MQ=255)
                                        ccTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAgcc                                                                                                               <  2:897131/150‑1 (MQ=255)
                                         cTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAgccg                                                                                                              >  1:371529/1‑150 (MQ=255)
                                           ttGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAgccgcc                                                                                                            >  2:302882/1‑150 (MQ=255)
                                             taaaTATTTCTGGTATAACTTTCTTGCACCGTCGGAGATGGGCAGGAGGCCCTCAACAATCAGGCGGATCATGCGCGGGTACAACACCATGATTGGCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAg                                                                                                          <  2:255197/149‑1 (MQ=255)
                                               aaTATATCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGAAGCAGGCCCTCAACAATAAGGCGGATCATGCGCGCGAAAAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTAAAAGCATACACCAGCCCCCAGGc                                                                                                        >  2:1146061/1‑150 (MQ=255)
                                                    ttCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCaaaa                                                                                                   >  2:412593/1‑150 (MQ=255)
                                                       tGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATc                                                                                                <  2:759625/150‑1 (MQ=255)
                                                         gaagaaCTTTCTTGGACCGTCGGAGATGGGCAGCAGGCCCGCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAg                                                                                              <  2:471489/150‑1 (MQ=255)
                                                         gaagaaCTTTCTTGCACCGTCGGAGATTGGCAGCAGGACCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGAAACTGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAg                                                                                              <  1:471198/150‑1 (MQ=255)
                                                         gaagaaCTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAg                                                                                              >  1:1408116/1‑150 (MQ=255)
                                                          aagaaCTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAgc                                                                                             >  1:1528204/1‑149 (MQ=255)
                                                           agaaCTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAgcc                                                                                            >  1:1141945/1‑148 (MQ=255)
                                                           agaaCTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGAc                                                                                            <  1:1538719/150‑1 (MQ=255)
                                                              aCTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGa                                                                                         >  1:216813/1‑150 (MQ=255)
                                                                tttCTTGCACCGTCGGAGATGGGCAGCAGGACCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGc                                                                                       <  1:608990/150‑1 (MQ=255)
                                                                        aCCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGcc                                                                               <  1:1508400/150‑1 (MQ=255)
                                                                        aCCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGcc                                                                               <  2:371529/150‑1 (MQ=255)
                                                                        aCCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGcc                                                                               <  1:631238/150‑1 (MQ=255)
                                                                             cGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATaa                                                                          <  1:412593/150‑1 (MQ=255)
                                                                                  aTGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCAcc                                                                     <  1:261213/150‑1 (MQ=255)
                                                                                   tGGGCAGCAGGCCCTCAACAATCAGGCGGAGCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCg                                                                    <  2:413840/150‑1 (MQ=255)
                                                                                    gggCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCgg                                                                   >  2:71594/1‑150 (MQ=255)
                                                                                    gggCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCATTAAAGCCTTCACCACCCGCCAGGCCAAAAATCAGACCCAGCACAACGCCGATAAACACCgg                                                                   >  2:130351/1‑150 (MQ=255)
                                                                                         gcagGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGc                                                              >  2:542448/1‑150 (MQ=255)
                                                                                            ggCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGa                                                           >  2:137039/1‑150 (MQ=255)
                                                                                              cccTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGAcc                                                         >  1:65483/1‑150 (MQ=255)
                                                                                              cccTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGAcc                                                         >  1:560446/1‑150 (MQ=255)
                                                                                               ccTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCa                                                        <  1:937216/150‑1 (MQ=255)
                                                                                               ccTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCa                                                        <  1:302882/150‑1 (MQ=255)
                                                                                                 tCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATg                                                      >  1:1278209/1‑150 (MQ=255)
                                                                                                  caacaaTCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGc                                                     >  1:1416301/1‑150 (MQ=255)
                                                                                                    acaaTCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCg                                                   >  1:606284/1‑150 (MQ=255)
                                                                                                    acaaTCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCg                                                   <  1:1378065/150‑1 (MQ=255)
                                                                                                      aaTCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTa                                                 <  1:633753/150‑1 (MQ=255)
                                                                                                      aaTCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTa                                                 >  2:447589/1‑150 (MQ=255)
                                                                                                       aTCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCAcc                                                                     <  2:1234627/129‑1 (MQ=255)
                                                                                                       aTCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCAcc                                                                     >  1:1234627/1‑129 (MQ=255)
                                                                                                       aTCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAg                                                >  2:592060/1‑150 (MQ=255)
                                                                                                         cAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCg                                              <  1:423830/150‑1 (MQ=255)
                                                                                                               gATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTg                                        <  1:1146061/150‑1 (MQ=255)
                                                                                                               gATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTg                                        <  1:414756/150‑1 (MQ=255)
                                                                                                               gATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTg                                        >  1:25761/1‑150 (MQ=255)
                                                                                                                aTCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTgg                                       <  2:216813/150‑1 (MQ=255)
                                                                                                                 tCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGa                                      <  1:1350014/150‑1 (MQ=255)
                                                                                                                 tCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCAGATCGCCGACCATGCCGTAGCGTTTCTGGa                                      >  1:1053428/1‑150 (MQ=255)
                                                                                                                   aTGCGCTGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCTCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGAtc                                    <  2:877940/150‑1 (MQ=255)
                                                                                                                   aTGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGAtc                                    <  2:390360/150‑1 (MQ=255)
                                                                                                                      cgcgGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTcc                                 >  1:78924/1‑150 (MQ=255)
                                                                                                                       gcgGGAACAACACGATGATTGCGGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACGAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATATTCACCGCATCGCCGACCATGCCGTAGCGTTTCTGGATTTCCt                                <  2:301147/150‑1 (MQ=255)
                                                                                                                         gggAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGCGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTgc                              >  1:694249/1‑149 (MQ=255)
                                                                                                                         gggAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGa                              <  1:795892/150‑1 (MQ=255)
                                                                                                                           gAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGc                            <  1:628733/150‑1 (MQ=255)
                                                                                                                              caacaCCATGATTGCCGCAACGGTAATCATCAGCCTGGCGCAGCCTTTAAAGCATGCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCCACGGATCGCCGACCATGGGGTAGCGTGTCTGGGTCTCCTGCGCATc                         >  2:502023/1‑150 (MQ=255)
                                                                                                                              caacaCCATGATTGCCGCAACGGTAATCATCAGACTGGAGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACTCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATc                         <  2:907420/150‑1 (MQ=255)
                                                                                                                               aacaCCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCg                        <  2:183/150‑1 (MQ=255)
                                                                                                                                 tctgaTGATTACCGTTGCGGCAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATaa                                                                          >  2:1261889/5‑98 (MQ=38)
                                                                                                                                 tctgaTGATTACCGTTGCGGCAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATaa                                                                          <  1:1261889/94‑1 (MQ=38)
                                                                                                                                 tctgaTGATTACCGTTGCGGCAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGat                      >  2:859028/5‑150 (MQ=255)
                                                                                                                                 cacCATGATTGACGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCCCCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGCCCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGat                      >  2:166723/1‑150 (MQ=255)
                                                                                                                                     aTGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTg                  <  2:1232645/150‑1 (MQ=255)
                                                                                                                                     aTGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTg                  <  1:1272875/150‑1 (MQ=255)
                                                                                                                                     aTGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTg                  <  2:78924/150‑1 (MQ=255)
                                                                                                                                     aTGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTg                  <  1:1198566/150‑1 (MQ=255)
                                                                                                                                      tGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGc                 <  2:1274129/150‑1 (MQ=255)
                                                                                                                                         ttGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGCGTc              <  1:130351/150‑1 (MQ=255)
                                                                                                                                              gCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGCGTCCCTTc         <  2:826592/150‑1 (MQ=255)
                                                                                                                                                aaCGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCCCCAGCCGCAAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGCGTCCCTTCAt       >  2:1212524/1‑150 (MQ=255)
                                                                                                                                                aaCGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGCGTCCCTTCAt       <  1:848798/150‑1 (MQ=255)
                                                                                                                                                 aCGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATTAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGCGTCCCTTCATg      <  1:1520795/150‑1 (MQ=255)
                                                                                                                                                   ggTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGCGTCCCTTCATGaa    <  2:1469794/150‑1 (MQ=255)
                                                                                                                                                     tAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGCGTCCCTTCATGAAAg  <  1:589115/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGCGTCCCTTCATGAAAG  >  NC_000913/4529353‑4529651

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: