Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 95,070 G→A 5.3% V141I (GTA→ATA) ‡ murF → UDP‑N‑acetylmuramoyl‑tripeptide:D‑alanyl‑D‑ alanine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091395,0700GA 5.3% 248.5 / 2.9 75V141I (GTA→ATA) ‡murFUDP‑N‑acetylmuramoyl‑tripeptide:D‑alanyl‑D‑ alanine ligase
Reads supporting (aligned to +/- strand):  ref base G (38/33);  new base A (2/2);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

ATGGGTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCA  >  NC_000913/94925‑95212
                                                                                                                                                 |                                                                                                                                              
aTGGGTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTAc                                                                                                                                             <  3:52318/149‑1 (MQ=255)
   ggTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGa                                                                                                                                          <  4:511992/149‑1 (MQ=255)
   ggTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGa                                                                                                                                          <  4:501330/149‑1 (MQ=255)
   ggTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGa                                                                                                                                          <  2:434693/149‑1 (MQ=255)
    gTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGAt                                                                                                                                         >  1:135729/1‑149 (MQ=255)
    gTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGAt                                                                                                                                         <  3:93203/149‑1 (MQ=255)
    gTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGAt                                                                                                                                         <  2:135729/149‑1 (MQ=255)
    gTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTACTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGAt                                                                                                                                         >  4:93203/1‑149 (MQ=255)
       cGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGAc                                                                                                                                      <  3:196911/149‑1 (MQ=255)
       cGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGAc                                                                                                                                      <  3:492718/149‑1 (MQ=255)
       cGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTACACCGATgtc                                                                                                                                      <  2:93474/149‑3 (MQ=255)
          cAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCt                                                                                                                                   >  4:241260/1‑149 (MQ=255)
           aGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTg                                                                                                                                  <  3:284228/149‑1 (MQ=255)
           aGCAAGTGCCGGCGCGCGTGGTTGCTCTGACGGGGTCTTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTg                                                                                                                                  >  2:554258/1‑149 (MQ=255)
            gCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGt                                                                                                                                 <  3:165334/149‑1 (MQ=255)
              aaGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTg                                                                                                                               >  2:488975/1‑149 (MQ=255)
              aaGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTg                                                                                                                               >  4:273014/1‑149 (MQ=255)
                 ttCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTgcgc                                                                                                                            >  2:499081/1‑149 (MQ=255)
                 ttCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTgcgc                                                                                                                            >  2:432578/1‑149 (MQ=255)
                   ccGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAGCGTCATCGGTACACCGATGTCGttgttgagatt                                                                                                                          <  4:298501/149‑6 (MQ=38)
                                   ctctGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGAt                                                                                                          >  1:320016/1‑149 (MQ=255)
                                    tctGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGAtt                                                                                                         >  4:21903/1‑149 (MQ=255)
                                      tGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTAc                                                                                                       >  3:469343/1‑149 (MQ=255)
                                      tGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTAc                                                                                                       <  3:241260/149‑1 (MQ=255)
                                          ggggTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAg                                                                                                   <  3:395981/149‑1 (MQ=255)
                                               cctccGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTAtt                                                                                              <  2:332290/149‑1 (MQ=255)
                                                ctctgacAAAACCTCCGTTAAAGAGATTACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTg                                                                                             <  4:490662/143‑1 (MQ=255)
                                                ctccGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTg                                                                                             >  4:295000/1‑149 (MQ=255)
                                                ctccGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTg                                                                                             >  2:21580/1‑149 (MQ=255)
                                                ctccGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTg                                                                                             <  2:81405/149‑1 (MQ=255)
                                                    ggCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACt                                                                                         >  4:123660/1‑149 (MQ=255)
                                                       aaaaCCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATGCGATTACGCAGTTATTGAACTTgg                                                                                      >  4:546757/1‑149 (MQ=255)
                                                       aaaaCCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTgg                                                                                      >  1:209745/1‑149 (MQ=255)
                                                            cTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGa                                                                                 >  1:103035/1‑149 (MQ=255)
                                                            cTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGa                                                                                 >  1:1379/1‑149 (MQ=255)
                                                               cGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAAcc                                                                              >  1:356949/1‑149 (MQ=255)
                                                                   aaaGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCa                                                                          <  1:488975/149‑1 (MQ=255)
                                                                       agaTGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGc                                                                      >  3:55285/1‑149 (MQ=255)
                                                                          tGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGaa                                                                   >  4:250329/1‑149 (MQ=255)
                                                                          tGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGAAGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGaa                                                                   >  4:397621/1‑149 (MQ=255)
                                                                           gACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGaaa                                                                  <  4:211483/149‑1 (MQ=255)
                                                                                cggcgATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGcc                                                             <  3:273014/149‑1 (MQ=255)
                                                                                cggcgATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGcc                                                             <  1:21580/149‑1 (MQ=255)
                                                                                  gcgATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTg                                                           <  2:581621/149‑1 (MQ=255)
                                                                                     aTTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGAc                                                        >  3:395781/1‑149 (MQ=255)
                                                                                      ttttAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATGGCCTGGACt                                                       >  3:427357/1‑149 (MQ=255)
                                                                                        ttAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACtgt                                                     >  1:140758/1‑149 (MQ=255)
                                                                                          aaGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGa                                                   <  4:2658/149‑1 (MQ=255)
                                                                                              cAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAAcc                                                                              <  2:340196/118‑1 (MQ=255)
                                                                                              cAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAAcc                                                                              >  1:340196/1‑118 (MQ=255)
                                                                                                 tGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGAc                                            <  2:488859/149‑1 (MQ=255)
                                                                                                 tGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGAc                                            >  1:199791/1‑149 (MQ=255)
                                                                                                       aaCACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGAc                                            <  3:510425/143‑1 (MQ=255)
                                                                                                       aaCACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGAc                                            >  4:510425/1‑143 (MQ=255)
                                                                                                       aaCACGCTTTATACGGCAGGCAATCTCAACAACGACATCGCTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTTTGACTCGccc                                      <  4:493202/149‑1 (MQ=255)
                                                                                                        acacGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCg                                     >  2:438265/1‑149 (MQ=255)
                                                                                                            gCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAg                                 >  2:274708/1‑149 (MQ=255)
                                                                                                              tttATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCt                               >  2:131234/1‑149 (MQ=255)
                                                                                                                  taCGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTgcgc                           <  1:438265/149‑1 (MQ=255)
                                                                                                                    cGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTg                         <  1:329179/149‑1 (MQ=255)
                                                                                                                    cGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTg                         >  3:519930/1‑149 (MQ=255)
                                                                                                                      gcaggcaATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTTTGACTCGCCCGGAAGCTGCGCTGGt                       <  4:469343/149‑1 (MQ=255)
                                                                                                                            aaTCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTcaacaa                 <  4:191736/149‑1 (MQ=255)
                                                                                                                              tctcAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAAcc               >  2:572777/1‑149 (MQ=255)
                                                                                                                              tctcAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAAcc               >  4:352925/1‑149 (MQ=255)
                                                                                                                               ctcAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCt              <  2:386451/149‑1 (MQ=255)
                                                                                                                                gcaacagCGTCATCGGTACACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTg             >  3:11829/2‑149 (MQ=255)
                                                                                                                                 caacaaCGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAAt                                                                                                               >  2:420241/1‑50 (MQ=255)
                                                                                                                                 caacaaCGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAAt                                                                                                               <  1:420241/50‑1 (MQ=255)
                                                                                                                                 caacaaCGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCTGAAGCTGCGCTGGTCAACAACCTgg            >  3:371443/1‑149 (MQ=255)
                                                                                                                                 caacaaCGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTgg            >  1:315011/1‑149 (MQ=255)
                                                                                                                                      acGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAgcg       <  2:320016/149‑1 (MQ=255)
                                                                                                                                          cATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGgcgc   <  2:199791/149‑1 (MQ=255)
                                                                                                                                          cATCGGTACACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGgcgc   >  3:264532/1‑149 (MQ=255)
                                                                                                                                           aTCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCa  >  4:166857/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
ATGGGTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCA  >  NC_000913/94925‑95212

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