Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 95,070 G→A 5.6% 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.6% 287.8 / 5.1 90V141I (GTA→ATA) ‡murFUDP‑N‑acetylmuramoyl‑tripeptide:D‑alanyl‑D‑ alanine ligase
Reads supporting (aligned to +/- strand):  ref base G (42/42);  new base A (2/3);  total (44/46)
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.93e-01

TGCATGGGTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCA  >  NC_000913/94922‑95212
                                                                                                                                                    |                                                                                                                                              
tGCATGGGTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGgtg                                                                                                                                                >  1:225971/1‑149 (MQ=255)
 gCATGGGTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGgtgt                                                                                                                                               >  2:477675/1‑149 (MQ=255)
   aTGGGTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTAc                                                                                                                                             >  4:93042/1‑149 (MQ=255)
     gggTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCg                                                                                                                                           <  3:104065/149‑1 (MQ=255)
      ggTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGa                                                                                                                                          >  3:428790/1‑149 (MQ=255)
       gTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGAt                                                                                                                                         <  1:218430/149‑1 (MQ=255)
       gTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGAt                                                                                                                                         <  4:347686/149‑1 (MQ=255)
       gTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGAt                                                                                                                                         >  3:948413/1‑149 (MQ=255)
         tCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGa                                                                                                                                       >  2:983830/1‑149 (MQ=255)
         tCGCCAGCAAGTTCCGGCGAGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAGCCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGAAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGa                                                                                                                                       >  4:889553/1‑149 (MQ=255)
          cGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGAc                                                                                                                                      <  3:199924/149‑1 (MQ=255)
           gCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTACACCGATGTCg                                                                                                                                     <  4:573931/149‑1 (MQ=255)
            ccAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAGCGTCATCGGTACACCGATgtcgt                                                                                                                                    >  3:565260/1‑145 (MQ=255)
            ccAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGc                                                                                                                                    >  3:541192/1‑149 (MQ=255)
             cAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAGCGTCATCGGTACACCGATgtcgtt                                                                                                                                   <  4:565260/149‑6 (MQ=255)
             cAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCt                                                                                                                                   >  3:541738/1‑149 (MQ=255)
             cAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTACACCGATGTCgtt                                                                                                                                   <  4:11455/149‑3 (MQ=255)
              aGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTg                                                                                                                                  >  2:692497/1‑149 (MQ=255)
                       cGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAGAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTAGACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTgcgc                                                                                                                            >  4:998903/1‑146 (MQ=255)
                       cGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTgcgc                                                                                                                            <  3:998903/146‑1 (MQ=255)
                         gcgcgcGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAAc                                                                                                                       <  1:859269/149‑1 (MQ=255)
                           gcgcgTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGc                                                                                                                     <  4:82113/149‑1 (MQ=255)
                              cgTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCgg                                                                                                                  <  3:72126/149‑1 (MQ=255)
                               gTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGa                                                                                                                 >  1:421714/1‑149 (MQ=255)
                                   ttGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATAc                                                                                                             >  2:18448/1‑149 (MQ=255)
                                     gCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGa                                                                                                           <  1:251990/149‑1 (MQ=255)
                                       tctGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGAtt                                                                                                         >  3:603761/1‑149 (MQ=255)
                                       tctGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGAtt                                                                                                         >  2:62272/1‑149 (MQ=255)
                                         tGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTAc                                                                                                       <  1:910957/149‑1 (MQ=255)
                                          gACGGGGTCCTCCGGCAAAAGCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCCGGCAATCTCAACAACGACATCGCTGTACCGATGACGATGTTGCGCTTAACGCCGGAATACGATTACg                                                                                                      >  1:850060/1‑149 (MQ=255)
                                             ggggTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAg                                                                                                   <  2:526222/149‑1 (MQ=255)
                                              gggTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAGGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGt                                                                                                  <  2:346711/149‑1 (MQ=255)
                                                 tcctccGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTAt                                                                                               <  3:474019/149‑1 (MQ=255)
                                                   ctccGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTg                                                                                             <  3:575465/149‑1 (MQ=255)
                                                     ccGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGaa                                                                                           <  3:649684/149‑1 (MQ=255)
                                                     ccGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGaa                                                                                           <  4:283291/149‑1 (MQ=255)
                                                      cGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAAc                                                                                          >  4:1040372/1‑149 (MQ=255)
                                                        gCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACtt                                                                                        >  2:842071/1‑149 (MQ=255)
                                                         cAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTg                                                                                       <  4:450419/149‑1 (MQ=255)
                                                               cTCCGTTAAAGAGATGACGGCGGCGATGTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGa                                                                                 >  2:463061/1‑149 (MQ=255)
                                                                tCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGaa                                                                                >  3:934536/1‑149 (MQ=255)
                                                                       aaGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCtt                                                                                                                          <  4:1018942/100‑1 (MQ=255)
                                                                       aaGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCtt                                                                                                                          >  3:1018942/1‑100 (MQ=255)
                                                                       aaGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAg                                                                         >  3:553102/1‑149 (MQ=255)
                                                                        agagATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGTGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAgg                                                                        >  3:678790/1‑149 (MQ=255)
                                                                        agagATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGACCTCGGCGCGAACCATCAgg                                                                        <  4:541738/149‑1 (MQ=255)
                                                                        agagATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAgg                                                                        <  2:81501/149‑1 (MQ=255)
                                                                        agagATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAgg                                                                        >  4:371740/1‑149 (MQ=255)
                                                                        agagATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATCACGCAGTTATTGAACTTGGCGCGAACCCTCAgg                                                                        <  1:304698/149‑1 (MQ=255)
                                                                         gagaTGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAggg                                                                       >  1:920774/1‑149 (MQ=255)
                                                                         gagaTGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAggg                                                                       >  4:699574/1‑149 (MQ=255)
                                                                           gaTGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCg                                                                     <  3:1042813/149‑1 (MQ=255)
                                                                              gACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACGACGACATCGGTCTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGaaa                                                                  <  4:959051/149‑1 (MQ=255)
                                                                              gACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGaaa                                                                  <  2:730745/149‑1 (MQ=255)
                                                                                    ggcgATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCt                                                            <  3:93042/149‑1 (MQ=255)
                                                                                       gATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGa                                                         >  4:640188/1‑149 (MQ=255)
                                                                                       gATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGa                                                         <  1:463061/149‑1 (MQ=255)
                                                                                         ttttAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACt                                                       <  4:938706/149‑1 (MQ=255)
                                                                                           ttAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACtgt                                                     <  1:357520/149‑1 (MQ=255)
                                                                                             aaGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGa                                                   <  2:790573/149‑1 (MQ=255)
                                                                                             aaGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGa                                                   <  2:965973/149‑1 (MQ=255)
                                                                                                ccAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGa                                                   <  2:921298/146‑1 (MQ=255)
                                                                                                ccAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGa                                                   >  1:921298/1‑146 (MQ=255)
                                                                                                 cAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCt                                               <  4:453841/149‑1 (MQ=255)
                                                                                                   gTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGa                                             >  1:20578/1‑149 (MQ=255)
                                                                                                          aaCACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGAc                                                        >  3:1013559/1‑131 (MQ=255)
                                                                                                          aaCACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGAc                                                        <  4:1013559/131‑1 (MQ=255)
                                                                                                           acacGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTAGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCg                                     >  2:872371/1‑149 (MQ=255)
                                                                                                            ctggCTTTATGCGGCAGGCAATCTCAACAACGGCATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCgg                                    >  2:24825/4‑149 (MQ=255)
                                                                                                             acGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGa                                   >  3:171421/1‑149 (MQ=255)
                                                                                                              cGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGccc                                      >  3:53384/1‑145 (MQ=255)
                                                                                                              cGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGccc                                      <  4:53384/145‑1 (MQ=255)
                                                                                                              cGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGaa                                  <  3:118446/149‑1 (MQ=255)
                                                                                                                 tttATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTAt                                                                                               >  3:928290/1‑85 (MQ=255)
                                                                                                                 tttATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTAt                                                                                               <  4:928290/85‑1 (MQ=255)
                                                                                                                 tttATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCt                               >  1:200894/1‑149 (MQ=255)
                                                                                                                 tttATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCt                               >  1:829125/1‑149 (MQ=255)
                                                                                                                   tataCGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTgc                             <  4:30968/149‑1 (MQ=255)
                                                                                                                      aCGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCt                          >  1:802726/1‑149 (MQ=255)
                                                                                                                       cGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTg                         <  1:1051023/149‑1 (MQ=255)
                                                                                                                       cGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTg                         <  4:920214/149‑1 (MQ=255)
                                                                                                                             gcaATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTcaac                   <  2:798104/149‑1 (MQ=255)
                                                                                                                                 tctcAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCTCGAACCAGCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAGcc               <  2:421714/149‑1 (MQ=255)
                                                                                                                                 tctcAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAAcc               >  2:630679/1‑149 (MQ=255)
                                                                                                                                 tctcAACAACGACATCGGTGTACCGATGACGCTGCTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGAGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAAcc               >  3:632297/1‑149 (MQ=255)
                                                                                                                                    caacaaCGACATCGGTGTACCGATAGCGCTGTTGCGCTTAACGCCAGAATACGATTACCCAGGTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTgg            >  2:426918/1‑149 (MQ=255)
                                                                                                                                        aacGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCGGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGGAgc        >  2:475269/1‑149 (MQ=255)
                                                                                                                                         acGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAgcg       <  1:716370/149‑1 (MQ=255)
                                                                                                                                           gtcATCGGTACACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAgcggc     >  1:622274/3‑149 (MQ=255)
                                                                                                                                           gACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGaa                                  <  4:967449/120‑1 (MQ=255)
                                                                                                                                           gACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGaa                                  >  3:967449/1‑120 (MQ=255)
                                                                                                                                              aTCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCa  >  1:787486/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
TGCATGGGTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCA  >  NC_000913/94922‑95212

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