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.5% 235.0 / NA 67M230L (ATG→CTG) frwCputative enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (41/25);  minor base T (1/0);  total (42/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGGTGATCATGGGTCTGATGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTC  >  NC_000913/4143069‑4143362
                                                                                                                                                    |                                                                                                                                                 
ggCGGTGATCATGGGTCTGATGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTg                                                                                                                                                 >  2:792163/1‑151 (MQ=255)
      gATCATGGGTCTGATGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTggggctggggct                                                                                                                                            <  1:405899/150‑1 (MQ=255)
      gATCATGGGTCTGATGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGgggctggggctg                                                                                                                                           <  2:775128/151‑1 (MQ=255)
        tcaggcgTCTGATGCTGGCGTTCGCAATCGGCGGTCGGGTTAACCAAGTGGCCTATGCCTTCATGCTGATTTGCGATGCACAGGGTGTATATACCGTGCTGGCCATCGCTGCCGAGGGTATTTGCATCCCGTCGCtggggctggggc                                                                                                                                             >  2:966104/7‑147 (MQ=255)
         cATGGGTCTGATGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGc                                                                                                                                         <  2:1058045/150‑1 (MQ=255)
         cATGGGTCTGATGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCg                                                                                                                                        <  1:48628/151‑1 (MQ=255)
             ggTCTGATGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGc                                                                                                                                    >  1:126985/1‑151 (MQ=255)
              gTCTGATGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCt                                                                                                                                   <  1:359308/151‑1 (MQ=255)
                  gATGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGa                                                                                                                                 <  1:171240/149‑1 (MQ=255)
                    tGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCGTCGCTGCCGTGGTTATTTGTATCCCGCCGCTGGGGCTGGGGCTTGCGACGCTGAtt                                                                                                                               <  2:958185/149‑1 (MQ=255)
                    tGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTgg                                                                                                                             >  1:871912/1‑151 (MQ=255)
                    tGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTg                                                                                                                              <  1:2385/150‑1 (MQ=255)
                    tGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTg                                                                                                                              <  1:473437/150‑1 (MQ=255)
                      cTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTc                                                                                                                           >  1:1008520/1‑151 (MQ=255)
                           gTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTaa                                                                                                                       >  1:1138222/1‑150 (MQ=255)
                            ttCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTaaa                                                                                                                      <  2:586133/150‑1 (MQ=255)
                               gATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAAtt                                                                                                                  >  1:455677/1‑151 (MQ=255)
                                  atGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTtct                                                                                                               >  1:54297/1‑151 (MQ=255)
                                   tGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAACTTtctc                                                                                                              >  2:588105/1‑151 (MQ=255)
                                   tGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTtctc                                                                                                              >  2:765931/1‑151 (MQ=255)
                                         gTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGc                                                                                                           >  2:836326/1‑148 (MQ=255)
                                            cGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCaga                                                                                                        >  2:8927/1‑148 (MQ=255)
                                                  aCAAAGTGGCCTATGCCTTCATGCTGATTTGCGCTGCTCAGGGAGTTTATAGCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGTTGGGCGTGGCTACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGcgcg                                                                                                >  2:776041/1‑150 (MQ=255)
                                                   cAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATCCCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGaa                                                                                              <  2:167229/151‑1 (MQ=255)
                                                       gTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCg                                                                                          >  1:68567/1‑151 (MQ=255)
                                                               tGCCTTCATGCTGATTTGCGTGGCTCGGGCTGTTGATACCGTCAGGGCCATCGCTGCCGCCGGCATTTGCATCCCGCAGCTGGGGCTGGCGCTGGCCGCCCTGATCGGTCGCAAAAACTTCTCGGCAGTGGAGCGCGaagatgacaaa                                                                                     >  2:56278/1‑139 (MQ=255)
                                                                 ccTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCg                                                                                  <  2:555374/149‑1 (MQ=255)
                                                                     cATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTg                                                                            <  2:413628/151‑1 (MQ=255)
                                                                       tGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTgg                                                                           >  2:401340/1‑150 (MQ=255)
                                                                         cTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTg                                                                         >  1:957020/1‑150 (MQ=255)
                                                                          tGATTTGCGTTGCTCAGGGTGTTTATACCGTGGGGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGAt                                                                       >  2:867458/1‑151 (MQ=255)
                                                                           gATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGAt                                                                       <  1:71285/150‑1 (MQ=255)
                                                                           gATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATg                                                                      >  2:1004137/1‑151 (MQ=255)
                                                                            aTTTGCGTTGCTCAGGGAGTTTATACCGTGGTGGCCATCGCTGCCGTGGGAATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATgg                                                                     >  2:64747/1‑151 (MQ=255)
                                                                                 cGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCtgt                                                                >  2:699336/1‑151 (MQ=255)
                                                                                   ttGCTCAGGGTGTTTATACCGTGGTGGCCACCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCt                                                                  <  2:1065593/147‑1 (MQ=255)
                                                                                    tGCTCAGGGTTTTTATCCCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCtgt                                                                <  1:945416/148‑1 (MQ=255)
                                                                                    tGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGtt                                                             >  1:341893/1‑151 (MQ=255)
                                                                                     gCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTGTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTg                                                            >  1:821822/1‑151 (MQ=255)
                                                                                     gCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTg                                                            <  2:276067/151‑1 (MQ=255)
                                                                                     gCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTg                                                            <  1:530155/151‑1 (MQ=255)
                                                                                     gCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTg                                                            >  2:837519/1‑151 (MQ=255)
                                                                                      cTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTgg                                                           >  1:926718/1‑151 (MQ=255)
                                                                                         aGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGt                                                        >  2:20756/1‑151 (MQ=255)
                                                                                         aGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGt                                                        >  2:841650/1‑151 (MQ=255)
                                                                                             tgtTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACt                                                    >  1:160437/1‑151 (MQ=255)
                                                                                               tttATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTAc                                                     >  2:655336/1‑148 (MQ=255)
                                                                                               tttATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTg                                                   >  2:692390/1‑150 (MQ=255)
                                                                                                 tataCCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGa                                                  >  1:250959/1‑149 (MQ=255)
                                                                                                 tataCCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGa                                                  >  2:613569/1‑149 (MQ=255)
                                                                                                     ccGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAggg                                              <  2:1008452/149‑1 (MQ=255)
                                                                                                      cGTGGTGGCTATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCg                                           <  1:379741/151‑1 (MQ=255)
                                                                                                       gtggtgGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGa                                          <  1:228184/151‑1 (MQ=255)
                                                                                                        tggtggCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGa                                          <  2:564729/150‑1 (MQ=255)
                                                                                                        tggtggCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGAt                                         >  2:662710/1‑151 (MQ=255)
                                                                                                          gtggCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTc                                       >  2:560462/1‑151 (MQ=255)
                                                                                                               cATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTc                                  >  1:1036623/1‑151 (MQ=255)
                                                                                                               cATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTc                                  <  1:921551/151‑1 (MQ=255)
                                                                                                                  ctctcccGGGGGTATTTGTATCCCGCCGGTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTGCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTcg                                 <  2:796934/144‑1 (MQ=255)
                                                                                                                  cGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTcgc                                <  1:166119/150‑1 (MQ=255)
                                                                                                                     ttccGTGGGTATTTGTATCCCGCAGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGTCTGTGTTGGCGTTACTGATGGGGCGATTCCTTTCGCCGCt                            >  2:919819/3‑151 (MQ=255)
                                                                                                                     tGCCGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTcgccgc                             >  2:194303/1‑150 (MQ=255)
                                                                                                                       ccGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCt                            <  2:159693/149‑1 (MQ=255)
                                                                                                                       ccGTGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGc                          >  2:677800/1‑151 (MQ=255)
                                                                                                                          tGGGTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCg                        >  2:275687/1‑150 (MQ=255)
                                                                                                                           gggTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGa                       <  1:320926/150‑1 (MQ=255)
                                                                                                                            ggTATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATc                     >  1:864559/1‑151 (MQ=255)
                                                                                                                                    tATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCgt             >  1:276004/1‑151 (MQ=255)
                                                                                                                                       cacgccgcTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTc          >  2:110419/3‑151 (MQ=255)
                                                                                                                                         cgccgcTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCAt        >  2:151250/1‑151 (MQ=255)
                                                                                                                                              ctggggctggggCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCAttcctt   >  1:1077976/1‑151 (MQ=255)
                                                                                                                                                ggggctggggctGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCAttccttc  <  1:820539/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGCGGTGATCATGGGTCTGATGCTGGCGTTCGATATGGGCGGTCCGGTTAACAAAGTGGCCTATGCCTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTC  >  NC_000913/4143069‑4143362

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