Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,529,492 G→A 6.7% A251V (GCA→GTA)  sgcC ← putative PTS system EIIC permease component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,529,4920GA 6.7% 372.1 / 9.3 118A251V (GCA→GTA) sgcCputative PTS system EIIC permease component
Reads supporting (aligned to +/- strand):  ref base G (62/48);  new base A (3/5);  total (65/53)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

Only 100 of 119 total aligned reads displayed.

GGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGCGTCCC  >  NC_000913/4529345‑4529641
                                                                                                                                                   |                                                                                                                                                     
ggATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGcc                                                                                                                                                     <  2:1511087/150‑1 (MQ=255)
ggATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGcc                                                                                                                                                     <  2:1517670/150‑1 (MQ=255)
 gATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAAAACCATGATTGCCg                                                                                                                                                    >  1:1057199/1‑150 (MQ=255)
     ccccAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAAc                                                                                                                                                >  2:1017917/1‑150 (MQ=255)
       ccAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACgg                                                                                                                                              >  2:1312838/1‑150 (MQ=255)
        cAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGt                                                                                                                                             <  2:72150/150‑1 (MQ=255)
        cAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTACCgttgcggc                                                                                                                                             >  2:1107876/1‑143 (MQ=255)
          gCGTCACAGCAGTATCCAGACCGATATACACTTCACGACCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCATCAGTCCATCAAAAAACAGGCGGATAAAGCCCGGGAACAACACCATGATTGCCGAAACGGaaa                                                                                                                                           >  2:1619382/1‑150 (MQ=255)
              cacaGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGACATTTCTTGCACCGACGTAGATGGGCAGCAGGCCCTCAACAAAAAGCCGTATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAatcat                                                                                                                                       >  1:397999/1‑150 (MQ=255)
              cacaGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAatcat                                                                                                                                       >  2:790748/1‑150 (MQ=255)
                 agcagTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTACCGTTGCGGCAatcatggt                                                                                                                                    <  1:589691/150‑4 (MQ=255)
                    agTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTACCGTTGCGGCAatcatggtgtt                                                                                                                                 <  2:1591555/150‑7 (MQ=255)
                     ttATCCATACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTg                                                                                                                                <  2:743204/149‑1 (MQ=255)
                      tATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTgg                                                                                                                               <  2:1222069/150‑1 (MQ=255)
                       aTCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCACCGGTAATCATCAGACTGgc                                                                                                                              >  2:57826/1‑150 (MQ=255)
                        tCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGgcg                                                                                                                             <  1:267212/150‑1 (MQ=255)
                        accAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGATATGGGCATCAGGCCCGCATCAATCAGGCGGATCATGCGCGGGAACAACTCCATTATTGCCGCAACGGTAATCATCAGACTGgcg                                                                                                                             <  1:382237/149‑1 (MQ=255)
                         ccAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGgcgc                                                                                                                            >  1:803649/1‑150 (MQ=255)
                         ccAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGgcgc                                                                                                                            >  1:1308587/1‑150 (MQ=255)
                         ccAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGgcgc                                                                                                                            >  1:466732/1‑150 (MQ=255)
                            gACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAAc                                                                                                                         >  2:654304/1‑150 (MQ=255)
                            gACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAAc                                                                                                                         >  2:949132/1‑150 (MQ=255)
                            gACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTACCGTTGCGGCAatcatg                                                                                                                                      >  1:1341984/1‑136 (MQ=255)
                              ccGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAAAAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCt                                                                                                                       >  1:555644/1‑150 (MQ=255)
                                gATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCttt                                                                                                                     <  2:803649/150‑1 (MQ=255)
                                     acacTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGc                                                                                                                >  1:554604/1‑150 (MQ=255)
                                         ttCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTc                                                                                                            <  1:790748/150‑1 (MQ=255)
                                         ttCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTAAAGCCTTCa                                                                                                           >  1:182603/1‑150 (MQ=255)
                                          tCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCa                                                                                                           <  2:301513/150‑1 (MQ=255)
                                           cACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGCTCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCAc                                                                                                          >  1:351021/1‑150 (MQ=255)
                                           cACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCAc                                                                                                          >  1:397835/1‑150 (MQ=255)
                                                ccTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAgcc                                                                                                     >  2:916120/1‑150 (MQ=255)
                                                ccTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAgcc                                                                                                     >  1:566015/1‑150 (MQ=255)
                                                ccTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAgcc                                                                                                     <  1:57826/150‑1 (MQ=255)
                                                 cTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAACGCCTTCACCAgccg                                                                                                    >  1:1015372/1‑150 (MQ=255)
                                                 cTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAgccg                                                                                                    >  2:446834/1‑150 (MQ=255)
                                                 cTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAgccg                                                                                                    <  1:1030374/150‑1 (MQ=255)
                                                  tttGAAATAATTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGAAGCAGGCCCTCAACACTCAGGCGGATCAAGCGCGGGAACCAAACCATGATTGCCGCAACGGTAATAATCAGACTGCCGCAACCTTTAAAGCCTTCACCAgcctc                                                                                                   >  1:1533702/1‑148 (MQ=255)
                                                   ttGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAgccgcc                                                                                                  >  2:1036558/1‑150 (MQ=255)
                                                    tGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCa                                                                                                 <  2:172273/150‑1 (MQ=255)
                                                      aaaTATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAgg                                                                                               <  2:566015/150‑1 (MQ=255)
                                                      aaaTATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAgg                                                                                               <  1:567357/150‑1 (MQ=255)
                                                          atTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCaa                                                                                           >  1:1678279/1‑150 (MQ=255)
                                                            ttCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCaaaa                                                                                         <  1:1017917/150‑1 (MQ=255)
                                                            atCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCATCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTAGCGCAACCTTAAAAGCCTTAAACAGCCGCCAGGCCaaaa                                                                                         >  1:413698/2‑150 (MQ=255)
                                                               tGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCCGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGGAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATc                                                                                      >  1:1545921/1‑150 (MQ=255)
                                                               tGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATc                                                                                      <  2:359546/150‑1 (MQ=255)
                                                                ggAAGAACTTTCTTGCACAGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCa                                                                                     >  1:764851/1‑150 (MQ=255)
                                                                 gaagaaCTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAg                                                                                    <  2:809135/150‑1 (MQ=255)
                                                                     aaCTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCg                                                                                <  2:1234071/150‑1 (MQ=255)
                                                                     aaCTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCg                                                                                <  2:1606937/150‑1 (MQ=255)
                                                                       cTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACACCACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCGACCTTTAACGCCGTCACCCGCCGCCAGGCCCAAACTCAGCCCGAg                                                                              >  2:225222/1‑150 (MQ=255)
                                                                       cTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAg                                                                              <  1:95478/150‑1 (MQ=255)
                                                                              gCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACg                                                                       >  1:1638827/1‑150 (MQ=255)
                                                                              gCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACg                                                                       >  1:312139/1‑150 (MQ=255)
                                                                                   gTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGAt                                                                  >  2:1421852/1‑150 (MQ=255)
                                                                                    tCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATAGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATa                                                                 <  2:66277/150‑1 (MQ=255)
                                                                                     cGGAGCTGGGCAGCAGGCCCTCAACAATCAGGCGGATCATTCGCTGGAACAACACCATGATTGACTCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATaa                                                                <  2:253333/150‑1 (MQ=255)
                                                                                     cGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGCCCGAGCACTACGCCGATaa                                                                >  1:1612219/1‑150 (MQ=255)
                                                                                       gagaTGGGCAGTAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCCGACTGGCGCCACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATc                                                              >  2:1304894/1‑150 (MQ=255)
                                                                                         gaTGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAAACCGAGCACTACGCCGATAATCAc                                                            >  2:1035459/1‑150 (MQ=255)
                                                                                             ggCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGa                                                        <  1:446834/150‑1 (MQ=255)
                                                                                                    ggCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGCCTGCCGCAACATTTAACGCCTGCCCCAGCCGCCAGGCCAAAATTCAGCCCGAGCAGTACGCCGATAATCACCGGAGCCCCgc                                                 >  1:678135/1‑149 (MQ=255)
                                                                                                     gCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGAc                                                >  2:1710794/1‑150 (MQ=255)
                                                                                                      cccTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGAcc                                               >  1:1268669/1‑150 (MQ=255)
                                                                                                          caacaaTCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGc                                           <  1:1312838/150‑1 (MQ=255)
                                                                                                           aacaaTCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGcc                                          <  2:1250221/150‑1 (MQ=255)
                                                                                                             caaTCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGt                                        <  2:397999/150‑1 (MQ=255)
                                                                                                              aaTCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTa                                       <  2:312139/150‑1 (MQ=255)
                                                                                                               aTCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAg                                      >  1:1183290/1‑150 (MQ=255)
                                                                                                                tCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGc                                     <  2:764851/150‑1 (MQ=255)
                                                                                                                tCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCATTCACAAGCCGACAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGACGTAGc                                     >  1:1045884/1‑150 (MQ=255)
                                                                                                                   ggCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCTGTAGCGtt                                  >  1:1583757/1‑150 (MQ=255)
                                                                                                                   ggCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGtt                                  >  2:936516/1‑150 (MQ=255)
                                                                                                                     cGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTc                                <  1:1670225/150‑1 (MQ=255)
                                                                                                                     cGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTc                                >  2:254677/1‑150 (MQ=255)
                                                                                                                      ggATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCt                               <  2:1352163/150‑1 (MQ=255)
                                                                                                                      ggAGCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACGATGCCGTTGAGtttgt                               <  2:1533702/150‑3 (MQ=255)
                                                                                                                        aTCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCTCCGGATCGCCGACCATGCCGTAGCGTTTCTgg                             <  1:321284/150‑1 (MQ=255)
                                                                                                                         tCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGa                            <  1:1036558/150‑1 (MQ=255)
                                                                                                                         gcATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGa                            <  2:1057199/149‑1 (MQ=255)
                                                                                                                           aTGCGCGGGAACAACACCATGATTGCCGCAGCGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAAGCACCGGATCGCCGCCCATGCCGTAGCGGTTCGGGAtc                          >  2:855310/1‑150 (MQ=255)
                                                                                                                             gcgcGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGCGCACTACGCCGATAATCCCCGGAGCGCCGACCCTGCCGTAGCGTTTCTGGAgctc                        >  2:525068/1‑150 (MQ=255)
                                                                                                                                  ggAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGCCCATGCCGTAGCGTTTCTGGATCTCCTGAg                   >  1:1183348/1‑150 (MQ=255)
                                                                                                                                  ggAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAg                   >  2:247249/1‑150 (MQ=255)
                                                                                                                                  ggAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCAGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAg                   >  1:1189562/1‑150 (MQ=255)
                                                                                                                                  ggAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCAGCCAGGCCAAAAATCAGACCGAGCACTACGCCTATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAg                   <  2:555644/150‑1 (MQ=255)
                                                                                                                                     acaacaCCATGATTGCCGCAACGGTAATCATCTGACTGGAGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCAt                >  2:212004/1‑150 (MQ=255)
                                                                                                                                       aacaCCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCg              >  1:1618788/1‑150 (MQ=255)
                                                                                                                                        acacCATGATTGCCGCAACGGTAATCATCAGAAAGGCGCAACCTTTAAAGCCTTCAAAATCCAAAAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGa             >  2:1074414/1‑150 (MQ=255)
                                                                                                                                         cacCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGCCCATGCCGTAGCGTTTCTGGCTCTCCTGAGCATCGat            >  1:1129358/1‑150 (MQ=255)
                                                                                                                                         cacCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGat            >  1:1135288/1‑150 (MQ=255)
                                                                                                                                            cATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATAtt         <  2:1308587/150‑1 (MQ=255)
                                                                                                                                             aTGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTg        <  1:254677/150‑1 (MQ=255)
                                                                                                                                             aTGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTg        <  1:1220031/150‑1 (MQ=255)
                                                                                                                                              tGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATCATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGc       <  1:916120/150‑1 (MQ=255)
                                                                                                                                              tGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGc       >  1:888094/1‑150 (MQ=255)
                                                                                                                                               gATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTATAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGCg      <  1:1425478/150‑1 (MQ=255)
                                                                                                                                                  tGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGAGCTCCTGAGCATCGATATTGCGTcc   <  1:936516/150‑1 (MQ=255)
                                                                                                                                                   gCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCCGGCCAAAACTCAGACCGAGAACTACGCCGATAATCACCGGATCGCCGAGCATGCGGGAGCGTTTCTGGCTCTCCTGAGCATCGAGATTGCGTccc  >  1:80836/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGATGCCCCAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGCGTCCC  >  NC_000913/4529345‑4529641

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: