Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,143,217 A→C M230L (ATG→CTG)  frwC → putative enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,143,2170AC98.7% 267.9 / ‑6.3 75M230L (ATG→CTG) frwCputative enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (45/29);  minor base T (0/1);  total (45/30)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GGTGATCATGGGTCTGATGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTC  >  NC_000913/4143072‑4143362
                                                                                                                                                 |                                                                                                                                                 
ggTGATCATGGGTCTGATGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCgctggggctggg                                                                                                                                               <  1:814419/150‑1 (MQ=255)
       aTGGGTCTGATGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCg                                                                                                                                        >  2:1102413/1‑150 (MQ=255)
                aTGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTg                                                                                                                              >  2:45738/1‑151 (MQ=255)
                aTGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTg                                                                                                                              >  2:231567/1‑151 (MQ=255)
                aTGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTg                                                                                                                              >  1:448501/1‑151 (MQ=255)
                aTGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGATGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTg                                                                                                                              >  1:1014900/1‑151 (MQ=255)
                       cGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTaa                                                                                                                       <  1:1084542/151‑1 (MQ=255)
                       cGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTa                                                                                                                        >  1:269230/1‑150 (MQ=255)
                        gTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTaa                                                                                                                       >  2:477808/1‑150 (MQ=255)
                         tccGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGTTTTGCGTTGCTCAGGGTGTTTATGCCGTGGTGGCCATCGCTGCCGTGGGTATTTTTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTaaa                                                                                                                      <  1:376400/148‑1 (MQ=255)
                           cGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAAt                                                                                                                   >  2:453912/1‑151 (MQ=255)
                             atatGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAAttt                                                                                                                 <  2:685892/151‑1 (MQ=255)
                             atatGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAAttt                                                                                                                 >  2:699405/1‑151 (MQ=255)
                              tatGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTtc                                                                                                                >  1:1011723/1‑151 (MQ=255)
                                tGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTtctc                                                                                                              >  1:77657/1‑151 (MQ=255)
                                  ggCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCg                                                                                                            <  1:419494/151‑1 (MQ=255)
                                  ggCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCg                                                                                                            <  1:1210079/151‑1 (MQ=255)
                                    cGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCa                                                                                                          <  2:95735/151‑1 (MQ=255)
                                     ggTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGc                                                                                                           >  1:510724/1‑149 (MQ=255)
                                            ttAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAgc                                                                                                   <  2:1053163/150‑1 (MQ=255)
                                             tAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAgcg                                                                                                  >  2:144578/1‑150 (MQ=255)
                                                cAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGa                                                                                               >  2:321728/1‑150 (MQ=255)
                                                 aaacgggCCTATGGCTTCGTGCTGTTTTGCGTTCCTGGGGGTGTTCATCGCAGGGTGCACATCGCTGGTGTTGGTATTTGTATCCCGCCGCTGGGGTTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGaa                                                                                              <  2:1029664/145‑1 (MQ=255)
                                                  aaGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAAc                                                                                            >  1:102126/1‑151 (MQ=255)
                                                        ccTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCa                                                                                       >  2:857736/1‑150 (MQ=255)
                                                          tATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCaa                                                                                      >  1:878255/1‑149 (MQ=255)
                                                           aTGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCaaa                                                                                     <  1:1023340/149‑1 (MQ=255)
                                                            tGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGc                                                                                   <  2:596856/150‑1 (MQ=255)
                                                              ccTTCATGCTGATTTGCGTTGCTCAGGGTTTTTATACCGTGGTGGCCTTCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGc                                                                                <  1:707645/151‑1 (MQ=255)
                                                              ccTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGc                                                                                >  2:1124808/1‑151 (MQ=255)
                                                               cTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGc                                                                                >  1:647282/1‑150 (MQ=255)
                                                                ttCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCAc                                                                              >  2:1149444/1‑151 (MQ=255)
                                                                   aTGCTGATTTGCGTTGCTCAGGTTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTACCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACt                                                                             <  2:534241/149‑1 (MQ=255)
                                                                   aTGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTgg                                                                           <  1:351173/151‑1 (MQ=255)
                                                                    tGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTgg                                                                           >  2:400816/1‑150 (MQ=255)
                                                                     gCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTg                                                                         <  2:879141/151‑1 (MQ=255)
                                                                      cTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTg                                                                         >  2:587192/1‑150 (MQ=255)
                                                                        gATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTTTCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGa                                                                        <  2:503417/149‑1 (MQ=255)
                                                                        gATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATg                                                                      <  1:455664/151‑1 (MQ=255)
                                                                        gATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATg                                                                      >  2:182760/1‑151 (MQ=255)
                                                                              cGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCtg                                                                 >  1:865072/1‑150 (MQ=255)
                                                                              cGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCtg                                                                 >  2:1007738/1‑150 (MQ=255)
                                                                              cGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTACCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCtgt                                                                <  2:257620/151‑1 (MQ=255)
                                                                                ttGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCtgt                                                                <  1:1082607/149‑1 (MQ=255)
                                                                                ttGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTgtgt                                                              >  2:1039817/1‑151 (MQ=255)
                                                                                      aGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGt                                                        >  1:234846/1‑151 (MQ=255)
                                                                                       gggTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGtt                                                       >  2:779409/1‑151 (MQ=255)
                                                                                       gggTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGtt                                                       <  1:340301/151‑1 (MQ=255)
                                                                                              tataCCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGa                                                  >  2:1181489/1‑149 (MQ=255)
                                                                                              tataCCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGa                                                  >  1:1128750/1‑149 (MQ=255)
                                                                                               ataCCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAgg                                               >  2:533700/1‑151 (MQ=255)
                                                                                               ataCCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAg                                                <  2:462516/150‑1 (MQ=255)
                                                                                                 aCCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAgg                                               >  2:916258/1‑149 (MQ=255)
                                                                                                  ccGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGc                                            >  2:878219/1‑151 (MQ=255)
                                                                                                    gtggtgGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCg                                           >  2:866491/1‑150 (MQ=255)
                                                                                                          gCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCtt                                    <  2:1083205/151‑1 (MQ=255)
                                                                                                           ccATCGCTGCCGTGGGTTTTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCttt                                   <  1:147630/151‑1 (MQ=255)
                                                                                                           ccATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCttt                                   <  1:788535/151‑1 (MQ=255)
                                                                                                            cATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTc                                  <  2:1227197/151‑1 (MQ=255)
                                                                                                            cATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTc                                  <  2:843031/151‑1 (MQ=255)
                                                                                                               cGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTcgc                                <  2:277183/150‑1 (MQ=255)
                                                                                                                 cTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTcgccgc                             >  2:366703/1‑151 (MQ=255)
                                                                                                                 cTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTcgccgc                             >  1:150965/1‑151 (MQ=255)
                                                                                                                  tGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCt                            >  2:655928/1‑151 (MQ=255)
                                                                                                                      gTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCg                        >  2:984795/1‑151 (MQ=255)
                                                                                                                      gTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCg                        >  2:157104/1‑151 (MQ=255)
                                                                                                                      gTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATAGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGcc                         <  2:547677/150‑1 (MQ=255)
                                                                                                                       tGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTcgccgc                             <  1:243037/145‑1 (MQ=255)
                                                                                                                           tATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCg                   >  2:186515/1‑151 (MQ=255)
                                                                                                                            aTTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATcc                    <  1:321929/149‑1 (MQ=255)
                                                                                                                            aTTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGc                  >  1:7801/1‑151 (MQ=255)
                                                                                                                              ttGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCt                 <  1:505548/150‑1 (MQ=255)
                                                                                                                                 tATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCg              >  2:25353/1‑150 (MQ=255)
                                                                                                                                  aTCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCgtg            >  1:1084083/1‑151 (MQ=255)
                                                                                                                                     ccgccgCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTc          >  2:852796/1‑150 (MQ=255)
                                                                                                                                       gccgcTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCAtt       >  1:526779/1‑151 (MQ=255)
                                                                                                                                          gctggggctgggGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCTGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCAttcct    <  1:744247/151‑1 (MQ=255)
                                                                                                                                          gctggggctgggGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCAttcct    >  1:777367/1‑151 (MQ=255)
                                                                                                                                             ggggctggggctGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCAttccttc  <  1:453900/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGTGATCATGGGTCTGATGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTC  >  NC_000913/4143072‑4143362

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: