Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 578,200 G→A 42.6% A32T (GCC→ACC)  rzpD → DLP12 prophage; putative murein endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913578,2000GA42.6% 57.8 / 66.2 75A32T (GCC→ACC) rzpDDLP12 prophage; putative murein endopeptidase
Reads supporting (aligned to +/- strand):  ref base G (20/23);  new base A (18/14);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-01

AAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGA  >  NC_000913/578052‑578345
                                                                                                                                                     |                                                                                                                                                 
aaGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCAGTCAATCATTACCGTGATAACG‑CAATCg                                                                                                                                                   >  7:220323/1‑149 (MQ=255)
   gCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATAATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGcc                                                                                                                                                 >  4:317541/1‑148 (MQ=255)
       gatgaCTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGAATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACa                                                                                                                                             >  6:187840/1‑148 (MQ=255)
        atgaCTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACaaa                                                                                                                                           <  6:123472/149‑1 (MQ=255)
               gtcgtGAGATTGAGCGTGAAGTCTCTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGATCAATAATTAACTTGATAAAT‑CAATAGCCTACAAAGTCAAgc                                                                                                                                    >  7:353498/1‑149 (MQ=255)
                        ttGAGCGTGAAGTCTGTTTGTGCCGCCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTCATTATCTGCATCATCGTCAGCCTGTCATGGGCGGGCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCACCGCGACaaa                                                                                                                            >  5:4801/1‑148 (MQ=37)
                               tGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGc                                                                                                                       >  3:255292/1‑146 (MQ=255)
                                 aaGTCTGTTTGTTGGGGCAACAGGGAGCAGAGTCACCGCGATTATATACGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCaga                                                                                                                   <  6:105064/148‑1 (MQ=255)
                                   gtattgttGTGGGGGCAACAGTGAGCAGAGTAACCGAGAGTATATCCGCTCTGATTATCTGCATCATCTTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGaa                                                                                                                <  2:199497/143‑1 (MQ=25)
                                              ggggCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGc                                                                                                     <  4:255292/149‑1 (MQ=255)
                                                 gCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGaa                                                                                                  <  5:148594/149‑1 (MQ=255)
                                                  gatgtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGCTTATCTGCATCATCGTCTGCCTGTCATGGTCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAAc                                                                                                 >  2:272482/7‑149 (MQ=11)
                                                   atgtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGc                                                                                                                       <  3:203957/121‑1 (MQ=11)
                                                   atgtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCa                                                                                                                     >  4:315251/6‑128 (MQ=11)
                                                   atgtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCa                                                                                                                     <  3:315251/123‑1 (MQ=11)
                                                    tgtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCCTCGTCTGCCTGTCATGGGCGGTTTATCATTCCCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAAcgc                                                                                               >  3:300436/5‑149 (MQ=11)
                                                     gtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCaga                                                                                                                   >  1:183910/4‑128 (MQ=11)
                                                     gtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCaga                                                                                                                   <  2:183910/125‑1 (MQ=11)
                                                     gtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAAc                                                                                                               <  1:217212/129‑1 (MQ=11)
                                                     gtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAAc                                                                                                               >  2:217212/4‑132 (MQ=11)
                                                     cAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAAcgcg                                                                                              >  7:344002/1‑149 (MQ=255)
                                                      tatGAGCAGAGTCACCGCGATTATCTCCGCTCTGTTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATACCG‑GCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAAcgcg                                                                                              <  6:191777/146‑1 (MQ=11)
                                                      tatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCa                                                                                                                     >  6:105395/3‑125 (MQ=11)
                                                      tatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCa                                                                                                                     <  5:105395/123‑1 (MQ=11)
                                                       gTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCgg                                                                                             >  6:174608/1‑148 (MQ=255)
                                                        tGAGCAGAGTCACCGCGATTATCTCCGCTCTGTTTATCTGCATCATCGTCTGCCTGTCATTGGCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGc                                                                                            <  6:109236/148‑1 (MQ=11)
                                                        tGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCagag                                                                                                                  >  3:311728/1‑126 (MQ=11)
                                                        tGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCagag                                                                                                                  <  4:311728/126‑1 (MQ=11)
                                                        tGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCtg                                                                                                       <  4:230302/137‑2 (MQ=11)
                                                        tGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCtg                                                                                                       >  3:230302/1‑136 (MQ=11)
                                                         gAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGa                                                                                                            >  8:210068/1‑131 (MQ=11)
                                                         gAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGa                                                                                                            <  7:210068/131‑1 (MQ=11)
                                                         gAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCtg                                                                                                       <  2:218203/136‑2 (MQ=11)
                                                          aGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCagag                                                                                                                  >  5:6182/1‑124 (MQ=11)
                                                          aGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCagag                                                                                                                  <  6:6182/124‑1 (MQ=11)
                                                          aGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCAGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAAt                                                                                         <  8:220323/149‑1 (MQ=255)
                                                           gCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAgcgc                                                                                                                                  >  8:65908/1‑107 (MQ=11)
                                                             agagTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCTGCCTGTCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGc                                                                                                     <  6:91822/134‑1 (MQ=255)
                                                             agagTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGc                                                                                                     >  5:91822/1‑134 (MQ=255)
                                                                    ccGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATg                                                                               >  2:96416/1‑149 (MQ=11)
                                                                     cgcgATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATg                                                                               >  8:169899/1‑148 (MQ=255)
                                                                              ctCCGATCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATGATTACCGAGATAACG‑ACATTACCTACAACGCCCAGCGCGCCAAAAATGCAAGAGAACTGAAGCTGGAGCACGCGGCAATTACTGACATGAAGATGCGt                                                                      >  6:236726/2‑148 (MQ=255)
                                                                                  gCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCTAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCg                                                                 <  8:103696/149‑1 (MQ=255)
                                                                                      tGATTATCTGCATCATCGTCACCCTGTCATGGGCGGTCAATCATTACCGTGATAACG‑GATTCGCCTACAATGTCCAGCGCGACAAAAATGCCAGAGATCTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGAt                                                             <  6:4801/149‑1 (MQ=17)
                                                                                        gtgaTCTTCATCATCGTCTGCCTGTCATGGGCGGGTAATCATTACCGTGATAACG‑CCATTACCTAAAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGt                                                           <  4:188766/146‑1 (MQ=11)
                                                                                        aTTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGt                                                           >  1:106550/1‑149 (MQ=255)
                                                                                                catcatCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgct                                                       <  5:174608/145‑1 (MQ=255)
                                                                                                catcatCGTCAGCATGTCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTgcgc                                                   >  8:153312/1‑149 (MQ=255)
                                                                                                 atcatcGTCAGCCTGTCATGGGCTGTAAATCATTACCGTGATAAAG‑CCATTACATACAAAGCCCATCGCGACAAAAATGCAAGAGAACTGAACCTGGCTAACACGGCAATAACTGACATGCAGATGCGTCAGAGTGATGTTGCTGCGCt                                                  >  8:249689/1‑149 (MQ=255)
                                                                                                  tcatcGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTc                                                 >  7:305095/1‑149 (MQ=255)
                                                                                                   catcGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCACAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCTGCGTGATGTTGCTGCGCTc                                                 >  6:382667/1‑148 (MQ=11)
                                                                                                   catcGTCTGCCTGTCATGGGCTGTTAATCATTAACGTGATAACG‑CCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGATAACTGAAGCTGTCTAACTAGGCAATTACTGACATTCAGATGAGTCAGAGTGATGTTTCTGAGCTc                                                 >  6:67882/1‑148 (MQ=255)
                                                                                                    atcGGCAGCCTGTCATGGGCGGTCTATCTTTACCGTGATCTGT‑CATTCGCGTACTTAGTCCAGCGCGTGATATATGCCAGTGAACTGATGCTAGCGATCGGGGCTATTACTGACATGCATATGCGTCAGCGTGATGTTGCTGCGCTCgt                                               <  6:71059/149‑2 (MQ=12)
                                                                                                     tcGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGAt                                              <  6:133371/149‑1 (MQ=255)
                                                                                                           gCCTGTCATGGGCGGTCAATCATTAGCGTGATAACG‑CAATCGCCTACAAAGTCCTGCGCGACAAAACTGCCAGAGCACTGAAGCTAGCGATCGCGGCAATTACTGACCTGCAGATGCGGCAGCGTGATGTTGCTGCGCTCGATGCaaaa                                        <  2:278439/149‑1 (MQ=17)
                                                                                                                tCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAg                                                                            >  6:347098/1‑108 (MQ=255)
                                                                                                                tCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAg                                                                            <  5:347098/108‑1 (MQ=255)
                                                                                                                tCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATac                                     <  5:187840/147‑1 (MQ=255)
                                                                                                                     ggatgttaaTCATTACCGTGATAACG‑CAATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAAATGAAGCTTGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAg                               <  1:97848/141‑1 (MQ=11)
                                                                                                                     ggCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAg                               <  2:106550/148‑1 (MQ=255)
                                                                                                                      gCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAg                               >  2:314392/1‑147 (MQ=255)
                                                                                                                        ggTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGAGACAAAAATGCCAGAGACCTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGt                           <  5:341392/149‑1 (MQ=255)
                                                                                                                          tCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTa                         <  7:352744/149‑1 (MQ=255)
                                                                                                                          tCAATCATTACCGTGATAACG‑CAATCGACTACAAAGTCCAGCGCGAAAAAAATGCCAGAGAACTGAAGATAGCGAACGCGGCAATTACTGACATGAAGATGCGTCAGCTTGATGTTGCTGCTCTCGATGCAAAATACACGaa                                >  8:217084/1‑142 (MQ=14)
                                                                                                                                aTTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAATCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGAt                   <  3:317541/149‑1 (MQ=255)
                                                                                                                                aTTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGa                    <  7:169899/148‑1 (MQ=255)
                                                                                                                                 ttACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGAAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATg                  >  4:359584/1‑149 (MQ=255)
                                                                                                                                 ttACAGTGATAACGACAAT‑ACATACAAATCCCAGCGAGACAAAAAAGCCAGAGAACTGAAGCTGGCGAACGAGGCAATTACTGACAAGCAGATGCGTCAGAGAGAAGTTGCTGCGATAGATGAAAAATACACGAAGGAGTTAGCTGAt                   >  4:314430/1‑148 (MQ=2)
                                                                                                                                  tACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGc                 >  6:149577/1‑149 (MQ=255)
                                                                                                                                      gTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGaaa             <  8:305095/149‑1 (MQ=255)
                                                                                                                                       tGATAACG‑CCATTACATACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTCGCGAACGCGGCAATTACTGACATGAAGATGCGTCAGCGTGATCTTCCTGCCCTCGAAGCAAAATCCACGAAGCAGTAAGCTGATGATAAAg            >  4:283889/1‑149 (MQ=11)
                                                                                                                                         aTAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCt          >  2:165247/1‑149 (MQ=255)
                                                                                                                                             cG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTa                         <  1:303938/130‑1 (MQ=255)
                                                                                                                                             cG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTa                         >  2:303938/1‑130 (MQ=255)
                                                                                                                                              g‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGaaaa     >  7:225440/1‑149 (MQ=255)
                                                                                                                                                  aTCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAatga  <  7:153312/149‑1 (MQ=39)
                                                                                                                                                     |                                                                                                                                                 
AAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACG‑CAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGA  >  NC_000913/578052‑578345

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