Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 569,160 G→T 11.6% D87Y (GAC→TAC) ‡ ybcK → DLP12 prophage; putative recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913569,1600GT11.6% 281.1 / 23.9 86D87Y (GAC→TAC) ‡ybcKDLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand):  ref base G (38/38);  new base T (5/5);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

TGATACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAA  >  NC_000913/569021‑569299
                                                                                                                                           |                                                                                                                                           
tGATACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTgg                                                                                                                                             >  1:1022465/1‑140 (MQ=255)
tGATACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTgg                                                                                                                                             <  2:1022465/140‑1 (MQ=255)
      cATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAaga                                                                                                                               >  4:10844/1‑148 (MQ=255)
                 aaGAGTTAGGATTAAGTGCATTCAAAGGAAAGTATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCgg                                                                                                                    >  1:320755/1‑148 (MQ=255)
                     tttAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGaa                                                                                                                <  2:320755/148‑1 (MQ=255)
                      ttAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCa                                                                                                                                  <  1:927658/129‑1 (MQ=255)
                      ttAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCa                                                                                                                                  >  2:927658/1‑129 (MQ=255)
                       tAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACttt                                                                                                                                    >  3:254121/1‑126 (MQ=255)
                       tAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACttt                                                                                                                                    <  4:254121/126‑1 (MQ=255)
                           aTTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACtt                                                                                                                                     >  1:247683/1‑121 (MQ=255)
                           aTTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACtt                                                                                                                                     <  2:247683/121‑1 (MQ=255)
                              aaGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGaa                                                                                                       <  3:10844/148‑1 (MQ=255)
                                      tCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACt                                                                                                                                      >  1:412931/1‑109 (MQ=255)
                                      tCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACt                                                                                                                                      <  2:412931/109‑1 (MQ=255)
                                       cAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCa                                                                                                                                  <  1:296669/112‑1 (MQ=255)
                                       cAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCa                                                                                                                                  >  2:296669/1‑112 (MQ=255)
                                         aaGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                                                                >  3:290225/1‑112 (MQ=255)
                                         aaGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                                                                <  4:290225/112‑1 (MQ=255)
                                         aaGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAAtt                                                                                            >  3:600624/1‑148 (MQ=255)
                                                 cATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAg                                                                                                                                         >  3:1113035/1‑95 (MQ=255)
                                                 cATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAg                                                                                                                                         <  4:1113035/95‑1 (MQ=255)
                                                   tGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGa                                                                                                            <  1:832341/122‑1 (MQ=255)
                                                   tGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGa                                                                                                            >  2:832341/1‑122 (MQ=255)
                                                                         gAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAgaga                                                                                                                             <  1:401741/83‑1 (MQ=255)
                                                                         gAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAgaga                                                                                                                             >  2:401741/1‑83 (MQ=255)
                                                                                        ataGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTAtt                                                                         >  3:156089/1‑120 (MQ=255)
                                                                                        ataGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTAtt                                                                         <  4:156089/120‑1 (MQ=255)
                                                                                         taGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTTATTTTGAATCACGGTATTGa                                                                       >  2:650915/1‑121 (MQ=255)
                                                                                         taGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGa                                                                       <  1:650915/121‑1 (MQ=255)
                                                                                             gCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGa                                                                       <  3:836923/117‑1 (MQ=255)
                                                                                             gCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGa                                                                       >  4:836923/1‑117 (MQ=255)
                                                                                                    tatataTTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAt                                                                                         <  4:1143514/92‑1 (MQ=255)
                                                                                                    tatataTTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAt                                                                                         >  3:1143514/1‑92 (MQ=255)
                                                                                                    tatataTTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtata                                                                <  3:587227/117‑1 (MQ=255)
                                                                                                    tatataTTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtata                                                                >  4:587227/1‑117 (MQ=255)
                                                                                                        tatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTc                                                                                                                      >  2:586524/1‑59 (MQ=255)
                                                                                                        tatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTc                                                                                                                      <  1:586524/59‑1 (MQ=255)
                                                                                                        tatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACa                                                   <  4:23335/126‑1 (MQ=255)
                                                                                                        tatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACa                                                   >  3:23335/1‑126 (MQ=255)
                                                                                                         atTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGa                                                                                                                 <  2:322789/63‑1 (MQ=255)
                                                                                                         atTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGa                                                                                                                 >  1:322789/1‑63 (MQ=255)
                                                                                                          ttGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGt                                                                                                                       >  3:731894/1‑56 (MQ=255)
                                                                                                          ttGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGt                                                                                                                       <  4:731894/56‑1 (MQ=255)
                                                                                                           tGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTg                                                                                           >  2:578455/1‑83 (MQ=255)
                                                                                                           tGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTg                                                                                           <  1:578455/83‑1 (MQ=255)
                                                                                                           tGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTAt                                                                          <  2:582146/100‑1 (MQ=255)
                                                                                                           tGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTAt                                                                          >  1:582146/1‑100 (MQ=255)
                                                                                                           tGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCttt                          >  2:691491/1‑148 (MQ=255)
                                                                                                                  aaCTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCtt                                                          <  1:21450/109‑1 (MQ=255)
                                                                                                                  aaCTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCtt                                                          >  2:21450/1‑109 (MQ=255)
                                                                                                                  aaCTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCtt                                                          <  1:101984/109‑1 (MQ=255)
                                                                                                                  aaCTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCtt                                                          >  2:101984/1‑109 (MQ=255)
                                                                                                                      aCATTGTTAGTTGAAAGTCTGGACAGTCTTTCAAGAGAAAAAGTCGGTGAAGCGa                                                                                                            <  3:384953/55‑1 (MQ=255)
                                                                                                                      aCATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGa                                                                                                            >  4:384953/1‑55 (MQ=255)
                                                                                                                      aCATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAt                                                                                                           <  1:2091/56‑1 (MQ=255)
                                                                                                                      aCATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAt                                                                                                           >  2:2091/1‑56 (MQ=255)
                                                                                                                      aCATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg                                                        <  3:1078246/107‑1 (MQ=255)
                                                                                                                      aCATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg                                                        >  4:1078246/1‑107 (MQ=255)
                                                                                                                      aCATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCat               <  1:691491/148‑1 (MQ=255)
                                                                                                                       cATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCtt                                                          <  3:847743/104‑1 (MQ=255)
                                                                                                                       cATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCtt                                                          >  4:847743/1‑104 (MQ=255)
                                                                                                                         ttGTTATTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGaa                                                                                                       <  1:29998/57‑1 (MQ=255)
                                                                                                                         ttGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACg                                                                              <  3:688215/82‑1 (MQ=255)
                                                                                                                         ttGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACg                                                                              >  4:688215/1‑82 (MQ=255)
                                                                                                                         ttGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCtata                                       <  4:225266/121‑1 (MQ=255)
                                                                                                                         ttGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCtata                                       >  3:225266/1‑121 (MQ=255)
                                                                                                                         ttGTTAGTTGAAAGTACGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGaa                                                                                                       >  2:29998/1‑57 (MQ=255)
                                                                                                                             tAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAg                                             >  1:945731/1‑111 (MQ=255)
                                                                                                                             tAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAg                                             <  2:945731/111‑1 (MQ=255)
                                                                                                                               cttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg                                                        <  4:782756/97‑1 (MQ=255)
                                                                                                                               cttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg                                                        >  3:782756/2‑98 (MQ=255)
                                                                                                                                ttGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACa                                                   <  2:640053/102‑1 (MQ=255)
                                                                                                                                ttGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACa                                                   >  1:640053/1‑102 (MQ=255)
                                                                                                                                 tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAt                                                                                         <  3:84970/63‑1 (MQ=255)
                                                                                                                                 tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAt                                                                                         >  4:84970/1‑63 (MQ=255)
                                                                                                                                 tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCttt                                                         <  4:695534/95‑1 (MQ=255)
                                                                                                                                 tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCttt                                                         >  3:695534/1‑95 (MQ=255)
                                                                                                                                 tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCt                                          <  4:886360/110‑1 (MQ=255)
                                                                                                                                 tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCt                                          >  3:886360/1‑110 (MQ=255)
                                                                                                                                   aaaGTCTGGACAGTCTTTCAAGAGAAAAAGTGGGTGAAGGGATTGAATGTCTGAAATTGATTTTGCATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATaaa  <  4:600624/148‑1 (MQ=255)
                                                                                                                                   aaaGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCtt                                                          <  2:981770/92‑1 (MQ=255)
                                                                                                                                   aaaGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCtt                                                          >  1:981770/1‑92 (MQ=255)
                                                                                                                                    aaGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATAc                                               <  2:23305/102‑1 (MQ=255)
                                                                                                                                    aaGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATAc                                               >  1:23305/1‑102 (MQ=255)
                                                                                                                                       tCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGa                                                                       <  4:50817/75‑1 (MQ=255)
                                                                                                                                       tCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGa                                                                       >  3:50817/1‑75 (MQ=255)
                                                                                                                                        cTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCa                                                                                <  2:956201/65‑1 (MQ=255)
                                                                                                                                        cTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCa                                                                                >  1:956201/1‑65 (MQ=255)
                                                                                                                                           gACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAAt                                                 <  4:1142909/93‑1 (MQ=255)
                                                                                                                                           gACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAAt                                                 >  3:1142909/1‑93 (MQ=255)
                                                                                                                                           |                                                                                                                                           
TGATACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAA  >  NC_000913/569021‑569299

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