Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,383,959 A→T 16.5% intergenic (+35/+36) elaD → / ← yfbK protease, capable of cleaving an AMC‑ubiquitin model substrate/Von Willebrand factor domain putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,383,9590AT16.5% 218.1 / 34.7 79intergenic (+35/+36)elaD/yfbKprotease, capable of cleaving an AMC‑ubiquitin model substrate/Von Willebrand factor domain putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (40/26);  new base T (8/5);  total (48/31)
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.65e-02

ATGCTGACAGATTTAAAAAATATTCAGCAGAAGAGCAGTACGAAATTGATTTTCAACATCGACTGGCAAACAGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCAGGAGGTATTGTTCAGGTATTCAGAACCGCGTAACTTTTGCCCATATGCAGCTACTGCTGCGCGAAAACGCATATCTTCAGAGGGCGCGTTTATTGTTG  >  NC_000913/2383740‑2384209
                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                          
aTGCTGACAGATTTAAAAAATATTCAGCAGAAGAGCAGTACGAAATTGATTTTCAACATCGACTGGCAAACAGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATgtcagtca                                                                                                                                                                                                     >  2:464955/1‑275 (MQ=255)
                                                                         aaaTTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATg                                                                                                                                                                                                                                                          <  8:823654/149‑1 (MQ=255)
                                                                          aaTTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATg                                                                                                                                                                                                                                                          <  5:847357/148‑1 (MQ=255)
                                                                             tGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGcc                                                                                                                        <  4:464953/275‑1 (MQ=255)
                                                                                            tatGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAg                                                                                                                                                                                                                                        >  5:622811/1‑148 (MQ=255)
                                                                                             atGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGt                                                                                                                                                                                                                                       <  6:622811/148‑1 (MQ=255)
                                                                                                             caatcaTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGtcatc                                                                                                                                                                                                                       >  6:520295/1‑148 (MQ=255)
                                                                                                                  attattATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGc                                                                                                                                                                                                                                                         <  7:222038/109‑1 (MQ=255)
                                                                                                                  attattATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGc                                                                                                                                                                                                                                                         >  8:222038/1‑109 (MQ=255)
                                                                                                                     attatAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATctgctg                                                                                                       >  1:134554/1‑252 (MQ=255)
                                                                                                                            aaaCTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCAGGAGGTATTGTTCAGGTATTCAGAAc                                                                          >  4:447798/1‑274 (MQ=255)
                                                                                                                              aCTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGtt                                                                                                                                                                                      >  2:255245/1‑164 (MQ=255)
                                                                                                                              aCTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGccc                                                                                                                                                              >  3:43029/1‑188 (MQ=255)
                                                                                                                              aCTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCctgctg                                                                                                                            >  2:247605/1‑222 (MQ=255)
                                                                                                                                  aGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAAt                                                                                                                                                                                  <  4:255244/164‑1 (MQ=255)
                                                                                                                                  aGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTa                                                                                                                                                          <  1:43030/188‑1 (MQ=255)
                                                                                                                                  aGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGcc                                                                                                                        <  4:247604/222‑1 (MQ=255)
                                                                                                                                     aaTAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACa                                                                                                                                                                                                                             <  8:514805/118‑1 (MQ=255)
                                                                                                                                     aaTAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACa                                                                                                                                                                                                                             >  7:514805/1‑118 (MQ=255)
                                                                                                                                          aaCACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGcc                                                                                                                                                                                          >  5:57885/1‑148 (MQ=255)
                                                                                                                                          aaCACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGc                                                                                                                                                                >  2:586070/1‑174 (MQ=255)
                                                                                                                                              cTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGtt                                                                                                                                                                                      >  3:169328/1‑148 (MQ=255)
                                                                                                                                               tCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCAGGAGGTATTGTTCAGGTATTCAGAACCGCGTAACTTTTGCCCatat                                                      >  2:599688/1‑275 (MQ=255)
                                                                                                                                                 acaACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAAt                                                                                                                                                                                  >  7:849234/1‑149 (MQ=255)
                                                                                                                                                  caACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAAt                                                                                                                                                                                  <  1:169329/148‑1 (MQ=255)
                                                                                                                                                      ccAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGc                                                                                                                                                                             >  7:499553/1‑149 (MQ=255)
                                                                                                                                                       cAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCg                                                                                                                                                                            >  7:45562/1‑149 (MQ=255)
                                                                                                                                                         aaaaTAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGaa                                                                                                                                                                          <  8:45562/149‑1 (MQ=255)
                                                                                                                                                         aaaaTAGAGCATCCGGAAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGa                                                                                                                                                                           >  5:346928/1‑148 (MQ=255)
                                                                                                                                                                agCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTg                                                                                                                                                                                                                                    >  2:574537/1‑84 (MQ=255)
                                                                                                                                                                agCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCt                                                                                                                                                                                                                                  >  2:562084/1‑86 (MQ=255)
                                                                                                                                                                 gCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTg                                                                                                                                                                                                                                    >  3:32699/1‑83 (MQ=255)
                                                                                                                                                                    tCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTaa                                                                                                                                                                                                                                <  4:574535/84‑1 (MQ=255)
                                                                                                                                                                    tCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAAc                                                                                                                                                                                                                              <  4:562082/86‑1 (MQ=255)
                                                                                                                                                                     ccGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTaa                                                                                                                                                                                                                                <  1:32700/83‑1 (MQ=255)
                                                                                                                                                                     ccGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGccc                                                                                                                                                              <  8:499553/149‑1 (MQ=255)
                                                                                                                                                                        gCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCAGGAGGTATTGTTCAGGTATTCAGAACCGCGTAACTTTTGCCCATATGCAGCTACTGCTGCGCGAAAACGCa                             <  2:447800/275‑1 (MQ=255)
                                                                                                                                                                           aaaGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTaa                                                                                                                                                         >  5:811428/1‑148 (MQ=255)
                                                                                                                                                                                gtgagtTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCtgt                                                                                                                                                   <  8:849234/149‑1 (MQ=255)
                                                                                                                                                                                 tgagtTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGcc                                                                                                                                                               >  4:806523/1‑136 (MQ=255)
                                                                                                                                                                                   agtTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCAGGAGGTATTGTTCAGGTATTCAGAACCGCGTAACTTTTGCCCATATGCAGCTACTGCTGCGCGAAAACGCATCTCTTCagag                  >  2:519169/1‑275 (MQ=255)
                                                                                                                                                                                     ttAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGt                                                                                                                                                           <  2:806525/136‑1 (MQ=255)
                                                                                                                                                                                         catcatTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCAGGAGGTATTGTTCAGGTATTCAGAACCGCGTAACTTTTGCCCATATGCAGCTACTGCTGCGCGAAAACGCATATCTTCAGAGGgcgc             >  4:620061/1‑274 (MQ=255)
                                                                                                                                                                                              ttGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCtt                                                                                                                                      >  5:778797/1‑148 (MQ=255)
                                                                                                                                                                                               tGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGACTAAATTCCACCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCGGCTGCCAGGAGGTATTGTTCAGGTATTCAGAACCGTGTAACTTTTGCCCATATGCAGCTACTGTTGCGCGAAATAGCATTTCTACAGAGGGCGTGTCTAtt      >  4:257024/1‑275 (MQ=255)
                                                                                                                                                                                               tGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAAcc                                                                                                                                                       >  3:83300/1‑130 (MQ=255)
                                                                                                                                                                                                   aattaattCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAAt                                                                                                                                                                                  >  3:42025/1‑99 (MQ=255)
                                                                                                                                                                                                   aattaattCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCtgt                                                                                                                                                   <  1:83301/130‑1 (MQ=255)
                                                                                                                                                                                                    attaattCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCAGGAGGTATTGTTCAGGTATTCAGAACCGCGTAACTTTTGCCCATATGCAGCTACTGCTGCGCGAAAACGCATATCTTCAGAGGGCGCGTTTAttgttg  >  2:5389/1‑274 (MQ=255)
                                                                                                                                                                                                       aattCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAgg                                                                                                                                                                              <  1:42026/99‑1 (MQ=255)
                                                                                                                                                                                                        attCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTc                                                                                                                                                                   >  4:745598/1‑109 (MQ=255)
                                                                                                                                                                                                        attCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGa                                                                                                                           >  7:423695/1‑149 (MQ=255)
                                                                                                                                                                                                            gACCAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCtgtg                                                                                                                                                  <  8:499870/122‑1 (MQ=255)
                                                                                                                                                                                                            gACCAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCtgtg                                                                                                                                                  >  7:499870/1‑122 (MQ=255)
                                                                                                                                                                                                            gACCAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTc                                                                                                                                           >  8:517866/1‑129 (MQ=255)
                                                                                                                                                                                                            gACCAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTc                                                                                                                                           <  7:517866/129‑1 (MQ=255)
                                                                                                                                                                                                            gACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGcc                                                                                                                                                               <  2:745600/109‑1 (MQ=255)
                                                                                                                                                                                                             aCCAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCtt                                                                                                                                            >  6:729938/1‑127 (MQ=255)
                                                                                                                                                                                                             aCCAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCtt                                                                                                                                            <  5:729938/127‑1 (MQ=255)
                                                                                                                                                                                                             aCCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCactgac                                                                                                                                                                                                                 >  3:48329/1‑58 (MQ=255)
                                                                                                                                                                                                              ccAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCttt                                                                                                                                     <  6:140758/133‑1 (MQ=255)
                                                                                                                                                                                                              ccAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCttt                                                                                                                                     >  5:140758/1‑133 (MQ=255)
                                                                                                                                                                                                              ccAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCAc                                                                                                                     >  8:626633/1‑149 (MQ=255)
                                                                                                                                                                                                               cAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGa                                                                                                                                                <  8:840213/121‑1 (MQ=255)
                                                                                                                                                                                                               cAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGa                                                                                                                                                >  7:840213/1‑121 (MQ=255)
                                                                                                                                                                                                                aGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACt                                                                                                                    >  5:355240/1‑148 (MQ=255)
                                                                                                                                                                                                                aGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACt                                                                                                                    >  6:702869/1‑148 (MQ=255)
                                                                                                                                                                                                                 gTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGAt                                                                                                                                                                                                             <  1:48330/58‑1 (MQ=255)
                                                                                                                                                                                                                 gTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTGGCCTGCTGAGCCCACTGTTTGNTCtgntgc                                                                                                      >  4:745476/1‑161 (MQ=255)
                                                                                                                                                                                                                  tCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGtt                                                                                                                 >  8:2490/1‑149 (MQ=255)
                                                                                                                                                                                                                  tCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGt                                                                                                                  >  5:785852/1‑148 (MQ=255)
                                                                                                                                                                                                                  tCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTcccc                                                                                                                                        <  8:47077/126‑1 (MQ=255)
                                                                                                                                                                                                                  tCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTcccc                                                                                                                                        >  7:47077/1‑126 (MQ=255)
                                                                                                                                                                                                                  tCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGt                                                                                                                  <  6:778797/148‑1 (MQ=255)
                                                                                                                                                                                                                    aGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTg                                                                                                                                                            >  8:336958/1‑104 (MQ=255)
                                                                                                                                                                                                                    aGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTg                                                                                                                                                            <  7:336958/104‑1 (MQ=255)
                                                                                                                                                                                                                     gCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGtttta                                                                                                              >  8:25690/1‑147 (MQ=255)
                                                                                                                                                                                                                     gCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCagg                                                                                                  <  2:745478/161‑1 (MQ=255)
                                                                                                                                                                                                                      cAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGa                                                                                                              >  6:876673/1‑148 (MQ=255)
                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                          
ATGCTGACAGATTTAAAAAATATTCAGCAGAAGAGCAGTACGAAATTGATTTTCAACATCGACTGGCAAACAGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCAGGAGGTATTGTTCAGGTATTCAGAACCGCGTAACTTTTGCCCATATGCAGCTACTGCTGCGCGAAAACGCATATCTTCAGAGGGCGCGTTTATTGTTG  >  NC_000913/2383740‑2384209

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