Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,507,501 A→C 5.3% I241L (ATT→CTT)  ydcP → putative peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,507,5010AC 5.3% 536.6 / 19.0 156I241L (ATT→CTT) ydcPputative peptidase
Reads supporting (aligned to +/- strand):  ref base A (75/71);  new base C (5/4);  total (81/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.21e-01

Only 100 of 160 total aligned reads displayed.

CGCGATTCACCAGGCTACGGACGCGACCATTGAATTCTTTATTCATGGGGCACTGTGCGTGGCCTATTCGGGTCAGTGCTACATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGCACAGATTATTTTGTGAATGCCCGTAAAGGCGATATTGGCGCGTTCGATTCGCCGAAATTTATCGGCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATC  >  NC_000913/1507230‑1507771
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                              
cgcgATTCACCAGGCTACGGACGCGACCATTGAATTCTTTATTCATGGGGCACTGTGCGTGGCCTATTCGGGTCAGTGCTACATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGGTGGTGTACGCTCCTTCAAGAtgt                                                                                                                                                                                                                                                                             >  4:109064/1‑273 (MQ=255)
            gtcgcagggcgcgACCCATGATGTCTTTATTAAAGGGGCACTGTGCGGGGACTAGCCGGTTCAGTGCTACATCTTTAATGCGCAAACAGGGCGTAGAGCCAACCGTGGCGATGGCTCGCAGGCGTGCCGTGTGCCAGACACATTGAAAGACGAGCAGGTGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGAGAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACa                                                                                                                                                                                                                                                                 <  4:15860/266‑1 (MQ=255)
            ggCTACGGACGCGACCATTGAATTCTTTATTCATGGGGCACTGTGCGTGGCCTATTCGGGTCAGTGCTACATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACa                                                                                                                                                                                                                                                                 <  4:50941/275‑1 (MQ=255)
               tACGGACGCGACCATTGAATTCTTTATTCATGGGGCACTGTGCGTGGCCTATTCGGGTCAGTGCTACATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGCCTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTCCGCTCCTTCAAGATTGAAGGGCGt                                                                                                                                                                                                                                                                     >  1:140812/1‑268 (MQ=255)
                 cGGACGCGACCATTGAATTCTTTATTCATGGGGCACTGTGCGTGGCCTATTCGGGTCAGTGCTACATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGat                                                                                                                                                                                                                                                            >  2:339302/1‑275 (MQ=255)
                      gcaaCCATTGAATACTTTATTCATGGGGAACTGTGCGTGGCCGATTCGGGTCATCGGTACGTTTATCATGCGCAAACAGGACGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCATTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAGAAACATCTGCTGTCGATGAAGGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAg                                                                                                                                                                                                                                                       <  4:259982/272‑1 (MQ=255)
                            aTTTAATTCTTTTTCATGGGGGACTGTGGGTGGCCTATTCGGGTCAGTGCTACATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGGGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTg                                                                                                                                                                                                                                                <  4:20150/275‑1 (MQ=255)
                              tGAATTCTTGCTTCATGGGGCACTGTGCGTGGCCGAGTCGGGTCAGTGCTAAATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGGGGCGATTGCTCGCAGGCGTGCCGGTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTg                                                                                                                                                                                                                                                <  4:244337/274‑1 (MQ=255)
                                                 ggactggggctGGCCGATTCCGGGCCGTGCTAAAGTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTGGGAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGTGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGTTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGCTATGAGCTACGTGAAGAATATCACCGTCCATTa                                                                                                                                                                                                                            <  4:192803/265‑1 (MQ=255)
                                                       tgCGTGGCCTATTCGGGTCAGTGCTACATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCCCCGCCCATTATCGCCa                                                                                                                                                                                                                      >  2:219477/1‑275 (MQ=255)
                                                                   tCGGGTCAGTGCTACATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGACTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGc                                                                                                                                                                                                          <  2:109065/275‑1 (MQ=255)
                                                                      ggTCAGTGCTACATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCa                                                                                                                                                                                                        <  2:148344/274‑1 (MQ=255)
                                                                            tGCTACATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTAACGGGCGTTACAAAGATATGAGCTACGTGCAGAATAACACCGCCCATTATCGCCAGATGCTTGCTGCCATTTTTg                                                                                                                                                                                                  >  4:73215/1‑274 (MQ=255)
                                                                             gCTGCATTTCGCATGCGCAAACAGGGCTTAGCGCCAACCGGGGCGAGTGCTCACAGGCGTGCCCTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTgaa                                                                                                                                                                                                <  4:482215/275‑1 (MQ=255)
                                                                                 cATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTCTGAAAAACATCTCCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCCCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACCAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGTTGCTTTATGCCATTTTTgaaaaa                                                                                                                                                                                             >  4:417405/1‑274 (MQ=255)
                                                                                      ctcATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCCGATTCTTGATGCCATTATTGAAAAACGTgcg                                                                                                                                                                                       >  4:52324/1‑273 (MQ=255)
                                                                                            cgcAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATAGCCCGATGCTTGATGCCATTATTGAAGAAAGTGGCGAACTg                                                                                                                                                                                 >  4:379768/1‑275 (MQ=255)
                                                                                                           gcgcCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGata                                                                                                                                                                                                                                                           >  1:79136/1‑186 (MQ=255)
                                                                                                                           ttGCTCGCAGGCATGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACa                                                                                                                                                    <  2:154472/273‑1 (MQ=255)
                                                                                                                                cGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGa                                                                                                                                                                                                                                                                            <  6:499189/148‑1 (MQ=255)
                                                                                                                                      cGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCg                                                                                                                                                                                                                                                                      <  5:872194/148‑1 (MQ=255)
                                                                                                                                        tGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGAACTGGGCGGCCCTTCATTGGGTCGTACTGAACCTTTCtctttatccc                                                                                                                                      >  4:74725/1‑265 (MQ=255)
                                                                                                                                           cGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGaa                                                                                                                                                                                                                                                                           >  5:1100057/1‑138 (MQ=255)
                                                                                                                                            gTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAgg                                                                                                                                                                <  2:74726/244‑1 (MQ=255)
                                                                                                                                              ttGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAg                                                                                                                                                                                                                                                              <  5:781580/148‑1 (MQ=255)
                                                                                                                                                      cacaTTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGa                                                                                                                                                                                                                                                                            >  8:391055/1‑126 (MQ=255)
                                                                                                                                                      cacaTTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGa                                                                                                                                                                                                                                                                            <  7:391055/126‑1 (MQ=255)
                                                                                                                                                          ttGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAgg                                                                                                                                                                                                                                                                         <  5:904823/125‑1 (MQ=255)
                                                                                                                                                          ttGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAgg                                                                                                                                                                                                                                                                         >  6:904823/1‑125 (MQ=255)
                                                                                                                                                            gAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTAAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGAAAAGACTTTCCAc                                                                                                                  <  2:73216/274‑1 (MQ=255)
                                                                                                                                                             aaaGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCCCTTCAATCTTGAAggagcgta                                                                                                                                                                                                                                                                   >  6:896859/1‑121 (MQ=255)
                                                                                                                                                              aaGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAg                                                                                                                                                                                                                                                       >  6:183985/1‑139 (MQ=255)
                                                                                                                                                              aaGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAg                                                                                                                                                                                                                                                       <  5:183985/139‑1 (MQ=255)
                                                                                                                                                                   gATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGata                                                                                                                                                                                                                                                           >  3:30068/1‑130 (MQ=255)
                                                                                                                                                                      cAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATTACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGTTGTACGCTCCTTCAAGATTGAAGGGCGTTAAAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGcac                                                                                                       <  1:70086/275‑1 (MQ=255)
                                                                                                                                                                       aGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAg                                                                                                                                                                                                                                                       <  1:30068/130‑1 (MQ=255)
                                                                                                                                                                         gggCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCCCTTCAATCTTGAAggagc                                                                                                                                                                                                                                                                      >  8:788636/1‑109 (MQ=255)
                                                                                                                                                                         gggCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCCCTTCAATCTTGAAggagc                                                                                                                                                                                                                                                                      <  7:788636/112‑4 (MQ=255)
                                                                                                                                                                               gtggtTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCa                                                                                                                                                                                                                               <  2:272900/146‑1 (MQ=255)
                                                                                                                                                                                     tCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACaa                                                                                                                                                                                                                                                                <  7:545203/107‑1 (MQ=255)
                                                                                                                                                                                     tCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACaa                                                                                                                                                                                                                                                                >  8:545203/1‑107 (MQ=255)
                                                                                                                                                                                      ccTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGat                                                                                                                                                                                                                                                            >  6:979422/1‑110 (MQ=255)
                                                                                                                                                                                      ccTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGat                                                                                                                                                                                                                                                            <  5:979422/110‑1 (MQ=255)
                                                                                                                                                                                          tGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaa                                                                                                                                                                                                                                              >  6:165876/1‑120 (MQ=255)
                                                                                                                                                                                          tGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaa                                                                                                                                                                                                                                              <  5:165876/120‑1 (MQ=255)
                                                                                                                                                                                             aaaaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGc                                                                                                                                                                                                                        <  7:573483/139‑1 (MQ=255)
                                                                                                                                                                                             aaaaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGc                                                                                                                                                                                                                        >  8:573483/1‑139 (MQ=255)
                                                                                                                                                                                               aaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGa                                                                                                                                                                                                             <  7:161892/148‑1 (MQ=255)
                                                                                                                                                                                               aaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGa                                                                                                                                                                                                             >  8:1012206/1‑148 (MQ=255)
                                                                                                                                                                                                   tCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGa                                                                                                                                                                                                                                                        <  7:885824/101‑1 (MQ=255)
                                                                                                                                                                                                   tCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGa                                                                                                                                                                                                                                                        >  8:885824/1‑101 (MQ=255)
                                                                                                                                                                                                      gctgTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGt                                                                                                                                                                                                                                                                     <  7:11041/85‑1 (MQ=255)
                                                                                                                                                                                                      gctgTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGt                                                                                                                                                                                                                                                                     >  8:11041/1‑85 (MQ=255)
                                                                                                                                                                                                          tCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGCACAGATTATTTTGTGAATGCCCGTAAAGGCGATATTgg                                                                   <  2:263735/275‑1 (MQ=255)
                                                                                                                                                                                                             aTGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGa                                                                                                                                                                                                                                                                            <  5:592837/71‑1 (MQ=255)
                                                                                                                                                                                                             aTGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGa                                                                                                                                                                                                                                                                            >  6:592837/1‑71 (MQ=255)
                                                                                                                                                                                                             aTGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGt                                                                                                                                                                                                                                                                     <  5:401218/78‑1 (MQ=255)
                                                                                                                                                                                                             aTGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGt                                                                                                                                                                                                                                                                     >  6:401218/1‑78 (MQ=255)
                                                                                                                                                                                                               gAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAg                                                                                                                                                                                                                                                       >  7:929447/1‑90 (MQ=255)
                                                                                                                                                                                                               gAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAg                                                                                                                                                                                                                                                       <  8:929447/90‑1 (MQ=255)
                                                                                                                                                                                                               gAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTgaagaa                                                                                                                                                                                             <  7:1012206/148‑1 (MQ=255)
                                                                                                                                                                                                                aaaGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGc                                                                                                                                                                                                                                  <  7:6314/110‑1 (MQ=255)
                                                                                                                                                                                                                aaaGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGc                                                                                                                                                                                                                                  >  8:6314/1‑110 (MQ=255)
                                                                                                                                                                                                                     tAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTgg                                                                                                                                                                                        <  5:473993/147‑1 (MQ=255)
                                                                                                                                                                                                                       aCGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaa                                                                                                                                                                                                                                              <  7:877607/91‑1 (MQ=255)
                                                                                                                                                                                                                       aCGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaa                                                                                                                                                                                                                                              >  8:877607/1‑91 (MQ=255)
                                                                                                                                                                                                                          aTCAGACTGCCAACCTCGGCTCGCTGATTCATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGTTGCCATTATTGAATAACGTGGCGATCTTCCGCGCCCTTCATCAGGTCGTACTGTACATTTCTTTGTTCCATCGCCGCAAAAGCCTTTCCACCGTGTGAGCTAAGCTTATGTTGTTAATGCCCGTAATGGCGATGTTGTCGTTGTCGCTTCGCCg                                                   >  4:209452/1‑275 (MQ=255)
                                                                                                                                                                                                                          aTCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTg                                                                                                                                                                                                              >  6:1083617/1‑120 (MQ=255)
                                                                                                                                                                                                                          aTCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCattat                                                                                                                                                                                                    >  6:1106294/1‑130 (MQ=255)
                                                                                                                                                                                                                          aTCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCattat                                                                                                                                                                                                    <  5:1106294/130‑1 (MQ=255)
                                                                                                                                                                                                                            cAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGc                                                                                                                                                                                                                        <  8:39311/108‑1 (MQ=255)
                                                                                                                                                                                                                            cAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGc                                                                                                                                                                                                                        >  7:39311/1‑108 (MQ=255)
                                                                                                                                                                                                                            cAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCatt                                                                                                                                                                                                      <  8:830042/126‑1 (MQ=255)
                                                                                                                                                                                                                            cAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAAGATCACCGCCCATTATCGCCAGATGCTTGATGCCatt                                                                                                                                                                                                      >  7:830042/1‑126 (MQ=255)
                                                                                                                                                                                                                             aGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCattat                                                                                                                                                                                                    >  8:19799/1‑127 (MQ=255)
                                                                                                                                                                                                                             aGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCattat                                                                                                                                                                                                    <  7:19799/127‑1 (MQ=255)
                                                                                                                                                                                                                                  gCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTg                                                                                                                                                                                         <  6:772256/133‑1 (MQ=255)
                                                                                                                                                                                                                                   ccAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCNGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGCACAGATTNAGTTGTGAATGCCCGTAAAGGCGAANTTGGGGCGTTCGATTCGCCGAAATTTTTc                                          >  3:91207/1‑275 (MQ=255)
                                                                                                                                                                                                                                    cAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCt                                                                                                                                                                                                                >  5:825248/1‑108 (MQ=255)
                                                                                                                                                                                                                                     aaCCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAg                                                                                                                                                                                                                                                                          <  6:104839/49‑1 (MQ=255)
                                                                                                                                                                                                                                     aaCCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAg                                                                                                                                                                                                                                                                          >  5:104839/1‑49 (MQ=255)
                                                                                                                                                                                                                                      aCCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTtc                                                                                                                                                                      >  8:281245/1‑148 (MQ=255)
                                                                                                                                                                                                                                       ccTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGAGTGAAGGGCGTTAAAAAGATATGAGCTACGTGAAGAATATTACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGCACAGATTATTTTGTGAATGCCCGTAAAGGCGATATTGGCGCGTTCGATTCGCCGAAATTTATCGGcc                                      <  2:379769/275‑1 (MQ=255)
                                                                                                                                                                                                                                         tCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTtcatc                                                                                                                                                                   <  7:281245/148‑1 (MQ=255)
                                                                                                                                                                                                                                                    attgatGCTGGTGGACGCACCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGCACAGATTATTTTGTGAATGCCCGTAAAGGCGATATTGGCGCGTTCGATTCGCCGAAATTTATCGGCCTGCCGGTAGGCGa                         <  2:247666/275‑1 (MQ=255)
                                                                                                                                                                                                                                                          gCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTg                                                                                                                                                                                         >  6:390934/1‑109 (MQ=255)
                                                                                                                                                                                                                                                           cTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTga                                                                                                                                                                                                 <  6:923970/100‑1 (MQ=255)
                                                                                                                                                                                                                                                             ggTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTg                                                                                                                                                        <  5:747669/139‑1 (MQ=255)
                                                                                                                                                                                                                                                             ggTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTg                                                                                                                                                        >  6:747669/1‑139 (MQ=255)
                                                                                                                                                                                                                                                             ggTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCt                                                                                                                                               >  5:36836/1‑148 (MQ=255)
                                                                                                                                                                                                                                                              gtgtACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGCACAGATTATTTTGTGAATGCCCGTAAAGGCGATATTGGCGCGTTCGATTCGCCGAAATTTATCGGCTTGCCGGTAGGCGAAGTAGTGaa                >  2:532505/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                    gCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTgaagaa                                                                                                                                                                                             >  6:823452/1‑95 (MQ=255)
                                                                                                                                                                                                                                                                      tCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCattat                                                                                                                                                                                                    <  6:662064/86‑1 (MQ=255)
                                                                                                                                                                                                                                                                      tCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCattat                                                                                                                                                                                                    >  5:662064/1‑86 (MQ=255)
                                                                                                                                                                                                                                                                        cTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGCACAGATTATTTTGTGAATGCCCGTAAAGGCGATATTGGCGCGTTCGATTCGCCGAAATTTATCGGCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAg     >  2:201190/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                        cTTCAAGATTGAAGGGAGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCCTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGca                                                                                                        >  4:248064/1‑176 (MQ=255)
                                                                                                                                                                                                                                                                          tCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCt                                                                                                                                               <  5:816476/135‑1 (MQ=255)
                                                                                                                                                                                                                                                                          tCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGAc                                                                                                                                  >  5:96661/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                            aaGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCatt                                                                                                                                                                                                      <  8:201028/78‑1 (MQ=255)
                                                                                                                                                                                                                                                                            aaGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCatt                                                                                                                                                                                                      >  7:201028/1‑78 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                              
CGCGATTCACCAGGCTACGGACGCGACCATTGAATTCTTTATTCATGGGGCACTGTGCGTGGCCTATTCGGGTCAGTGCTACATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGCACAGATTATTTTGTGAATGCCCGTAAAGGCGATATTGGCGCGTTCGATTCGCCGAAATTTATCGGCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATC  >  NC_000913/1507230‑1507771

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