Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 95,070 G→A 5.2% 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.2% 391.2 / 7.7 117V141I (GTA→ATA) ‡murFUDP‑N‑acetylmuramoyl‑tripeptide:D‑alanyl‑D‑ alanine ligase
Reads supporting (aligned to +/- strand):  ref base G (56/54);  new base A (4/2);  total (61/56)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.29e-01

Only 100 of 120 total aligned reads displayed.

GCATGGGTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCT  >  NC_000913/94923‑95215
                                                                                                                                                   |                                                                                                                                                 
  aTGGGTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTAc                                                                                                                                                >  1:260879/1‑149 (MQ=255)
  aTGGGTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTAc                                                                                                                                                <  2:676135/149‑1 (MQ=255)
   tGGGTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTAc                                                                                                                                                >  3:56492/1‑148 (MQ=255)
      gTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGAt                                                                                                                                            <  2:652293/149‑1 (MQ=255)
       ttCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGAt                                                                                                                                            >  3:221704/1‑148 (MQ=255)
           ccAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTACACCGATGTCgt                                                                                                                                       <  2:992211/149‑2 (MQ=255)
                aaGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTg                                                                                                                                  <  2:774011/149‑1 (MQ=255)
                 aGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTg                                                                                                                                  >  3:816944/1‑148 (MQ=255)
                   ttCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGAt                                                                                                                                            <  1:804315/136‑1 (MQ=255)
                   ttCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGAt                                                                                                                                            >  2:804315/1‑136 (MQ=255)
                       ggCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTa                                                                                                                            <  3:754677/148‑1 (MQ=255)
                        gcgcgcGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACTGATGACGCTGTTGCGCTTaa                                                                                                                           >  4:743762/1‑148 (MQ=255)
                        gcgcgcGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTaa                                                                                                                           >  4:493502/1‑148 (MQ=255)
                          gcgcgTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACg                                                                                                                         <  3:947047/148‑1 (MQ=255)
                              gTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGa                                                                                                                    >  2:669645/1‑149 (MQ=255)
                               tGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGa                                                                                                                    <  4:40666/148‑1 (MQ=255)
                                ggTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGaa                                                                                                                   >  3:124824/1‑148 (MQ=255)
                                      tctGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGAtt                                                                                                            >  2:881788/1‑149 (MQ=255)
                                      tctGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTACACCGATg                                                                                                                                           >  2:993188/1‑118 (MQ=255)
                                        tGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTa                                                                                                           <  4:816944/148‑1 (MQ=255)
                                              ggTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCa                                                                                                       <  2:761741/146‑1 (MQ=255)
                                              ggTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCa                                                                                                       >  1:761741/1‑146 (MQ=255)
                                                 cctccGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTAt                                                                                                  <  4:814606/148‑1 (MQ=255)
                                                    ccGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGa                                                                                               >  3:766214/1‑148 (MQ=255)
                                                       gCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGCCAACACGCTTTTAACGGCAGCCATCCTCACCAACGCCATCGGTTTACCGATGACGCGGTTGCGCTTACCGCGGAATTACGATTAGGCAGTTATTAAACtt                                                                                           >  1:981967/1‑149 (MQ=255)
                                                        cAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACtt                                                                                           <  4:124824/148‑1 (MQ=255)
                                                        cAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTg                                                                                          >  2:970174/1‑149 (MQ=255)
                                                          aaaCCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGcc                                                                                                                       >  1:655060/1‑118 (MQ=255)
                                                          aaaCCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGa                                                                                                              <  1:342548/127‑1 (MQ=255)
                                                          aaaCCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGa                                                                                                              >  2:342548/1‑127 (MQ=255)
                                                          aaaCCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGgc                                                                                        <  1:637946/149‑1 (MQ=255)
                                                          aaaCCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGgc                                                                                        >  1:39777/1‑149 (MQ=255)
                                                               tCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGaa                                                                                   >  2:72553/1‑149 (MQ=255)
                                                               tCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGaa                                                                                   >  1:469730/1‑149 (MQ=255)
                                                                  gTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCa                                                                                >  1:983078/1‑149 (MQ=255)
                                                                    tAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCAt                                                                               <  4:411776/148‑1 (MQ=255)
                                                                    tAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCAt                                                                               >  3:943792/1‑148 (MQ=255)
                                                                    tAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCAt                                                                               >  3:396215/1‑148 (MQ=255)
                                                                     aaaGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCa                                                                             >  1:24791/1‑149 (MQ=255)
                                                                      aaGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCa                                                                             <  3:146955/148‑1 (MQ=255)
                                                                      aaGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAg                                                                            >  2:184501/1‑149 (MQ=255)
                                                                       agagATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAgg                                                                           <  2:29227/149‑1 (MQ=255)
                                                                       agagATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAg                                                                            >  3:960323/1‑148 (MQ=255)
                                                                        gagaTGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAggg                                                                          <  1:824578/149‑1 (MQ=255)
                                                                         agaTGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGc                                                                         >  1:529828/1‑149 (MQ=255)
                                                                          gaTGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGc                                                                         >  4:168505/1‑148 (MQ=255)
                                                                               cggcggCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAAt                                                                    <  3:3455/148‑1 (MQ=255)
                                                                               cggcggCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATa                                                                   <  2:976249/149‑1 (MQ=255)
                                                                                 gcggcgATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAg                                                                  <  3:550354/148‑1 (MQ=255)
                                                                                 gcggcgATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAg                                                                  <  3:355561/148‑1 (MQ=255)
                                                                                 gcggcgATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAg                                                                  <  3:352905/148‑1 (MQ=255)
                                                                                 gcggcgATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGc                                                                 <  2:287606/149‑1 (MQ=255)
                                                                                  cggcgATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGcc                                                                <  1:335014/149‑1 (MQ=255)
                                                                                    gcgATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCt                                                               >  3:805962/1‑148 (MQ=255)
                                                                                    gcgATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTg                                                              <  1:773124/149‑1 (MQ=255)
                                                                                     cgAATTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTgg                                                             <  2:981967/149‑1 (MQ=255)
                                                                                      gATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGa                                                            >  2:2748/1‑149 (MQ=255)
                                                                                        ttttAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACt                                                          >  1:70627/1‑149 (MQ=255)
                                                                                         tttAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACtg                                                         >  2:434486/1‑149 (MQ=255)
                                                                                              gCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGt                                                    <  1:57857/149‑1 (MQ=255)
                                                                                              gCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGt                                                    <  2:1028562/149‑1 (MQ=255)
                                                                                              gCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGt                                                    >  1:534257/1‑149 (MQ=255)
                                                                                                cAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTc                                                   >  3:993075/1‑148 (MQ=255)
                                                                                                cAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCt                                                  <  2:512105/149‑1 (MQ=255)
                                                                                                       gCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCg                                            >  4:300058/1‑148 (MQ=255)
                                                                                                               cTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTg                                                                                          <  3:521893/94‑1 (MQ=255)
                                                                                                               cTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAAt                                                                    >  3:917153/1‑116 (MQ=255)
                                                                                                               cTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAAt                                                                    <  4:917153/116‑1 (MQ=255)
                                                                                                                     aCGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTgcgc                              >  4:585861/1‑148 (MQ=255)
                                                                                                                      cGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTg                            <  1:184501/149‑1 (MQ=255)
                                                                                                                         caggcaATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGt                          <  4:221704/148‑1 (MQ=255)
                                                                                                                          aggcaATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTca                        >  2:301192/1‑149 (MQ=255)
                                                                                                                           ggcaATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTcaa                       <  2:24791/149‑1 (MQ=255)
                                                                                                                           ggcaATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTca                        >  4:623771/1‑148 (MQ=255)
                                                                                                                           ggcaATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATAAGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTca                        >  4:139564/1‑148 (MQ=255)
                                                                                                                            gcaATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTcaac                      <  1:134922/149‑1 (MQ=255)
                                                                                                                            gcaATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTcaac                      <  1:775937/149‑1 (MQ=255)
                                                                                                                                tctcAACAACGACATCGGTGTACTGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCaacaac                   <  3:743762/148‑1 (MQ=255)
                                                                                                                                tctcAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCt                                                               <  3:284213/104‑1 (MQ=255)
                                                                                                                                tctcAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCt                                                               >  4:284213/1‑104 (MQ=255)
                                                                                                                                 cgcaacagCGTCATCGGTACACCGATGTCGTTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAAcc                  >  4:837793/3‑148 (MQ=38)
                                                                                                                                  tcAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTg                >  1:140395/1‑149 (MQ=255)
                                                                                                                                   caacaaCGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTg                >  3:632480/1‑148 (MQ=255)
                                                                                                                                   caacaaCGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTg                >  4:572306/1‑148 (MQ=255)
                                                                                                                                   caacaaCGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTg                >  3:330891/1‑148 (MQ=255)
                                                                                                                                    aacaacGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTgg               <  3:623771/148‑1 (MQ=255)
                                                                                                                                    aacaacGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGc              >  2:637487/1‑149 (MQ=255)
                                                                                                                                     acaacGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGc              >  4:260107/1‑148 (MQ=255)
                                                                                                                                     acaacGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCa             >  2:197036/1‑149 (MQ=255)
                                                                                                                                      caacGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCa             >  4:887010/1‑148 (MQ=255)
                                                                                                                                       aacGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAgc           >  2:774138/1‑149 (MQ=255)
                                                                                                                                          gACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAgcgg         >  2:786626/1‑148 (MQ=255)
                                                                                                                                           tcATCGGTACACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAgcggcg       >  1:940345/2‑149 (MQ=255)
                                                                                                                                           aCATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAgcggc        >  4:530077/1‑148 (MQ=255)
                                                                                                                                            cATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGgcgc      >  2:220380/1‑149 (MQ=255)
                                                                                                                                              tCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACtgtg                                                       <  3:263634/98‑1 (MQ=255)
                                                                                                                                                ggTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCt  <  2:529828/149‑1 (MQ=255)
                                                                                                                                                ggTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCt  <  1:434486/149‑1 (MQ=255)
                                                                                                                                                 gtgtACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGa                                                >  1:765360/1‑102 (MQ=255)
                                                                                                                                                 gtgtACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGa                                                <  2:765360/102‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCATGGGTTCGCCAGCAAGTTCCGGCGCGCGTGGTTGCTCTGACGGGGTCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCT  >  NC_000913/94923‑95215

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: