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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,758,3590CG 6.7% 318.3 / 11.8 90G123R (GGT→CGT) ‡ldtEmurein L,D‑transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (44/40);  new base G (3/3);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.88e-01

AGGAACGACTGGCGGTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTGGCAC  >  NC_000913/1758219‑1758503
                                                                                                                                            |                                                                                                                                                
aGGAACGACTGGCGGTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGAt                                                                                                                                          >  4:715464/1‑149 (MQ=255)
         tGGCGGTAATTTAATGCCACGCAAAAGCGAACGCTGACGAATGCCTGCCGTAGGCGTTCAGGATGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGca                                                                                                                                 >  2:26164/1‑149 (MQ=255)
           gCGGTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGcaca                                                                                                                               >  3:968192/1‑149 (MQ=255)
              gTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTgg                                                                                                                                            <  3:276301/133‑1 (MQ=255)
              gTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTgg                                                                                                                                            >  4:276301/1‑133 (MQ=255)
                aaTTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATAtt                                                                                                                          <  4:1456237/149‑1 (MQ=255)
                   ttAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTtc                                                                                                                       <  4:968192/149‑1 (MQ=255)
                   ttAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTtc                                                                                                                       >  3:375460/1‑149 (MQ=255)
                       tGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTccc                                                                                                                   >  1:1247090/1‑149 (MQ=255)
                          cACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCcggc                                                                                                                >  4:681084/1‑149 (MQ=255)
                           aCGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATAGGTATTGGAtt                                                                                                                                         <  3:977796/123‑2 (MQ=255)
                           aCGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATAGGTATTGGAtt                                                                                                                                         >  4:977796/1‑122 (MQ=255)
                           aCGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCcggcg                                                                                                               >  3:1392820/1‑149 (MQ=255)
                             gCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTTCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCGCCCGGCGGa                                                                                                             <  4:1392820/149‑1 (MQ=255)
                              cTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATAtttt                                                                                                                        >  1:794515/1‑137 (MQ=255)
                              cTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATAtttt                                                                                                                        <  2:794515/137‑1 (MQ=255)
                              cTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGat                                                                                                            <  3:715464/149‑1 (MQ=255)
                               tCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTtc                                                                                                                       >  4:1417297/1‑137 (MQ=255)
                               tCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTtc                                                                                                                       <  3:1417297/137‑1 (MQ=255)
                                ccAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGataa                                                                                                          <  4:892670/149‑1 (MQ=255)
                                  aGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGataata                                                                                                        <  3:1301910/149‑1 (MQ=255)
                                   gCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACGATATTTTCTCCCGGCGGataataa                                                                                                       <  4:1251528/149‑1 (MQ=255)
                                     gAACGCTGACGAATGCCTGCCGTAGGCGTCTAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAata                                                                                                                            >  4:982356/1‑126 (MQ=255)
                                     gAACGCTGACGAATGCCTGCCGTAGGCGTCTAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAata                                                                                                                            <  3:982356/126‑1 (MQ=255)
                                     gAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGataataata                                                                                                     >  4:176818/1‑149 (MQ=255)
                                         gCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAgg                                                                                                 <  2:1247090/149‑1 (MQ=255)
                                         gCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAgg                                                                                                 >  1:753774/1‑149 (MQ=255)
                                           tGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGgcg                                                                                               >  3:664909/1‑149 (MQ=255)
                                           tGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGgcg                                                                                               >  3:1167102/1‑149 (MQ=255)
                                            gACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATAGGTATTGGAtt                                                                                                                                         >  4:509829/1‑105 (MQ=255)
                                            gACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATAGGTATTGGAtt                                                                                                                                         <  3:509829/106‑2 (MQ=255)
                                               gAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATAtttt                                                                                                                        >  1:490942/1‑120 (MQ=255)
                                               gAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATAtttt                                                                                                                        <  2:490942/120‑1 (MQ=255)
                                                     ctgccGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGataat                                                                                                         >  3:596521/1‑129 (MQ=255)
                                                     ctgccGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGataat                                                                                                         <  4:596521/129‑1 (MQ=255)
                                                       gccGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAg                                                                                   <  4:375460/149‑1 (MQ=255)
                                                        ccGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAgg                                                                                  <  2:856702/149‑1 (MQ=255)
                                                        ccGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAgg                                                                                  <  2:753774/149‑1 (MQ=255)
                                                                gTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGAt                                                                          >  2:665319/1‑149 (MQ=255)
                                                                 tCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGataataa                                                                                                       <  2:154423/119‑1 (MQ=255)
                                                                 tCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGataataa                                                                                                       >  1:154423/1‑119 (MQ=255)
                                                                 tCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATa                                                                         >  1:793516/1‑149 (MQ=255)
                                                                    aGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATc                                                                      <  4:1167102/149‑1 (MQ=255)
                                                                      gTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCt                                                                                                                                    <  4:392200/85‑1 (MQ=255)
                                                                      gTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCt                                                                                                                                    >  3:392200/1‑85 (MQ=255)
                                                                      gTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATccc                                                                    <  1:665319/149‑1 (MQ=255)
                                                                           gTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCcggcg                                                                                                               <  1:441542/101‑1 (MQ=255)
                                                                           gTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCcggcg                                                                                                               >  2:441542/1‑101 (MQ=255)
                                                                            ttAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGAcg                                                              >  3:1121788/1‑149 (MQ=255)
                                                                                  aTTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGc                                                        >  4:531832/1‑149 (MQ=255)
                                                                                  aTTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGc                                                        >  1:61969/1‑149 (MQ=255)
                                                                                    tttCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCAt                                                      <  1:26164/149‑1 (MQ=255)
                                                                                     ttCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATc                                                     <  4:664909/149‑1 (MQ=255)
                                                                                        tGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAgg                                                  <  3:531832/149‑1 (MQ=255)
                                                                                            cAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATCATAATAAAGgc                                                                                                <  3:1157366/99‑1 (MQ=255)
                                                                                            cAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGgc                                                                                                >  4:1157366/1‑99 (MQ=255)
                                                                                            cAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGtaat                                              <  3:176818/149‑1 (MQ=255)
                                                                                              acacGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGtaataa                                            <  3:681084/149‑1 (MQ=255)
                                                                                                acGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACa                                          >  2:1088910/1‑149 (MQ=255)
                                                                                                acGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACa                                          >  4:569177/1‑149 (MQ=255)
                                                                                                  gcgTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGctct                                                                                        <  2:150310/101‑1 (MQ=255)
                                                                                                  gcgTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGctct                                                                                        >  1:150310/1‑101 (MQ=255)
                                                                                                     tttCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTtgt                                     >  2:870223/1‑149 (MQ=255)
                                                                                                        ccATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGaa                                  <  1:870223/149‑1 (MQ=255)
                                                                                                          aTCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAgg                                >  3:809578/1‑149 (MQ=255)
                                                                                                          aTCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAgg                                >  1:1326014/1‑149 (MQ=255)
                                                                                                               cGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAAc                                                                             >  2:597179/1‑99 (MQ=255)
                                                                                                               cGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAAc                                                                             <  1:597179/99‑1 (MQ=255)
                                                                                                                ggCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGt                          >  3:1352447/1‑149 (MQ=255)
                                                                                                                   gTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTAt                       >  3:1019143/1‑149 (MQ=255)
                                                                                                                    tttCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTAtc                      <  3:569177/149‑1 (MQ=255)
                                                                                                                       ccAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGt                                                          >  3:895430/1‑110 (MQ=255)
                                                                                                                       ccAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGt                                                          <  4:895430/110‑1 (MQ=255)
                                                                                                                          gCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCt                                                                  >  2:294746/1‑99 (MQ=255)
                                                                                                                          gCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCt                                                                  <  1:294746/99‑1 (MQ=255)
                                                                                                                          gCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATa                            >  2:583863/1‑137 (MQ=255)
                                                                                                                          gCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATa                            <  1:583863/137‑1 (MQ=255)
                                                                                                                            cccTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTAcc              >  3:743755/1‑149 (MQ=255)
                                                                                                                              cTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAg            >  1:802556/1‑149 (MQ=255)
                                                                                                                                  taTCCAATAGGTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAg                                                   >  1:170728/2‑106 (MQ=255)
                                                                                                                                  taTCCAATAGGTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAg                                                   <  2:170728/105‑1 (MQ=255)
                                                                                                                                  aaTCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGttt        <  2:61969/149‑1 (MQ=255)
                                                                                                                                   aTCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGtttt       >  1:122928/1‑149 (MQ=255)
                                                                                                                                   aTCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACAAGGtttt       <  2:802556/149‑1 (MQ=255)
                                                                                                                                    tCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTg      <  2:1326014/149‑1 (MQ=255)
                                                                                                                                    tCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGATAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTg      >  4:1269231/1‑149 (MQ=255)
                                                                                                                                      cAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGt                                                 <  4:439413/104‑1 (MQ=255)
                                                                                                                                      cAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGt                                                 >  3:439413/1‑104 (MQ=255)
                                                                                                                                       aaTACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTGGCa   >  3:708305/1‑149 (MQ=255)
                                                                                                                                        aTACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTGGCAc  <  2:793516/149‑1 (MQ=255)
                                                                                                                                            cTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTgg     >  3:989943/1‑142 (MQ=255)
                                                                                                                                            cTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTgg     <  4:989943/142‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                
AGGAACGACTGGCGGTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTGGCAC  >  NC_000913/1758219‑1758503

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