Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,529,492 G→A 5.9% A251V (GCA→GTA)  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,4920GA 5.9% 318.9 / 6.2 101A251V (GCA→GTA) sgcCputative PTS enzyme IIC component SgcC
Reads supporting (aligned to +/- strand):  ref base G (48/47);  new base A (2/4);  total (50/51)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

Only 100 of 101 total aligned reads displayed.

GATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGCGTC  >  NC_000913/4529346‑4529639
                                                                                                                                                  |                                                                                                                                                   
gATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCg                                                                                                                                                  <  2:810283/150‑1 (MQ=255)
   gCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCaa                                                                                                                                               >  1:759625/1‑150 (MQ=255)
   gCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCaa                                                                                                                                               >  2:198815/1‑150 (MQ=255)
    ccccAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAAc                                                                                                                                              >  2:937216/1‑150 (MQ=255)
    ccccAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAAc                                                                                                                                              >  2:628733/1‑150 (MQ=255)
      ccAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACgg                                                                                                                                            <  2:1537788/150‑1 (MQ=255)
       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                                                                   <  1:631238/150‑1 (MQ=255)
                                                                               aCCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGcc                                                                   <  2:371529/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:302882/150‑1 (MQ=255)
                                                                                                      ccTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCa                                            <  1:937216/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                                                         >  1:1234627/1‑129 (MQ=255)
                                                                                                              aTCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCAcc                                                         <  2:1234627/129‑1 (MQ=255)
                                                                                                              aTCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAg                                    >  2:592060/1‑150 (MQ=255)
                                                                                                                cAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCg                                  <  1:423830/150‑1 (MQ=255)
                                                                                                                      gATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTg                            <  1:414756/150‑1 (MQ=255)
                                                                                                                      gATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTg                            <  1:1146061/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                                                              <  1:1261889/94‑1 (MQ=38)
                                                                                                                                        tctgaTGATTACCGTTGCGGCAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATaa                                                              >  2:1261889/5‑98 (MQ=38)
                                                                                                                                        tctgaTGATTACCGTTGCGGCAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGat          >  2:859028/5‑150 (MQ=255)
                                                                                                                                        cacCATGATTGACGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCCCCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGCCCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGat          >  2:166723/1‑150 (MQ=255)
                                                                                                                                            aTGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTg      <  1:1272875/150‑1 (MQ=255)
                                                                                                                                            aTGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTg      <  2:78924/150‑1 (MQ=255)
                                                                                                                                            aTGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTg      <  2:1232645/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)
                                                                                                                                                  |                                                                                                                                                   
GATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGCGTC  >  NC_000913/4529346‑4529639

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: