Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,758,359 C→G 6.3% G123R (GGT→CGT) ‡ ldtE ← L,D‑transpeptidase LdtE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,758,3590CG 6.3% 299.5 / 7.4 90G123R (GGT→CGT) ‡ldtEL,D‑transpeptidase LdtE
Reads supporting (aligned to +/- strand):  ref base C (42/42);  new base G (2/4);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TGGTCCGGCAGGAACGACTGGCGGTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTGGC  >  NC_000913/1758210‑1758501
                                                                                                                                                     |                                                                                                                                              
tGGTCCGGCAGGAACGACTGGCGGTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCATGCAATCCAATAcc                                                                                                                                                <  2:988090/150‑1 (MQ=255)
   tCCGGCAGGAACGACTGGCGGTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTAt                                                                                                                                             >  1:1958084/1‑150 (MQ=255)
    ccGGCAGGAACGACTGGCGGTAATTTAATGCAACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTAtt                                                                                                                                            >  2:1805985/1‑150 (MQ=255)
     cGGCAGGAACGACTGGCGGTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTg                                                                                                                                           <  1:704601/150‑1 (MQ=255)
          ggAACGACTGGCGGTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATAGGTATTGGAttg                                                                                                                                      <  1:1095638/150‑3 (MQ=255)
           gAACGACTGGCGGTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATAGGTATTGGAttgc                                                                                                                                     <  2:941907/150‑4 (MQ=255)
             aCGACTGGCGGTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATAGGTATTGGAttgcag                                                                                                                                   <  1:1655849/150‑6 (MQ=255)
             aCGACTGGCGGTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAAcc                                                                                                                                   >  2:1545901/1‑150 (MQ=255)
             aCGACTGGCGGTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAAcc                                                                                                                                   <  1:346328/150‑1 (MQ=255)
             aCGACTGGCGGTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTGGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAAcc                                                                                                                                   <  2:297633/150‑1 (MQ=255)
                  tGGCGGTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGcac                                                                                                                              <  1:1877209/150‑1 (MQ=255)
                   ggCGGTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGcaca                                                                                                                             >  1:1653466/1‑150 (MQ=255)
                    gCGGTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTATGCCCACCCCGCGTTTCCATCACCGGCTTTTCCAGCCAATGCAATCCCATACCTATTGGATAAACATGCACaa                                                                                                                            >  1:984408/1‑150 (MQ=255)
                      ggTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAata                                                                                                                          >  1:1040139/1‑150 (MQ=255)
                      ggTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAata                                                                                                                          >  1:448190/1‑150 (MQ=255)
                          aTTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATAtttt                                                                                                                      >  1:909406/1‑150 (MQ=255)
                           tttAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTtc                                                                                                                     <  1:132916/150‑1 (MQ=255)
                           tttAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTtc                                                                                                                     <  2:1762760/150‑1 (MQ=255)
                            ttAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTtct                                                                                                                    >  2:1086029/1‑150 (MQ=255)
                             tAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTtctc                                                                                                                   >  1:276178/1‑150 (MQ=255)
                               aTGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTccc                                                                                                                 >  1:1787494/1‑150 (MQ=255)
                                  ccACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAAcc                                                                                                                                   <  2:1250752/129‑1 (MQ=255)
                                  ccACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAAcc                                                                                                                                   >  1:1250752/1‑129 (MQ=255)
                                  ccACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCcggc                                                                                                              >  1:182245/1‑150 (MQ=255)
                                          cAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGataata                                                                                                      >  1:1277471/1‑150 (MQ=255)
                                               aaCGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATaaa                                                                                                 <  1:1166535/150‑1 (MQ=255)
                                                aCGCTGACGAATGCCTGCCGTAGGCGTCGAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAg                                                                                                <  2:100389/150‑1 (MQ=255)
                                                  gCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGTACAATATTTTCTCCCGGCGGATAATAATAAAGgc                                                                                              <  1:1354245/150‑1 (MQ=255)
                                                    tGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGgcgc                                                                                            <  2:909406/150‑1 (MQ=255)
                                                     gACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGAGGATAATAATAAAGGCGCa                                                                                           >  2:625196/1‑150 (MQ=255)
                                                             cctgccGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCa                                                                                   >  1:1176862/1‑150 (MQ=255)
                                                               tgccGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAg                                                                                 >  1:1398804/1‑150 (MQ=255)
                                                                gccGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTtct                                                                                                                    >  2:1063666/1‑114 (MQ=255)
                                                                gccGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTtct                                                                                                                    <  1:1063666/114‑1 (MQ=255)
                                                                    tAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATaaa                                                                                                                                     <  1:538886/93‑1 (MQ=255)
                                                                    tAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATaaa                                                                                                                                     >  2:538886/1‑93 (MQ=255)
                                                                      ggCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACg                                                                          >  1:332846/1‑150 (MQ=255)
                                                                       gCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGa                                                                         <  2:276178/150‑1 (MQ=255)
                                                                        cGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGAt                                                                        >  1:831347/1‑150 (MQ=255)
                                                                        cGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGAt                                                                        >  1:1060852/1‑150 (MQ=255)
                                                                        cGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGAt                                                                        <  2:259863/150‑1 (MQ=255)
                                                                         gTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATATTAATAAAGGCGCAGCTCGGCAAGGTTAACGATa                                                                       >  2:1392213/1‑150 (MQ=255)
                                                                              ggTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATccc                                                                  <  2:1958084/150‑1 (MQ=255)
                                                                                 ggggTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTGCCAGCCCCTGCAATCCACTACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTg                                                               <  2:1699083/149‑1 (MQ=255)
                                                                                   ggTTAGTGATTTTCTGCCCTACACGCGGTTCCATCACCGGCGTTTCCAGCCCCTGCCATCCTATACGTATTGGATATAGCTGCACATTTTTTTCTCCCGGCGGCTAATAATTAAGGCGCAGCTCTGCAAGGTTAACGATCATCCCCTGAc                                                             <  2:857151/150‑1 (MQ=255)
                                                                                   ggTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGAc                                                             >  1:323186/1‑150 (MQ=255)
                                                                                    gTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGAcg                                                            >  1:1062873/1‑150 (MQ=255)
                                                                                       aGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCgg                                                         >  1:338422/1‑150 (MQ=255)
                                                                                            ttttCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCAt                                                    >  2:1970032/1‑150 (MQ=255)
                                                                                             tttCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCACCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATc                                                   <  1:1230156/150‑1 (MQ=255)
                                                                                              ttCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCa                                                  <  2:1210564/150‑1 (MQ=255)
                                                                                              ttCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCGCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCa                                                  <  2:1545603/150‑1 (MQ=255)
                                                                                                cTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAgg                                                >  1:125429/1‑150 (MQ=255)
                                                                                                 tGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGt                                               <  2:1653466/150‑1 (MQ=255)
                                                                                                  gCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGta                                              >  1:1102391/1‑150 (MQ=255)
                                                                                                  gCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGta                                              >  2:218908/1‑150 (MQ=255)
                                                                                                     cAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGtaata                                           <  1:1154085/150‑1 (MQ=255)
                                                                                                      aaCACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGtaataa                                          >  1:288042/1‑150 (MQ=255)
                                                                                                       acacGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAAc                                         <  2:338422/150‑1 (MQ=255)
                                                                                                       acacGCGGTTCCATCACCGGCGTTTCCAGCAACTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCAGGCGGATAATAATAAAGGCGCATAGCTGCAAGGTTAACGATAATCCCCTGACGCGGGGCATCAGGTAATAAc                                         <  2:984408/150‑1 (MQ=255)
                                                                                                           gcgTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATGAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGtt                                     <  2:1063957/150‑1 (MQ=255)
                                                                                                           gcgTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATGAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGtt                                     <  2:1056988/150‑1 (MQ=255)
                                                                                                              tttCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTtgtg                                  <  2:831347/150‑1 (MQ=255)
                                                                                                              tttCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTtgtg                                  <  2:448190/150‑1 (MQ=255)
                                                                                                               ttCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACTCGGTGCATCAGGTAATAACAGTTGTGa                                 <  2:1398804/150‑1 (MQ=255)
                                                                                                               ttCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGa                                 <  1:1282278/150‑1 (MQ=255)
                                                                                                                 ccATCACCGGCGTTTCCAGTCCCTGCAATCCAATACCTATTGGATAAACCTGCACAAAATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCATGACGAGGTGCATCAGCTAATAACAGTTGTGAAg                               >  2:937235/1‑150 (MQ=255)
                                                                                                                 ccATCACCGGCGTTTCCAGCCCATGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATTATAAAGGCGCAGGTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGGATCATGTAATAACAGTTGTGAAg                               <  2:1043698/150‑1 (MQ=255)
                                                                                                                  cATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAgg                              <  2:332846/150‑1 (MQ=255)
                                                                                                                  cATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAAGATTTTCTCCCGGCGGATTATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAGTCCCCTGACGCGGTGCATCAGGTAATAACAGTTGGGAAgg                              >  2:1799409/1‑150 (MQ=255)
                                                                                                                     cACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAAt                           >  1:740729/1‑150 (MQ=255)
                                                                                                                      aCCGGCTTTTCCAGTCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATa                          <  1:937235/150‑1 (MQ=255)
                                                                                                                      aCCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTTATAACAGTTGTGAAGGAATa                          >  2:1384801/1‑150 (MQ=255)
                                                                                                                            gTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTAtc                    >  1:1679293/1‑150 (MQ=255)
                                                                                                                               tCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTAtcgtc                 >  1:73872/1‑150 (MQ=255)
                                                                                                                                ccAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTAtcgtcg                <  1:625196/150‑1 (MQ=255)
                                                                                                                                 cAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATcgtcgt               <  1:1794704/150‑1 (MQ=255)
                                                                                                                                     cccTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCa           >  1:255800/1‑150 (MQ=255)
                                                                                                                                      aggtttaTCCAATAGGTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAg          >  2:1852553/7‑150 (MQ=255)
                                                                                                                                       ggtttaTCCAATAGGTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAgg         >  1:1467590/6‑150 (MQ=255)
                                                                                                                                       cTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAgg         >  1:1327936/1‑150 (MQ=255)
                                                                                                                                       cTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAgg         <  2:182245/150‑1 (MQ=255)
                                                                                                                                        tGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGt        >  2:1444723/1‑150 (MQ=255)
                                                                                                                                          cAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGttt      >  1:1852757/1‑150 (MQ=255)
                                                                                                                                           aaTCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGtttt     <  2:1062873/150‑1 (MQ=255)
                                                                                                                                            aTCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTg    <  2:288042/150‑1 (MQ=255)
                                                                                                                                             tCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTgg   <  1:218908/150‑1 (MQ=255)
                                                                                                                                              ccAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTg    >  2:1369548/1‑148 (MQ=255)
                                                                                                                                              ccAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTg    <  1:1369548/148‑1 (MQ=255)
                                                                                                                                              ccAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTGGc  <  1:1805985/150‑1 (MQ=255)
                                                                                                                                              ccAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTGGc  <  2:323186/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
TGGTCCGGCAGGAACGACTGGCGGTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTGGC  >  NC_000913/1758210‑1758501

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