Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 578,175 C→T 42.2% V23V (GTC→GTT rzpD → DLP12 prophage; putative murein endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913578,1750CT42.2% 55.7 / 61.4 72V23V (GTC→GTTrzpDDLP12 prophage; putative murein endopeptidase
Reads supporting (aligned to +/- strand):  ref base C (21/20);  new base T (16/14);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-02

ACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTA  >  NC_000913/578027‑578322
                                                                                                                                                    |                                                                                                                                                   
aCATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTc                                                                                                                                                     <  2:9859/149‑1 (MQ=255)
  aTACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGCCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCaa                                                                                                                                                   >  7:103696/1‑149 (MQ=255)
   tACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTTTCATGGGCGGTCAAt                                                                                                                                                  >  1:199497/1‑149 (MQ=255)
             gTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGt                                                                                                                                         >  5:123472/1‑148 (MQ=255)
                  cAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAAc                                                                                                                                   >  5:133371/1‑149 (MQ=255)
                       gAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAAt                                                                                                                              <  8:274063/149‑1 (MQ=255)
                         aaGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCAGTCAATCATTACCGTGATAACGCAATCg                                                                                                                            >  7:220323/1‑149 (MQ=255)
                            gCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATAATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGcc                                                                                                                          >  4:317541/1‑148 (MQ=255)
                                gatgaCTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGAATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACa                                                                                                                      >  6:187840/1‑148 (MQ=255)
                                 atgaCTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACaaa                                                                                                                    <  6:123472/149‑1 (MQ=255)
                                        gtcgtGAGATTGAGCGTGAAGTCTCTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGATCAATAATTAACTTGATAAATCAATAGCCTACAAAGTCAAgc                                                                                                             >  7:353498/1‑149 (MQ=255)
                                                 ttGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCa                                                                                                                                                <  5:319191/105‑1 (MQ=255)
                                                 ttGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCa                                                                                                                                                >  6:319191/1‑105 (MQ=255)
                                                 ttGAGCGTGAAGTCTGTTTGTGCCGCCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTCATTATCTGCATCATCGTCAGCCTGTCATGGGCGGGCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCACCGCGACaaa                                                                                                     >  5:4801/1‑148 (MQ=37)
                                                        tGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGc                                                                                                >  3:255292/1‑146 (MQ=255)
                                                          aaGTCTGTTTGTTGGGGCAACAGGGAGCAGAGTCACCGCGATTATATACGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCaga                                                                                            <  6:105064/148‑1 (MQ=255)
                                                            gtattgttGTGGGGGCAACAGTGAGCAGAGTAACCGAGAGTATATCCGCTCTGATTATCTGCATCATCTTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGaa                                                                                         <  2:199497/143‑1 (MQ=25)
                                                                       ggggCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGc                                                                              <  4:255292/149‑1 (MQ=255)
                                                                          gCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGaa                                                                           <  5:148594/149‑1 (MQ=255)
                                                                           gatgtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGCTTATCTGCATCATCGTCTGCCTGTCATGGTCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAAc                                                                          >  2:272482/7‑149 (MQ=11)
                                                                            atgtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGc                                                                                                <  3:203957/121‑1 (MQ=11)
                                                                            atgtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCa                                                                                              <  3:315251/123‑1 (MQ=11)
                                                                            atgtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCa                                                                                              >  4:315251/6‑128 (MQ=11)
                                                                             tgtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCCTCGTCTGCCTGTCATGGGCGGTTTATCATTCCCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAAcgc                                                                        >  3:300436/5‑149 (MQ=11)
                                                                              gtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCaga                                                                                            >  1:183910/4‑128 (MQ=11)
                                                                              gtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCaga                                                                                            <  2:183910/125‑1 (MQ=11)
                                                                              gtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAAc                                                                                        <  1:217212/129‑1 (MQ=11)
                                                                              gtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAAc                                                                                        >  2:217212/4‑132 (MQ=11)
                                                                              cAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAAcgcg                                                                       >  7:344002/1‑149 (MQ=255)
                                                                               tatGAGCAGAGTCACCGCGATTATCTCCGCTCTGTTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATACCGGCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAAcgcg                                                                       <  6:191777/146‑1 (MQ=11)
                                                                               tatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCa                                                                                              <  5:105395/123‑1 (MQ=11)
                                                                               tatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCa                                                                                              >  6:105395/3‑125 (MQ=11)
                                                                                gTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCgg                                                                      >  6:174608/1‑148 (MQ=255)
                                                                                 tGAGCAGAGTCACCGCGATTATCTCCGCTCTGTTTATCTGCATCATCGTCTGCCTGTCATTGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGc                                                                     <  6:109236/148‑1 (MQ=11)
                                                                                 tGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCagag                                                                                           >  3:311728/1‑126 (MQ=11)
                                                                                 tGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCagag                                                                                           <  4:311728/126‑1 (MQ=11)
                                                                                 tGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCtg                                                                                <  4:230302/137‑2 (MQ=11)
                                                                                 tGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCtg                                                                                >  3:230302/1‑136 (MQ=11)
                                                                                  gAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGa                                                                                     >  8:210068/1‑131 (MQ=11)
                                                                                  gAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGa                                                                                     <  7:210068/131‑1 (MQ=11)
                                                                                  gAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCtg                                                                                <  2:218203/136‑2 (MQ=11)
                                                                                   aGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGc                                                                                                                                 >  3:222005/1‑86 (MQ=11)
                                                                                   aGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGc                                                                                                                                 <  4:222005/86‑1 (MQ=11)
                                                                                   aGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCagag                                                                                           <  6:6182/124‑1 (MQ=11)
                                                                                   aGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCagag                                                                                           >  5:6182/1‑124 (MQ=11)
                                                                                   aGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCAGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAAt                                                                  <  8:220323/149‑1 (MQ=255)
                                                                                    gCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAgcgc                                                                                                           >  8:65908/1‑107 (MQ=11)
                                                                                      agagTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCTGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGc                                                                              <  6:91822/134‑1 (MQ=255)
                                                                                      agagTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGc                                                                              >  5:91822/1‑134 (MQ=255)
                                                                                             ccGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATg                                                        >  2:96416/1‑149 (MQ=11)
                                                                                              cgcgATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATg                                                        >  8:169899/1‑148 (MQ=255)
                                                                                                       ctCCGATCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATGATTACCGAGATAACGACATTACCTACAACGCCCAGCGCGCCAAAAATGCAAGAGAACTGAAGCTGGAGCACGCGGCAATTACTGACATGAAGATGCGt                                               >  6:236726/2‑148 (MQ=255)
                                                                                                           gCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCTAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCg                                          <  8:103696/149‑1 (MQ=255)
                                                                                                               tGATTATCTGCATCATCGTCACCCTGTCATGGGCGGTCAATCATTACCGTGATAACGGATTCGCCTACAATGTCCAGCGCGACAAAAATGCCAGAGATCTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGAt                                      <  6:4801/149‑1 (MQ=17)
                                                                                                                 gtgaTCTTCATCATCGTCTGCCTGTCATGGGCGGGTAATCATTACCGTGATAACGCCATTACCTAAAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGt                                    <  4:188766/146‑1 (MQ=11)
                                                                                                                 aTTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGt                                    >  1:106550/1‑149 (MQ=255)
                                                                                                                         catcatCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgct                                <  5:174608/145‑1 (MQ=255)
                                                                                                                         catcatCGTCAGCATGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTgcgc                            >  8:153312/1‑149 (MQ=255)
                                                                                                                          atcatcGTCAGCCTGTCATGGGCTGTAAATCATTACCGTGATAAAGCCATTACATACAAAGCCCATCGCGACAAAAATGCAAGAGAACTGAACCTGGCTAACACGGCAATAACTGACATGCAGATGCGTCAGAGTGATGTTGCTGCGCt                           >  8:249689/1‑149 (MQ=255)
                                                                                                                           tcatcGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTc                          >  7:305095/1‑149 (MQ=255)
                                                                                                                            catcGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCACAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCTGCGTGATGTTGCTGCGCTc                          >  6:382667/1‑148 (MQ=11)
                                                                                                                            catcGTCTGCCTGTCATGGGCTGTTAATCATTAACGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGATAACTGAAGCTGTCTAACTAGGCAATTACTGACATTCAGATGAGTCAGAGTGATGTTTCTGAGCTc                          >  6:67882/1‑148 (MQ=255)
                                                                                                                             atcGGCAGCCTGTCATGGGCGGTCTATCTTTACCGTGATCTGTCATTCGCGTACTTAGTCCAGCGCGTGATATATGCCAGTGAACTGATGCTAGCGATCGGGGCTATTACTGACATGCATATGCGTCAGCGTGATGTTGCTGCGCTCgt                        <  6:71059/149‑2 (MQ=12)
                                                                                                                              tcGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGAt                       <  6:133371/149‑1 (MQ=255)
                                                                                                                                    gCCTGTCATGGGCGGTCAATCATTAGCGTGATAACGCAATCGCCTACAAAGTCCTGCGCGACAAAACTGCCAGAGCACTGAAGCTAGCGATCGCGGCAATTACTGACCTGCAGATGCGGCAGCGTGATGTTGCTGCGCTCGATGCaaaa                 <  2:278439/149‑1 (MQ=17)
                                                                                                                                         tCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAg                                                     <  5:347098/108‑1 (MQ=255)
                                                                                                                                         tCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAg                                                     >  6:347098/1‑108 (MQ=255)
                                                                                                                                         tCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATac              <  5:187840/147‑1 (MQ=255)
                                                                                                                                              ggCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAg        <  2:106550/148‑1 (MQ=255)
                                                                                                                                               gCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAg        >  2:314392/1‑147 (MQ=255)
                                                                                                                                                 ggTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGAGACAAAAATGCCAGAGACCTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGt    <  5:341392/149‑1 (MQ=255)
                                                                                                                                                   tCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTa  <  7:352744/149‑1 (MQ=255)
                                                                                                                                                   tCAATCATTACCGTGATAACGCAATCGACTACAAAGTCCAGCGCGAAAAAAATGCCAGAGAACTGAAGATAGCGAACGCGGCAATTACTGACATGAAGATGCGTCAGCTTGATGTTGCTGCTCTCGATGCAAAATACACGaa         >  8:217084/1‑142 (MQ=14)
                                                                                                                                                    |                                                                                                                                                   
ACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTA  >  NC_000913/578027‑578322

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: