Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,699,865 C→G 5.1% P171A (CCG→GCG)  malX → maltose and glucose‑specific PTS enzyme IIB component and IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,699,8650CG 5.1% 295.7 / 4.7 92P171A (CCG→GCG) malXmaltose and glucose‑specific PTS enzyme IIB component and IIC component
Reads supporting (aligned to +/- strand):  ref base C (41/46);  new base G (2/3);  total (43/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

CTGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGTCTGGCCGATTTTCGCCATGGGTATTAGCGGCTTGGGCCATATGATAAACAGCGCGGG  >  NC_000913/1699718‑1700013
                                                                                                                                                   |                                                                                                                                                    
cTGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATACTCGGTGCGGTGATCGCCGGTATTATCGTCTGGAAGCTGCATGAGCGTTACAATAATATCCGCCTGCCg                                                                                                                                                    >  2:990529/1‑150 (MQ=255)
 tGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGTGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCgg                                                                                                                                                   <  2:314736/150‑1 (MQ=255)
  gCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGa                                                                                                                                                  <  2:1445874/150‑1 (MQ=255)
  gCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGa                                                                                                                                                  >  1:1688183/1‑150 (MQ=255)
   ccAACCACGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGAt                                                                                                                                                 >  2:1288397/1‑150 (MQ=255)
       ccACGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATgcgc                                                                                                                                             >  1:1170904/1‑150 (MQ=255)
        cACGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGGCAGGCGGAtattat                                                                                                                                            <  2:711107/150‑6 (MQ=255)
         aCGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGGCAGGCGGAtattatg                                                                                                                                           <  2:1048021/150‑7 (MQ=255)
          cGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTgg                                                                                                                                          <  1:1344788/150‑1 (MQ=255)
          cGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTgg                                                                                                                                          <  1:853405/150‑1 (MQ=255)
          cGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGAGCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTgg                                                                                                                                          <  1:1352159/150‑1 (MQ=255)
           ggATGCCGCGGTTCTGAAAGACAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGc                                                                                                                                         <  2:72969/150‑1 (MQ=255)
            gATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCa                                                                                                                                        >  1:328566/1‑150 (MQ=255)
              tGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCAtt                                                                                                                                      <  2:1668699/150‑1 (MQ=255)
                ccGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCCGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCAttct                                                                                                                                    >  2:550566/1‑150 (MQ=255)
                     gTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGAGCGATACCGGGATCCTCGGGGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTcggc                                                                                                                               <  2:278133/150‑1 (MQ=255)
                       tCTGATAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTcggcgg                                                                                                                             >  1:219159/1‑150 (MQ=255)
                       tCTGAAAGCCATTAACATCCAGCGCATTCTTGGGATCCTGTCGATCTATACCGCGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTcggcgg                                                                                                                             <  1:128012/150‑1 (MQ=255)
                       tCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCTGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTcggcgg                                                                                                                             >  1:213422/1‑150 (MQ=255)
                         tGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGAGTGTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTa                                                                                                                           >  1:913204/1‑150 (MQ=255)
                          gAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGCTGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGtgc                                                                                                                          >  2:275927/1‑148 (MQ=255)
                          gAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTAc                                                                                                                          <  2:915050/150‑1 (MQ=255)
                             aGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCTATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTAcgcg                                                                                                                       <  2:54291/150‑1 (MQ=255)
                                cAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCtt                                                                                                                    >  2:1712310/1‑150 (MQ=255)
                                cAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCtt                                                                                                                    >  1:1414265/1‑150 (MQ=255)
                                cAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCtt                                                                                                                    >  2:61089/1‑150 (MQ=255)
                                 aaTAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGAGCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTc                                                                                                                   <  1:1051929/150‑1 (MQ=255)
                                      cATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTAcc                                                                                                              >  2:1321863/1‑150 (MQ=255)
                                       aTCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCa                                                                                                             >  2:634550/1‑150 (MQ=255)
                                        tCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCCCTGGCATTCATCGGCGGTACGAGCTTCTTACCaa                                                                                                            >  2:89274/1‑150 (MQ=255)
                                         cctgagCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTAGCAAt                                                                                                           <  2:1598333/147‑1 (MQ=255)
                                          cAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAAtt                                                                                                          <  1:1359604/150‑1 (MQ=255)
                                             agCATTCTTGTGATCCAGTCGATCGATACCGGGATGCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCTTTTCCATAATATCCGCCTGGCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCTTACCAATTAta                                                                                                       <  1:596531/150‑2 (MQ=255)
                                               cATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTAtctc                                                                                                     >  1:308041/1‑150 (MQ=255)
                                               cATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTAtctc                                                                                                     <  1:1089074/150‑1 (MQ=255)
                                                 ttCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATctcct                                                                                                   >  1:642078/1‑150 (MQ=255)
                                                   cTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCg                                                                                                 <  1:590983/150‑1 (MQ=255)
                                                      gggATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTg                                                                                              >  2:1068585/1‑150 (MQ=255)
                                                       ggATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTgg                                                                                             <  1:1287628/150‑1 (MQ=255)
                                                       ggATCCAGTCGATCGATACAGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGGTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACAAATTATCTCCTCGCTgg                                                                                             >  1:1654695/1‑150 (MQ=255)
                                                        gATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGt                                                                                            >  1:1822569/1‑150 (MQ=255)
                                                        gATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGt                                                                                            >  1:3517/1‑150 (MQ=255)
                                                           ccAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGAt                                                                                         >  2:181858/1‑150 (MQ=255)
                                                            cAGTCGATCGATACCGGGATACTCGGTGCGGTGATAGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCAGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATg                                                                                        >  1:836373/1‑150 (MQ=255)
                                                               tcgatcgaTACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGc                                                                                     <  1:440712/150‑1 (MQ=255)
                                                                cgatcgatACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGcc                                                                                    <  1:1576314/150‑1 (MQ=255)
                                                                cgatcgatACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGcc                                                                                    >  1:1754519/1‑150 (MQ=255)
                                                                 gatcgatACCGGGATCCTCGTTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCt                                                                                   <  1:292823/150‑1 (MQ=255)
                                                                     gatACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGTTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTc                                                                               >  2:639559/1‑150 (MQ=255)
                                                                         ccGGGATCCGCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGcc                                                                           <  2:80636/150‑1 (MQ=255)
                                                                          cGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCt                                                                          <  1:1809746/150‑1 (MQ=255)
                                                                          cGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCt                                                                          <  2:1170904/150‑1 (MQ=255)
                                                                          cGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCt                                                                          <  1:1563841/150‑1 (MQ=255)
                                                                           gggATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTg                                                                         <  2:1513068/150‑1 (MQ=255)
                                                                           gggATCCTCGGTCCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCATCGCTGGTGATGGGCCTTGTCGGCCTg                                                                         >  1:472502/1‑150 (MQ=255)
                                                                               tCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGa                                                                     >  2:10285/1‑150 (MQ=255)
                                                                                ccTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGAt                                                                    >  2:1592286/1‑150 (MQ=255)
                                                                                 cTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGAtt                                                                   <  1:1010216/150‑1 (MQ=255)
                                                                                 cTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGAtt                                                                   <  2:1688183/150‑1 (MQ=255)
                                                                                      tgcggtgATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCAtt                                                              <  1:1608162/150‑1 (MQ=255)
                                                                                       gcggtgATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCAttt                                                             >  2:1719572/1‑149 (MQ=255)
                                                                                       gcggtgATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATACCATTa                                                             <  1:897569/150‑1 (MQ=255)
                                                                                        cggtgATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCCTGTCGGCCTGGTGATTCCATTAg                                                            <  2:1414265/150‑1 (MQ=255)
                                                                                         cgtgATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGt                                                           >  1:1660145/2‑150 (MQ=255)
                                                                                          gtgATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGTc                                                          >  1:131125/1‑150 (MQ=255)
                                                                                            gATCGACGGTATTATCGTCTGCATGCTGCAAGAGCGTTTCCATAATATCCGCCTGCCGGATGAGCTGGCATTCTTCGGCGCTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTTATTCCATTAGTCtt                                                        >  2:1541917/1‑149 (MQ=255)
                                                                                             aTCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGTCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGTCTgg                                                       >  1:1501115/1‑150 (MQ=255)
                                                                                                 ccGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGTCTGGCCGa                                                   >  2:11486/1‑150 (MQ=255)
                                                                                                          tCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGTCTGGCCGATTTTCTCCa                                          >  1:221442/1‑150 (MQ=255)
                                                                                                              cTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGTGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGTCTGGCCGATTTTCGCCATggg                                      <  2:449683/150‑1 (MQ=255)
                                                                                                              cTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGTCTGGCCGATTTTCGCCATggg                                      <  2:1270727/150‑1 (MQ=255)
                                                                                                              cTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCGGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGACTGGCCGATTTTCGCCATggg                                      <  1:1783148/150‑1 (MQ=255)
                                                                                                                ggATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGTCTGGCCGATTTTCGCCATGGGTa                                    <  2:1271900/150‑1 (MQ=255)
                                                                                                                ggATGCTGCATGAGCGTTTCCATAATATCCGCCTGCAGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGTCTGGCCGATTTTCGCCATGGGTa                                    >  1:1326243/1‑150 (MQ=255)
                                                                                                                ggATGCTGCATGAGCG‑TTCCATAATATCCGCCTGCCGGATGCGGTGGCATTATTAGGCGGTACGCGCTTCGTACCAATTAACTCCTACCTGGTGATGGGCCTTGTCGGTCTGGTGATTCCATTAGTCAGGCCGATTTTCGCCATGGGTAt                                   >  1:1230670/1‑150 (MQ=255)
                                                                                                                 gATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGTCTGGCCGATTTTCGCCATggggta                                   <  1:1397056/150‑4 (MQ=255)
                                                                                                                  aTGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGTCTGGCCGATTTTCGCCATGGGTAtt                                  <  1:61089/150‑1 (MQ=255)
                                                                                                                     ctgcATGAGCGTTTAGATAATATCCGCCTGCCGGATGCGCGTGCATTCTTCGGCGGTACGCGCTTCGTACCATTTATCTCTTCGCTGGTGATGGGCCTTGGTGGCCTTGGTATTCCATTTGTGTGTCAGATTTTCTCCATGGGTATTatc                               <  2:1133194/150‑3 (MQ=255)
                                                                                                                      tgcATGAGCGTTTCCATAATATCCGCCTGCGGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGTCTGGCCGATTTTCGCCATGGGTATTAGCg                              <  1:665274/150‑1 (MQ=255)
                                                                                                                      tgcATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCATTTTCTGCCTGTTGATTCCATTAGTCTGTCCGATTTTCGACATGGGTATTAGCg                              >  2:299654/1‑150 (MQ=255)
                                                                                                                           gAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGTCTGGCCGATTTTCGCCATGGGTATTAGCGGCTTg                         <  2:1268219/150‑1 (MQ=255)
                                                                                                                             gCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCTGAGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTAGTCGGCCTGGTGATTCCATCAGTCTGGCCGATTTTCGCCATGGGTATTAGCGGCTTggg                       <  2:1339452/150‑1 (MQ=255)
                                                                                                                                tttCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGTCTGGCCGATTTTCGCCATGGGTATTAGCGGCTTGGGCCa                    <  2:1385327/150‑1 (MQ=255)
                                                                                                                                   ccATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGTCTGGCCGATTTTCGCCATGGGTATTAGCGGCTTGGGCCatat                 <  2:1724466/150‑1 (MQ=255)
                                                                                                                                    cATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGGCGGCCTGGTGATTCCATTAGTCTGGCCGATTTTCGCCATGGGTATTAGCGGCTTGGGCCATATg                <  2:73876/150‑1 (MQ=255)
                                                                                                                                     ataataTCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGATGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGTCTGGCCTATTTTCGCCATGGGTATTAGCGGCTTGGGCCATATGa               >  1:299680/1‑150 (MQ=255)
                                                                                                                                     agcgcatCCGGCAGGCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGTCTGGCCGATTTTCGCCATGGGTATTAGCGGCTTGGGCCATATGa               >  1:68129/6‑150 (MQ=255)
                                                                                                                                     agcgcatCCGGCAGGCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGTCTGGCCGATTTTCGCCATGGGTATTAGCGGCTTGGGCCATATGa               >  1:1566962/6‑150 (MQ=255)
                                                                                                                                        atatCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGTCTGGCCGATTTTCGCCATGGGTATTAGCGGCTTGGGCCATATGATaa            <  1:573185/150‑1 (MQ=255)
                                                                                                                                          atCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGGCGGCCTGGTGATTCCATTAGTCTGGCCGATTTTCGCCATGGGTATTAGCGGCTTGGGCCATATGATAAAc          <  2:1449098/150‑1 (MQ=255)
                                                                                                                                           tCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGTCTGGCCGATTTTCGCCATGGGTATTAGCGGCTTGGGCCATATGATAAACa         >  2:558916/1‑150 (MQ=255)
                                                                                                                                                  gCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGTCTGGCCGATTTTCGCCATGGGTATTAGCGGCTTGGGCCATATGATAAACAGCGCggg  <  1:550566/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATCCTCGGTGCGGTGATCGCCGGTATTATCGTCTGGATGCTGCATGAGCGTTTCCATAATATCCGCCTGCCGGATGCGCTGGCATTCTTCGGCGGTACGCGCTTCGTACCAATTATCTCCTCGCTGGTGATGGGCCTTGTCGGCCTGGTGATTCCATTAGTCTGGCCGATTTTCGCCATGGGTATTAGCGGCTTGGGCCATATGATAAACAGCGCGGG  >  NC_000913/1699718‑1700013

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