Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,383,959 A→T 16.3% 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.3% 181.4 / 25.0 66intergenic (+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 (30/25);  new base T (5/6);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCTGACAGCGCCCCCTATACGTTATTAAAAAACCATGCTGACAGATTTAAAAAATATTCAGCAGAAGAGCAGTACGAAATTGATTTTCAACATCGACTGGCAAACAGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCAGGAGGTATTGTTCAGGTATTCAGAACCGCGTAACTTTTGCCCATATGCAGCTACTGCTGCGCGAAAACGCATATCTTCAGAGGGCGCGTTTATTGTTGGCCCCAGCGGGAATT  >  NC_000913/2383704‑2384224
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                         
aTTCTGACAGCGCCCCCTATACGTTATTAAAAAACCATGCTGACAGATTTAAAAAATATTCAGCAGAAGAGCAGTACGAAATTGATTTTCAACATCGACTGGCAAACAGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATa                                                                                                                                                                                                                                                        <  2:250046/275‑1 (MQ=255)
aTTCTAACAGCGCCCCCTATACGTTATTAAAAAACCATGCTGACAGATTTAAAAAATATTCAGCAGAAGAGCAGTACGAAATTGATTTTCAACATCGACTGGCAAACAGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTAAAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATa                                                                                                                                                                                                                                                        <  1:53344/275‑1 (MQ=255)
                         attaAAAAACCATGCTGACAGATTTAAAAAATATTCAGCAGAAGAGCAGTACGAAATTGATTTTCAACATCGACTGGCAAACAGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTCGCAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCCAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAAAAGTCATCactgac                                                                                                                                                                                                                                >  2:107077/1‑274 (MQ=255)
                         attaAAAAACCATGCTGACAGATTTAAAAAATATTCAGCAGAAGAGCAGTACGAAATTGATTTTCAACATCGACTGGCAAACAGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGAAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCactgac                                                                                                                                                                                                                                >  2:28706/1‑274 (MQ=255)
                         attaAAAAACCATGCTGACAGATTTAAAAAATATTCAGCAGAAGAGCAGTACGAAATTGATTTTCAACATCGACTGGCAAACAGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCAGTGTTAATTAATTCGACCAGTCAGCAAAAATGCTAACTGGTCGAACATAGTTTCCTAACAAGTAAtcactcac                                                                                                                                                                                                                                >  4:324540/1‑274 (MQ=255)
                                                                                                                                            aaaatcaatcaTGATGTTAGAAAATTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGGGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCAGGAGGTa                                                                                                            <  4:225937/272‑1 (MQ=255)
                                                                                                                                               atcaatcaTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCAGGAGGTATTg                                                                                                         >  2:283008/1‑275 (MQ=255)
                                                                                                                                                   atcaTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGa                                                                                                                                                                                                                                                                     >  8:1043829/1‑115 (MQ=255)
                                                                                                                                                   atcaTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGa                                                                                                                                                                                                                                                                     <  7:1043829/115‑1 (MQ=255)
                                                                                                                                                    tcaTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAatat                                                                                                                                                                                                                                                                          <  8:1093339/109‑1 (MQ=255)
                                                                                                                                                    tcaTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAatat                                                                                                                                                                                                                                                                          >  7:1093339/1‑109 (MQ=255)
                                                                                                                                                    tcaTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCAGGAGGTATTGTTCAg                                                                                                    <  2:486462/275‑1 (MQ=255)
                                                                                                                                                                  aCTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGcc                                                                                                                                       >  2:110115/1‑226 (MQ=255)
                                                                                                                                                                   cTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACt                                                                                                                                                                                                                                                                   >  8:927522/1‑101 (MQ=255)
                                                                                                                                                                   cTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACt                                                                                                                                                                                                                                                                   <  7:927522/101‑1 (MQ=255)
                                                                                                                                                                   attAGAAATAAAACACTCACAACCCAAAAAAAGACCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGAAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCAGGAGGTATTGTTCAGGTATTCAGAACcgcg                                                                                     <  2:324541/274‑1 (MQ=255)
                                                                                                                                                                      aGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACt                                                                                                                                   <  4:110114/226‑1 (MQ=255)
                                                                                                                                                                            aaaaCACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGc                                                                                                                                                                                                                                                                        <  8:633580/87‑1 (MQ=255)
                                                                                                                                                                            aaaaCACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGc                                                                                                                                                                                                                                                                        >  7:633580/1‑87 (MQ=255)
                                                                                                                                                                              aaCACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCAGGAGGTATTGTTCAGGTATTCAGAACCGCGTAACTTTTGcc                                                                          <  4:283007/275‑1 (MQ=255)
                                                                                                                                                                                    cacaACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCt                                                                                                                                                                                                                                                                       >  7:388979/1‑80 (MQ=255)
                                                                                                                                                                                    cacaACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCt                                                                                                                                                                                                                                                                       <  8:388979/80‑1 (MQ=255)
                                                                                                                                                                                           cAAAAATAGAGCATCCGGCAAAAGGGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAGATTTg                                                                                                                                                                                                                                                                         >  6:390108/1‑71 (MQ=255)
                                                                                                                                                                                           cAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTg                                                                                                                                                                                                                                                                         <  5:390108/71‑1 (MQ=255)
                                                                                                                                                                                             aaaaTAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGACt                                                                                                                                                                                                                                                                   >  7:1137920/1‑75 (MQ=255)
                                                                                                                                                                                             aaaaTAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGACt                                                                                                                                                                                                                                                                   <  8:1137920/75‑1 (MQ=255)
                                                                                                                                                                                                   gagCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCAGGAGGTATTGTTCAGGTATTCAGAACCGCGTAACTTTTGCCCATATGCAGCTACTGCTgcgc                                                     >  4:41906/1‑275 (MQ=255)
                                                                                                                                                                                                    agCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATCAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCAGGAGGTATTGTTCAGGTATTCAGAACCGCGTAACTTTTGCCCATATGCAGCTACTGCTGcgcg                                                    >  4:129483/1‑275 (MQ=255)
                                                                                                                                                                                                    agCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGa                                                                                                                                                               >  4:115100/1‑168 (MQ=255)
                                                                                                                                                                                                    agCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGNTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTTATCCCACTGTNTGACCTGCTTCCAGGCGTTATTGTTAGGTTTTTCAGACCCGCGTANANTTTTCCCATATGCattcacttttgccg                                                     >  3:50743/1‑261 (MQ=255)
                                                                                                                                                                                                    agCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACAATCCGACAGTTCAATCAGGCGAATAACTTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCCGGAGGTATTGTTCAGGTATTCAGCACCGCGTAACTTTTCCCCATATGCCGCTATTGCTccgc                                                     >  4:153773/1‑274 (MQ=255)
                                                                                                                                                                                                     gCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCAGGCGGTATTGTTCAGGTATTCAGAACCGCGTAACTTTTGCCCATATGCAGCTACTGCTGCGCGa                                                   >  4:366956/1‑275 (MQ=255)
                                                                                                                                                                                                      cATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCAGGAGGTATTGTTCAGGTATTCAGAACCGCGTAACTTTTGCCCATATGCAGCTACTGCTGCGCGaa                                                  >  1:67178/1‑275 (MQ=255)
                                                                                                                                                                                                       aTCCGGCAACAGAGTGAGTTAACATAATTGTTAATTAATTCGACCAGGCAGCAAATATGCTGAc                                                                                                                                                                                                                                                                    >  6:661008/1‑64 (MQ=255)
                                                                                                                                                                                                       aTCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGAc                                                                                                                                                                                                                                                                    <  5:661008/64‑1 (MQ=255)
                                                                                                                                                                                                        tCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCtt                                                                                                                                                           <  2:115101/168‑1 (MQ=255)
                                                                                                                                                                                                            gaaaaaagattagtTAACATCTATGTTAATTAATTCGACCAGTCAGCAAATTATCTGACTGGTCGACAATATTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGTCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCAGGAGGTATTGTTCAGGTATTCAGAACCGCGTAACTTTTGCCCATATGCAGCTACTGCTGCGCGAAAACGCa                                            <  4:28705/264‑1 (MQ=255)
                                                                                                                                                                                                                aaGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCt                                                                                                                                                                            >  2:439163/1‑143 (MQ=255)
                                                                                                                                                                                                                    gtgagtTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTaa                                                                                                                                                                        <  4:439162/143‑1 (MQ=255)
                                                                                                                                                                                                                         gtAAAATCAATGTTAAATAAATCGAACAGTCAGCAAATATGCTGAATGGTCGAAAATAGTTTGCTAAAAAGTAATCAATTACTGATGTCAGTACCACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTGGCCTGCTGAGCCCACTGTTTGATCTGCTGCCAGGAGGTATTGTTCAGGTATTCAGAACCGCGTAACTTTGGCCCATATGCAGCTACTGCTGCGCGAAAACGCATATCTTCAGAggg                               <  4:214773/274‑1 (MQ=255)
                                                                                                                                                                                                                                   tGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAgg                                                                                                                                                                                             >  3:154511/1‑107 (MQ=255)
                                                                                                                                                                                                                                   tGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAAt                                                                                                                                                                                    >  2:544509/1‑116 (MQ=255)
                                                                                                                                                                                                                                       aattaattCGACCAGTCAGCAAATATGCTGACTGTGCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGaa                                                                                                                                                                                         <  1:154511/107‑1 (MQ=255)
                                                                                                                                                                                                                                       aattaattCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCg                                                                                                                                                                                <  4:544508/116‑1 (MQ=255)
                                                                                                                                                                                                                                                    aGTCAGTAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCttt                                                                                                                                                    <  5:409007/131‑1 (MQ=255)
                                                                                                                                                                                                                                                    aGTCAGTAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCttt                                                                                                                                                    >  6:409007/1‑131 (MQ=255)
                                                                                                                                                                                                                                                    aGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTcac                                                                                                                                                                                                                   <  5:850618/68‑1 (MQ=255)
                                                                                                                                                                                                                                                    aGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTcac                                                                                                                                                                                                                   >  6:850618/1‑68 (MQ=255)
                                                                                                                                                                                                                                                      tCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTcc                                                                                                                                                         <  6:1080017/124‑1 (MQ=255)
                                                                                                                                                                                                                                                      tCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTcc                                                                                                                                                         >  5:1080017/1‑124 (MQ=255)
                                                                                                                                                                                                                                                      tCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCAGGAGGTATTGTTCAGGTATTCAGAACCGCGTAACTTTTGCCCATATGCAGCTACTGCTGCGCGAAAACGCATATCTTCAGAGGGCGCGTTTATTGTTGGCCCCAGCGGGAAtt  >  2:464491/1‑275 (MQ=255)
                                                                                                                                                                                                                                                       cAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAgg                                                                                                                                                                                             >  5:1089170/1‑87 (MQ=255)
                                                                                                                                                                                                                                                       cAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAgg                                                                                                                                                                                             <  6:1089170/87‑1 (MQ=255)
                                                                                                                                                                                                                                                       cAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGc                                                                                                                                                                               <  8:648311/101‑1 (MQ=255)
                                                                                                                                                                                                                                                       cAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGc                                                                                                                                                                               >  7:648311/1‑101 (MQ=255)
                                                                                                                                                                                                                                                       cAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTg                                                                                                                                                                           >  8:764026/1‑105 (MQ=255)
                                                                                                                                                                                                                                                       cAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTg                                                                                                                                                                           <  7:764026/105‑1 (MQ=255)
                                                                                                                                                                                                                                                         gCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGAt                                                                                                                                                              <  5:1094281/116‑1 (MQ=255)
                                                                                                                                                                                                                                                         gCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGAt                                                                                                                                                              >  6:1094281/1‑116 (MQ=255)
                                                                                                                                                                                                                                                         gCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCAACACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCtgt                                                                                                                                                                  >  8:465555/1‑112 (MQ=255)
                                                                                                                                                                                                                                                         gCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCAACACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCtgt                                                                                                                                                                  <  7:465555/112‑1 (MQ=255)
                                                                                                                                                                                                                                                          cAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGt                                                                                                                                                                                                      >  7:854523/1‑75 (MQ=255)
                                                                                                                                                                                                                                                          cAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGt                                                                                                                                                                                                      <  6:1060350/75‑1 (MQ=255)
                                                                                                                                                                                                                                                          cAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGt                                                                                                                                                                                                      >  5:1060350/1‑75 (MQ=255)
                                                                                                                                                                                                                                                          cAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGt                                                                                                                                                                                                      <  8:854523/75‑1 (MQ=255)
                                                                                                                                                                                                                                                          cAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCtt                                                                                                                                                           <  8:344881/118‑1 (MQ=255)
                                                                                                                                                                                                                                                          cAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCtt                                                                                                                                                           >  7:344881/1‑118 (MQ=255)
                                                                                                                                                                                                                                                          cAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGa                                                                                                                             >  8:295756/1‑148 (MQ=255)
                                                                                                                                                                                                                                                               atGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCttt                                                                                                                                                    <  8:610538/120‑1 (MQ=255)
                                                                                                                                                                                                                                                               atGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCtgt                                                                                                                                                    >  7:610538/1‑118 (MQ=255)
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                         
ATTCTGACAGCGCCCCCTATACGTTATTAAAAAACCATGCTGACAGATTTAAAAAATATTCAGCAGAAGAGCAGTACGAAATTGATTTTCAACATCGACTGGCAAACAGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTGCCAGGAGGTATTGTTCAGGTATTCAGAACCGCGTAACTTTTGCCCATATGCAGCTACTGCTGCGCGAAAACGCATATCTTCAGAGGGCGCGTTTATTGTTGGCCCCAGCGGGAATT  >  NC_000913/2383704‑2384224

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