Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,507,500 G→T 5.1% K240N (AAG→AAT ydcP → putative peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,507,5000GT 5.1% 564.5 / 15.1 156K240N (AAG→AATydcPputative peptidase
Reads supporting (aligned to +/- strand):  ref base G (77/71);  new base T (4/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 = 4.70e-01

Only 100 of 157 total aligned reads displayed.

CGCGATTCACCAGGCTACGGACGCGACCATTGAATTCTTTATTCATGGGGCACTGTGCGTGGCCTATTCGGGTCAGTGCTACATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGCACAGATTATTTTGTGAATGCCCGTAAAGGCGATATTGGCGCGTTCGATTCGCCGAAATTTATCGGCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATC  >  NC_000913/1507230‑1507771
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                               
            gtcgcagggcgcgACCCATGATGTCTTTATTAAAGGGGCACTGTGCGGGGACTAGCCGGTTCAGTGCTACATCTTTAATGCGCAAACAGGGCGTAGAGCCAACCGTGGCGATGGCTCGCAGGCGTGCCGTGTGCCAGACACATTGAAAGACGAGCAGGTGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGAGAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACa                                                                                                                                                                                                                                                                 <  4:15860/266‑1 (MQ=255)
            ggCTACGGACGCGACCATTGAATTCTTTATTCATGGGGCACTGTGCGTGGCCTATTCGGGTCAGTGCTACATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACa                                                                                                                                                                                                                                                                 <  4:50941/275‑1 (MQ=255)
               tACGGACGCGACCATTGAATTCTTTATTCATGGGGCACTGTGCGTGGCCTATTCGGGTCAGTGCTACATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGCCTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTCCGCTCCTTCAAGATTGAAGGGCGt                                                                                                                                                                                                                                                                     >  1:140812/1‑268 (MQ=255)
                      gcgACCATTGAATTCTTTATTCATGGGGCACTGTGCGTGGCCTATTCGGGTCAGTGCTACATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAg                                                                                                                                                                                                                                                       <  4:491580/275‑1 (MQ=255)
                                              ggggCACTGTGCGTGGCCTATTCGGGTCAGTGCTACATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGcc                                                                                                                                                                                                                                 >  2:533467/1‑273 (MQ=255)
                                                   aCTGTGCGTGGCCTATTCGGGTCAGTGCTACATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGAACAGACTGCCCACCTCGGCGCGCTTATTGATGCTGGGTTACGGTCCTTCCAGATTGAAGGGGGTTAAAAAAATATGAGCTACGTGAAAAAAATAACAGcccgtgact                                                                                                                                                                                                                          >  4:1107/1‑269 (MQ=255)
                                                       tgCGTGGCCTATTCGGGTCAGTGCTACATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaaga                                                                                                                                                                                                                                            <  2:1107/253‑1 (MQ=255)
                                                       tgCGTGGCCTATTCGGGTCAGTGCTACATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCCCCGCCCATTATCGCCa                                                                                                                                                                                                                      >  2:219477/1‑275 (MQ=255)
                                                                 aTTCGGGTCAGTGCTACATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGACCAGGGGCGGGTGGTTTCCTCTGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGCTTGAAGGGCGTTACAAAGATAGGAGCTACGGCAAGAATATCCCAGCCAATTATCCCCAGATGCTTTAt                                                                                                                                                                                                            >  4:199299/1‑275 (MQ=255)
                                                                   tCGGGTCAGTGCTACATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGACTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGc                                                                                                                                                                                                          <  2:109065/275‑1 (MQ=255)
                                                                      ggTCAGTGCTACATTGCTCATGCGCAAACAGGGCTGAGCGCGAACCGTGGCGATTGCTCGCAGGCGGGCCGGTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCGTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCa                                                                                                                                                                                                        <  2:487330/274‑1 (MQ=255)
                                                                       gTCAGTGCTACATTTCTCATGCTCAAACAGGGCGTAGCGCCAACCGTGGCGCTTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGCAAGATAACGATCAGACTGCCAACCTCGGCGCCACGCATGATGCTTGTGTACGCTCCTTCACGGTTGAAGGGCCTTACAAAGCTATGGGCTACGTTAAGAATTACCCTGCCCGTTTACGCCAGGTACTTgtgtaccat                                                                                                                                                                                                      >  4:94566/1‑267 (MQ=255)
                                                                            tGCTACATTTCTCATGCGCAAACAGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTAACGGGCGTTACAAAGATATGAGCTACGTGCAGAATAACACCGCCCATTATCGCCAGATGCTTGCTGCCATTTTTg                                                                                                                                                                                                  >  4:73215/1‑274 (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)
                                                                                                                                           cGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGaa                                                                                                                                                                                                                                                                           >  5:1100057/1‑138 (MQ=255)
                                                                                                                                            gTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAgg                                                                                                                                                                <  2:74726/244‑1 (MQ=255)
                                                                                                                                              ttGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAg                                                                                                                                                                                                                                                              <  5:781580/148‑1 (MQ=255)
                                                                                                                                                      cacaTTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGa                                                                                                                                                                                                                                                                            <  7:391055/126‑1 (MQ=255)
                                                                                                                                                      cacaTTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGa                                                                                                                                                                                                                                                                            >  8:391055/1‑126 (MQ=255)
                                                                                                                                                          ttGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAgg                                                                                                                                                                                                                                                                         <  5:904823/125‑1 (MQ=255)
                                                                                                                                                          ttGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAgg                                                                                                                                                                                                                                                                         >  6:904823/1‑125 (MQ=255)
                                                                                                                                                             aaaGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCCCTTCAATCTTGAAggagcgta                                                                                                                                                                                                                                                                   >  6:896859/1‑121 (MQ=255)
                                                                                                                                                             aaaGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCCCTTCAATCTTGAAggagcgta                                                                                                                                                                                                                                                                   <  7:387115/127‑7 (MQ=255)
                                                                                                                                                             aaaGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCCCTTCAATCTTGAAggagcgta                                                                                                                                                                                                                                                                   >  8:387115/1‑121 (MQ=255)
                                                                                                                                                              aaGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAg                                                                                                                                                                                                                                                       >  6:183985/1‑139 (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)
                                                                                                                                                                           gCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCg                                                                                                                                                                                                                                   >  4:272899/1‑146 (MQ=255)
                                                                                                                                                                               gtggtTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCa                                                                                                                                                                                                                               <  2:272900/146‑1 (MQ=255)
                                                                                                                                                                                tggtTTCCTATGAAAAACATCTGCTGTCGAGGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGCACAGattatt                                                                                               <  4:533466/273‑1 (MQ=255)
                                                                                                                                                                                     tCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATTCCATTATTGAAGAACGTGGCGTACTGGCGCGCGCTTCATCAGGTCGTACTGAGCATTTCTTTGTTACATCGACTGAAAAGACTNNNCACCGTGGTAACAAAGATTAtttttttaa                                                                                        >  3:145293/1‑270 (MQ=255)
                                                                                                                                                                                      ccTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGat                                                                                                                                                                                                                                                            >  6:979422/1‑110 (MQ=255)
                                                                                                                                                                                      ccTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGat                                                                                                                                                                                                                                                            <  5:979422/110‑1 (MQ=255)
                                                                                                                                                                                         aTGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAggg                                                                                                                                                                                                                                                                        <  8:788030/95‑1 (MQ=255)
                                                                                                                                                                                         aTGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAggg                                                                                                                                                                                                                                                                        >  7:788030/1‑95 (MQ=255)
                                                                                                                                                                                         aTGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaagaa                                                                                                                                                                                                                                           <  5:285331/124‑1 (MQ=255)
                                                                                                                                                                                          tGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaa                                                                                                                                                                                                                                              <  5:165876/120‑1 (MQ=255)
                                                                                                                                                                                            aaaaaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAATCTTGAAgg                                                                                                                                                                                                                                                                         <  7:87318/91‑1 (MQ=255)
                                                                                                                                                                                            aaaaaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAATCTTGAAgg                                                                                                                                                                                                                                                                         >  8:87318/1‑91 (MQ=255)
                                                                                                                                                                                             aaaaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGc                                                                                                                                                                                                                        >  8:573483/1‑139 (MQ=255)
                                                                                                                                                                                             aaaaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGc                                                                                                                                                                                                                        <  7:573483/139‑1 (MQ=255)
                                                                                                                                                                                               aaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGa                                                                                                                                                                                                             >  8:1012206/1‑148 (MQ=255)
                                                                                                                                                                                               aaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGa                                                                                                                                                                                                             <  7:161892/148‑1 (MQ=255)
                                                                                                                                                                                                      gctgTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGt                                                                                                                                                                                                                                                                     >  8:11041/1‑85 (MQ=255)
                                                                                                                                                                                                      gctgTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGt                                                                                                                                                                                                                                                                     <  7:11041/85‑1 (MQ=255)
                                                                                                                                                                                                             aTGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGa                                                                                                                                                                                                                                                                            <  5:592837/71‑1 (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)
                                                                                                                                                                                                               gAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTga                                                                                                                                                                                                                                               >  7:577113/1‑98 (MQ=255)
                                                                                                                                                                                                               gAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTga                                                                                                                                                                                                                                               <  8:577113/98‑1 (MQ=255)
                                                                                                                                                                                                               gAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTgaagaa                                                                                                                                                                                             <  7:1012206/148‑1 (MQ=255)
                                                                                                                                                                                                                     tAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTgg                                                                                                                                                                                        <  5:473993/147‑1 (MQ=255)
                                                                                                                                                                                                                       aCGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaa                                                                                                                                                                                                                                              >  8:877607/1‑91 (MQ=255)
                                                                                                                                                                                                                       aCGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaa                                                                                                                                                                                                                                              <  7:877607/91‑1 (MQ=255)
                                                                                                                                                                                                                          aTCAGACTGCCAACCTCGGCTCGCTGATTCATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGTTGCCATTATTGAATAACGTGGCGATCTTCCGCGCCCTTCATCAGGTCGTACTGTACATTTCTTTGTTCCATCGCCGCAAAAGCCTTTCCACCGTGTGAGCTAAGCTTATGTTGTTAATGCCCGTAATGGCGATGTTGTCGTTGTCGCTTCGCCg                                                   >  4:209452/1‑275 (MQ=255)
                                                                                                                                                                                                                          aTCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTg                                                                                                                                                                                                              >  6:1083617/1‑120 (MQ=255)
                                                                                                                                                                                                                          aTCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCattat                                                                                                                                                                                                    >  6:1106294/1‑130 (MQ=255)
                                                                                                                                                                                                                             aGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCattat                                                                                                                                                                                                    <  7:19799/127‑1 (MQ=255)
                                                                                                                                                                                                                             aGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCattat                                                                                                                                                                                                    >  8:19799/1‑127 (MQ=255)
                                                                                                                                                                                                                                  gCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCt                                                                                                                                                                                                                <  6:553888/110‑1 (MQ=255)
                                                                                                                                                                                                                                  gCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTg                                                                                                                                                                                         <  6:772256/133‑1 (MQ=255)
                                                                                                                                                                                                                                   ccAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCCAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCNGGCGCGCGCTTCATCAGGTCGTACTTAACATTTCTTTGTTCCTTCGACAGACACGACTTTCCACCGTGGTAGCACAGATTNTTTTTTGAATGCCCTTAAAGTCGCCNTTGTCTCGTTCGATTCTCCTAACttttt                                           >  3:99536/1‑272 (MQ=255)
                                                                                                                                                                                                                                   ccAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCNGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGCACAGATTNAGTTGTGAATGCCCGTAAAGGCGAANTTGGGGCGTTCGATTCGCCGAAATTTTTc                                          >  3:91207/1‑275 (MQ=255)
                                                                                                                                                                                                                                    cAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCt                                                                                                                                                                                                                >  5:825248/1‑108 (MQ=255)
                                                                                                                                                                                                                                    cAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCt                                                                                                                                                                                                                <  6:825248/108‑1 (MQ=255)
                                                                                                                                                                                                                                     aaCCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAg                                                                                                                                                                                                                                                                          <  6:104839/49‑1 (MQ=255)
                                                                                                                                                                                                                                     aaCCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAg                                                                                                                                                                                                                                                                          >  5:104839/1‑49 (MQ=255)
                                                                                                                                                                                                                                       ccTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGAGTGAAGGGCGTTAAAAAGATATGAGCTACGTGAAGAATATTACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGCACAGATTATTTTGTGAATGCCCGTAAAGGCGATATTGGCGCGTTCGATTCGCCGAAATTTATCGGcc                                      <  2:379769/275‑1 (MQ=255)
                                                                                                                                                                                                                                             cgcgcTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGACACGACTTTCCACCGTGGTAGCACAGATTATTTTGTGAATGCCCGTAAAGGCGATATTGGCGCTTTCGATTCGCAAAATTTTCTCGGCCTGcacg                                >  4:247665/1‑272 (MQ=255)
                                                                                                                                                                                                                                                    attgatGCTGGTGGACGCACCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGCACAGATTATTTTGTGAATGCCCGTAAAGGCGATATTGGCGCGTTCGATTCGCCGAAATTTATCGGCCTGCCGGTAGGCGa                         <  2:247666/275‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttgatGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTAt                                                                                                                                                                                                                           <  8:543310/80‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttgatGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTAt                                                                                                                                                                                                                           >  7:543310/1‑80 (MQ=255)
                                                                                                                                                                                                                                                          gCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTg                                                                                                                                                                                         >  6:390934/1‑109 (MQ=255)
                                                                                                                                                                                                                                                           cTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTga                                                                                                                                                                                                 <  6:923970/100‑1 (MQ=255)
                                                                                                                                                                                                                                                             ggTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTg                                                                                                                                                        <  5:747669/139‑1 (MQ=255)
                                                                                                                                                                                                                                                             ggTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCt                                                                                                                                               >  5:36836/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                   cGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGATCATTTCTTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGCACAGATTATTTTGTGAATGCCCGTAAGGGCGATATTGGCCCGTTCGATTCGCCGAAATTTATCTGCCTGCCGGTATGCGAGTTTGTGaacttgtc          >  4:425684/1‑269 (MQ=255)
                                                                                                                                                                                                                                                                   cGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGCACAGATTATTTTGTGAATGCCCGTAAAGGCGATATTGGCGCGTTCGATTCGCCGAAATTTATCGGCCTGCCGGTAGGCGAATTAGTGAAAGTgg           >  2:306795/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                    gCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTgaagaa                                                                                                                                                                                             >  6:823452/1‑95 (MQ=255)
                                                                                                                                                                                                                                                                     cTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTg                                                                                                                                            >  5:670742/1‑143 (MQ=255)
                                                                                                                                                                                                                                                                     cTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTg                                                                                                                                            <  6:670742/143‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGtt                                                                                                                                          <  5:113838/145‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGtt                                                                                                                                          >  6:113838/1‑145 (MQ=255)
                                                                                                                                                                                                                                                                      tCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCattat                                                                                                                                                                                                    <  6:662064/86‑1 (MQ=255)
                                                                                                                                                                                                                                                                        cTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGCACAGATTATTTTGTGAATGCCCGTAAAGGCGATATTGGCGCGTTCGATTCGCCGAAATTTATCGGCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAg     >  4:237092/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)
                                                                                                                                                                                                                                                                           cAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGCACAGATTATTTTGTGCATGCCCGTAAAGGCGATATTGGAGCGTTCGATTCGCCGAAATTTATCGGCCTGCCGGTAGGCGAAGTAGTGTAAGTGGCGTAAGatc  >  2:161426/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                            aaGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCatt                                                                                                                                                                                                      >  7:201028/1‑78 (MQ=255)
                                                                                                                                                                                                                                                                            aaGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCatt                                                                                                                                                                                                      <  8:201028/78‑1 (MQ=255)
                                                                                                                                                                                                                                                                            aaGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCCTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGCACAGa                                                                                                    <  2:248065/176‑1 (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: