Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,445,973 G→C 6.9% P138A (CCG→GCG)  mepA ← peptidoglycan DD‑endopeptidase/peptidoglycan LD‑endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,445,9730GC 6.9% 241.1 / 7.0 72P138A (CCG→GCG) mepApeptidoglycan DD‑endopeptidase/peptidoglycan LD‑endopeptidase
Reads supporting (aligned to +/- strand):  ref base G (32/35);  new base C (3/2);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

GCAAAGTTGTTGTTTAATCGCCGGATTAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCGGGCATCCCCATATCGCCAATCAGCACCGTACCCATGCCCA  >  NC_000913/2445824‑2446122
                                                                                                                                                     |                                                                                                                                                     
gCAAAGTTGTTGTTTAATCGCCGGATTAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCgg                                                                                                                                                       <  2:1112830/150‑1 (MQ=255)
 cAAAGTTGTTGTTTAATCGCCGGATTAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCggg                                                                                                                                                      <  1:630656/150‑1 (MQ=255)
 cAAAGTTGTTGTTTAATCGCCGGATTAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCggg                                                                                                                                                      <  1:949722/150‑1 (MQ=255)
    aGTTGTTGTTTAATCGCCGGATTAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGgcgc                                                                                                                                                   <  2:516995/150‑1 (MQ=255)
     gttgttgttTAATCGCCGGATTAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCa                                                                                                                                                  >  1:1391350/1‑150 (MQ=255)
       tgttgttTAATCGCCGGATTAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAag                                                                                                                                                <  1:794280/150‑1 (MQ=255)
         ttgttTAATCGCCGGATTAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAagag                                                                                                                                              <  2:520559/150‑1 (MQ=255)
          tgttTAATCGCCGGATTAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTTCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGc                                                                                                                                             <  1:789049/150‑1 (MQ=255)
          tgttTAATCGCCGGATTAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGATACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGc                                                                                                                                             <  2:541113/150‑1 (MQ=255)
          tgttTAATCGCCGGATTAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGc                                                                                                                                             <  2:1541219/150‑1 (MQ=255)
          tgttTAATCGCCGGATTAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGc                                                                                                                                             <  1:630513/150‑1 (MQ=255)
          tgttTAATCGCCGGATTAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGCTCTTGCGCCCGCAagcac                                                                                                                                             >  2:1419894/1‑147 (MQ=255)
           ggttAATCGCCGGATTAACAAAAATGCGCGGGACGTCTTTGTCCTGGGCGGCGAGTTTGTTCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTGTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGCCCGCAagcact                                                                                                                                            <  2:391159/148‑5 (MQ=255)
             ttAATCGCCGGATTAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGCTCTTGCGCCCGCAagcactgg                                                                                                                                          >  2:638842/1‑144 (MQ=255)
             ttAATCGCCGGATTAACAAAAATGAGCGTGACGTCTTTGTCCTTGGCGTCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTAGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTgc                                                                                                                                          >  2:1616461/1‑150 (MQ=255)
              tAATCGCCGGATTAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGCTCTTGCGCCCGCAagcactgga                                                                                                                                         <  1:359880/150‑8 (MQ=255)
                aTCGCCGGATTAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGcgcg                                                                                                                                       >  1:1553952/1‑150 (MQ=255)
                  cGCCGGATTAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGa                                                                                                                                     >  1:277690/1‑150 (MQ=255)
                   gCCGGATTAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAg                                                                                                                                    >  1:231433/1‑150 (MQ=255)
                    ccGGATTAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAgg                                                                                                                                   <  2:1433727/150‑1 (MQ=255)
                       gATTAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTcc                                                                                                                                >  2:655772/1‑150 (MQ=255)
                        aTTAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCa                                                                                                                               <  1:818550/150‑1 (MQ=255)
                          tAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGc                                                                                                                             >  1:1208012/1‑150 (MQ=255)
                           aaCAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCg                                                                                                                            <  1:899665/150‑1 (MQ=255)
                           aaCAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCg                                                                                                                            >  1:527808/1‑150 (MQ=255)
                           aaCAAAAATGCGCGTGACGGCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCg                                                                                                                            <  2:1411509/150‑1 (MQ=255)
                            aCAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGa                                                                                                                           >  2:194044/1‑150 (MQ=255)
                             cAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAg                                                                                                                          >  1:740920/1‑150 (MQ=255)
                                   tGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCg                                                                                                                    <  1:655772/150‑1 (MQ=255)
                                     cgcgTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGc                                                                                                                  <  2:229280/150‑1 (MQ=255)
                                         tGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGtt                                                                                                              >  1:1123012/1‑150 (MQ=255)
                                              tCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGa                                                                                                         <  1:955720/150‑1 (MQ=255)
                                               cTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGaa                                                                                                        >  2:39554/1‑150 (MQ=255)
                                                  tGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGa                                                                                                     <  2:1581964/150‑1 (MQ=255)
                                                   gTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGat                                                                                                    >  1:1381950/1‑150 (MQ=255)
                                                             gcgAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATcc                                                                                          <  1:760205/150‑1 (MQ=255)
                                                              cgAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCa                                                                                                                                                  <  1:1450068/93‑1 (MQ=255)
                                                              cgAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCa                                                                                                                                                  >  2:1450068/1‑93 (MQ=255)
                                                                 gTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGc                                                                                                                  >  1:1627730/1‑122 (MQ=255)
                                                                 gTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGc                                                                                                                  <  2:1627730/122‑1 (MQ=255)
                                                                  tttGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTcc                                                                                     <  2:1582836/150‑1 (MQ=255)
                                                                         aaGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCtgg                                                                              <  2:1208012/150‑1 (MQ=255)
                                                                            cTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCtggtgg                                                                           <  1:655439/150‑1 (MQ=255)
                                                                                      cTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCAATTGGTTGCGGGCGAAAGAGATGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAATATATCGACATCCCGTCCGGTATGTTGGCTGGCATGAc                                                                 >  2:147496/1‑150 (MQ=255)
                                                                                        ggCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGAccg                                                               <  2:277690/150‑1 (MQ=255)
                                                                                        ggCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGAccg                                                               >  1:726447/1‑150 (MQ=255)
                                                                                          cccccACAGCCTCGAGACAACCTGTTTACCGTCGCGGCAAACTAATTCCAGTGCTTGCGGGCTCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGAccgcc                                                             <  2:726447/147‑1 (MQ=255)
                                                                                          cTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGAccgcc                                                             >  2:1551506/1‑150 (MQ=255)
                                                                                            tCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGt                                                           <  1:407352/150‑1 (MQ=255)
                                                                                               acaGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGGCCGGTCTGGTGGCTGGCATGACCGCCGTTGa                                                        <  2:731511/150‑1 (MQ=255)
                                                                                                  gCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAAc                                                     <  2:1553952/150‑1 (MQ=255)
                                                                                                  gCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAAc                                                     <  1:763385/150‑1 (MQ=255)
                                                                                                     tGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGcc                                                  >  1:1166928/1‑150 (MQ=255)
                                                                                                               gtgtTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCggg                                        <  2:1391350/150‑1 (MQ=255)
                                                                                                               gtgtTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCggg                                        <  2:474558/150‑1 (MQ=255)
                                                                                                                   ttACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCGGGCATc                                    >  1:296130/1‑150 (MQ=255)
                                                                                                                    tACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCGGGCATcc                                   >  2:1307437/1‑150 (MQ=255)
                                                                                                                    tACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCGGGCATcc                                   >  1:205182/1‑150 (MQ=255)
                                                                                                                      ccGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCGGGCATcccc                                 <  1:1307437/150‑1 (MQ=255)
                                                                                                                            cgGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCGGGCATCCCCATATCg                           >  2:1587732/1‑150 (MQ=255)
                                                                                                                            cgGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCGGGCATCCCCATATCg                           >  1:1190622/1‑150 (MQ=255)
                                                                                                                            cgGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCGGGCATCCCCATATCg                           >  1:156433/1‑150 (MQ=255)
                                                                                                                             gggAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCa                                                                     >  1:706106/1‑107 (MQ=255)
                                                                                                                             gggAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCa                                                                     <  2:706106/107‑1 (MQ=255)
                                                                                                                                 aaCTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCGGGCATCCCCATATCGCCAAt                      <  2:740920/150‑1 (MQ=255)
                                                                                                                                  aCTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCGGGCATCCCCATATCGCCAATc                     <  1:39554/150‑1 (MQ=255)
                                                                                                                                     aaGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCGGGCATCCCCATATCGCCAATCAGc                  >  2:13856/1‑150 (MQ=255)
                                                                                                                                       gTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCGGGCATCCCCATATCGCCAATCAGCAc                <  1:539288/150‑1 (MQ=255)
                                                                                                                                       gTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCGGGCATCCCCATATCGCCAATCAGCAc                >  1:1390863/1‑150 (MQ=255)
                                                                                                                                        tCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCGGGCATCCCCATATCGCCAATCAGCAcc               >  2:239959/1‑150 (MQ=255)
                                                                                                                                           agctcTTGCGCCCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCGGGCATCCCCATATCGCCAATCAGCACCGTa            >  2:894486/5‑150 (MQ=255)
                                                                                                                                               cTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCGGGCATCCCCATATCGCCAATCAGCACCGTACCCa        >  1:786486/1‑150 (MQ=255)
                                                                                                                                                 tGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCGGGCATCCCCATATCGCCAATCAGCACCGTACCCATg      <  2:231433/150‑1 (MQ=255)
                                                                                                                                                 tGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCGGGCATCCCCATATCGCCAATCAGCACCGTACCCATg      >  2:947598/1‑150 (MQ=255)
                                                                                                                                                  gCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGCCATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCCGCGGGCATCCCCATATCGCCAATCAGCCCCGTACCCATGc     >  1:1216419/1‑150 (MQ=255)
                                                                                                                                                     ggCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCGGGCATCCCCATATCGCCAATCAGCAcc               >  1:1599789/1‑137 (MQ=255)
                                                                                                                                                     ggCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCGGGCATCCCCATATCGCCAATCAGCAcc               <  2:1599789/137‑1 (MQ=255)
                                                                                                                                                     ggCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCAGGACCGCCGTTGAAACGCCCACCAGCGGGCATCCCCATATCGCCAATCAGCACCGTACCCATGCCCa  <  2:197395/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAAAGTTGTTGTTTAATCGCCGGATTAACAAAAATGCGCGTGACGTCTTTGTCCTGGGCGGCGAGTTTGATCAAGCTGAAAATTTCTGGCTTCCACAGCGTGGAGACAACGTGTTTACCGTCGCGGGAAACTAAGTCCAGTGCTTGCGGGCGCAAGAGCTGCGCGGAGGTCCAGCGAGTTTTCGGCAGTTGCAGAAAGATATCGACATCCAGTCCGGTCTGGTGGCTGGCATGACCGCCGTTGAAACGCCCACCAGCGGGCATCCCCATATCGCCAATCAGCACCGTACCCATGCCCA  >  NC_000913/2445824‑2446122

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