Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 578,172 G→T 42.2% A22A (GCG→GCT rzpD → DLP12 prophage; putative murein endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913578,1720GT42.2% 44.1 / 72.0 70A22A (GCG→GCTrzpDDLP12 prophage; putative murein endopeptidase
Reads supporting (aligned to +/- strand):  ref base G (20/19);  new base T (16/13);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 8.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.16e-02

TGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGT  >  NC_000913/578024‑578320
                                                                                                                                                    |                                                                                                                                                    
tGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCg                                                                                                                                                      >  8:352744/1‑149 (MQ=255)
   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                                                                                         <  4:311728/126‑1 (MQ=11)
                                                                                    tGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCagag                                                                                         >  3:311728/1‑126 (MQ=11)
                                                                                    tGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCtg                                                                              <  4:230302/137‑2 (MQ=11)
                                                                                    tGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCtg                                                                              >  3:230302/1‑136 (MQ=11)
                                                                                     gAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGa                                                                                   <  7:210068/131‑1 (MQ=11)
                                                                                     gAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGa                                                                                   >  8:210068/1‑131 (MQ=11)
                                                                                     gAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCtg                                                                              <  2:218203/136‑2 (MQ=11)
                                                                                      aGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGc                                                                                                                               >  3:222005/1‑86 (MQ=11)
                                                                                      aGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGc                                                                                                                               <  4:222005/86‑1 (MQ=11)
                                                                                      aGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCagag                                                                                         >  5:6182/1‑124 (MQ=11)
                                                                                      aGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCagag                                                                                         <  6:6182/124‑1 (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)
                                                                                                                                                    |                                                                                                                                                    
TGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGT  >  NC_000913/578024‑578320

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: