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 system EIIC permease component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,529,5020TC 6.0% 339.9 / 9.4 117T248A (ACC→GCC) sgcCputative PTS system EIIC permease component
Reads supporting (aligned to +/- strand):  ref base T (60/49);  new base C (2/5);  total (63/54)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

Only 100 of 120 total aligned reads displayed.

CAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGCGTCCCTTCATGAAAG  >  NC_000913/4529353‑4529651
                                                                                                                                                     |                                                                                                                                                     
cAGCGTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGt                                                                                                                                                       <  2:72150/150‑1 (MQ=255)
  gCGTCACAGCAGTATCCAGACCGATATACACTTCACGACCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCATCAGTCCATCAAAAAACAGGCGGATAAAGCCCGGGAACAACACCATGATTGCCGAAACGGaaa                                                                                                                                                     >  2:1619382/1‑150 (MQ=255)
    gTCACAGCAGTATCCAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAatc                                                                                                                                                   >  2:1030374/1‑150 (MQ=255)
     gtacaGCAGGATCCAGACCTATATACAATTCAAGTCCTTTGAAAACTTACTGGAAGCACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCATCAATCAGGAGTATCATGCTCTGGAACAACAGCATGATTACCGTTGCGGCAatca                                                                                                                                                  <  2:293238/148‑1 (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:1308587/1‑150 (MQ=255)
                 ccAGACCGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGgcgc                                                                                                                                      >  1:803649/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                                                                                                                                                <  2:1341984/137‑2 (MQ=255)
                      ccGATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAAAAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCt                                                                                                                                 >  1:555644/1‑150 (MQ=255)
                        gATATACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCttt                                                                                                                               <  2:803649/150‑1 (MQ=255)
                         atatACACTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTa                                                                                                                              <  2:156110/150‑1 (MQ=255)
                             acacTTCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGc                                                                                                                          >  1:554604/1‑150 (MQ=255)
                                 ttCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTc                                                                                                                      <  1:790748/150‑1 (MQ=255)
                                  tCACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCa                                                                                                                     <  2:301513/150‑1 (MQ=255)
                                   cACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGCTCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCAc                                                                                                                    >  1:351021/1‑150 (MQ=255)
                                   cACGTCCTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCAc                                                                                                                    >  1:397835/1‑150 (MQ=255)
                                        ccTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAgcc                                                                                                               <  1:57826/150‑1 (MQ=255)
                                        ccTTTGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAgcc                                                                                                               >  2:916120/1‑150 (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)
                                            tGAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCa                                                                                                           <  2:172273/150‑1 (MQ=255)
                                             gAAATATTTCTGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAg                                                                                                          <  1:26370/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)
                                                       tGGAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCCGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGGAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATc                                                                                                >  1:1545921/1‑150 (MQ=255)
                                                        ggAAGAACTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCa                                                                                               >  1:1396150/1‑150 (MQ=255)
                                                         gaagaaCTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAg                                                                                              <  2:809135/150‑1 (MQ=255)
                                                             aaCTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCg                                                                                          <  2:1606937/150‑1 (MQ=255)
                                                              aCTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGa                                                                                         >  1:467943/1‑150 (MQ=255)
                                                               cTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACACCACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCGACCTTTAACGCCGTCACCCGCCGCCAGGCCCAAACTCAGCCCGAg                                                                                        >  2:225222/1‑150 (MQ=255)
                                                               cTTTCTTGCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAg                                                                                        <  1:95478/150‑1 (MQ=255)
                                                                      gCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACg                                                                                 >  1:312139/1‑150 (MQ=255)
                                                                      gCACCGTCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACg                                                                                 >  1:1638827/1‑150 (MQ=255)
                                                                            tCGGAGATGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATAGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATa                                                                           <  2:66277/150‑1 (MQ=255)
                                                                               gagaTGGGCAGTAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCCGACTGGCGCCACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATc                                                                        >  2:1304894/1‑150 (MQ=255)
                                                                                 gaTGGGCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAAACCGAGCACTACGCCGATAATCAc                                                                      >  2:1035459/1‑150 (MQ=255)
                                                                                     ggCAGCAGGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGa                                                                  <  1:446834/150‑1 (MQ=255)
                                                                                           agGCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCACCCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAACGCAGACCGCGCACTTCGCCGATAATCACCGGATCGCCg                                                            >  1:868564/1‑150 (MQ=255)
                                                                                            ggCCCTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGCCTGCCGCAACATTTAACGCCTGCCCCAGCCGCCAGGCCAAAATTCAGCCCGAGCAGTACGCCGATAATCACCGGAGCCCCgc                                                           >  1:678135/1‑149 (MQ=255)
                                                                                              cccTCAACAATCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGAcc                                                         >  1:1268669/1‑150 (MQ=255)
                                                                                                  caacaaTCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGc                                                     <  1:1312838/150‑1 (MQ=255)
                                                                                                   aacaaTCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGcc                                                    <  2:1250221/150‑1 (MQ=255)
                                                                                                     caaTCAGGCGGATCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGt                                                  >  1:1197367/1‑150 (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)
                                                                                                                 tCATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGa                                      <  1:1036558/150‑1 (MQ=255)
                                                                                                                 gcATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGa                                      <  2:1057199/149‑1 (MQ=255)
                                                                                                                  cATGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGAt                                     <  2:1183290/150‑1 (MQ=255)
                                                                                                                   aTGCGCGGGAACAACACCATGATTGCCGCAGCGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAAGCACCGGATCGCCGCCCATGCCGTAGCGGTTCGGGAtc                                    >  2:855310/1‑150 (MQ=255)
                                                                                                                    tGCGCGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGAtct                                   <  1:654304/150‑1 (MQ=255)
                                                                                                                     gcgcGGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGCGCACTACGCCGATAATCCCCGGAGCGCCGACCCTGCCGTAGCGTTTCTGGAgctc                                  >  2:525068/1‑150 (MQ=255)
                                                                                                                      cgcgGGAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGCCCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTcc                                 >  2:185333/1‑150 (MQ=255)
                                                                                                                          ggAACAACACCATGATTGCCGCTACGGTTATCATCAGACTGTCGCAACCTTTATAGCCTTCAGCAGCCGCCAGGCCAAAAATCAGACCGAGCTCTACGCCGATAATCTCCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAg                             <  1:1619382/150‑1 (MQ=255)
                                                                                                                          ggAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGCCCATGCCGTAGCGTTTCTGGATCTCCTGAg                             >  1:1183348/1‑150 (MQ=255)
                                                                                                                          ggAACAACACCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAg                             >  2:247249/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)
                                                                                                                               aacaCCATGATTGCCGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCg                        >  1:1120934/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)
                                                                                                                                      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)
                                                                                                                                            ccGCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGCGTCCCt           >  1:362130/1‑150 (MQ=255)
                                                                                                                                              gCAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTATGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGCGTcc             >  2:1107100/1‑146 (MQ=255)
                                                                                                                                               cAACGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGCGTCCCTTCa        <  2:467943/150‑1 (MQ=255)
                                                                                                                                                  cGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCTTCGATATTGCGTCCCTTCATGa     <  2:397835/150‑1 (MQ=255)
                                                                                                                                                  cGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGCGTCCCTTCATGa     <  2:351021/150‑1 (MQ=255)
                                                                                                                                                  cGGTAATCATCAGACTGGCGCAACCTTTAAAGCCTTCACCAGCCGCCAGGCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGCGTCCCTTCATGa     >  1:268017/1‑150 (MQ=255)
                                                                                                                                                     tAATCATCAGACTGGCTCAACCTTTAAAGCCTTCACCAGCCGCCAGTCCAAAAATCAGACCGAGCACTACGCCGATAATCACCGGATCGCCGACCATGCCGTAGCGTTTCTGGATCTCCTGAGCATCGATATTGCGTCCCTTCATGAAAg  <  1:949132/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 < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: