Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 566,742 A→C 100% pseudogene (59/93 nt) renD → pseudogene, DLP12 prophage;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913566,7420AC97.9% 5762.2 / 46.4 1996pseudogene (59/93 nt)renDpseudogene, DLP12 prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base A (22/22);  new base C (1060/889);  total (1083/913)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.55e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 1996 total aligned reads displayed.

ATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAG  >  NC_000913/566598‑566883
                                                                                                                                                |                                                                                                                                             
aTGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGa                                                                                                                                            <  5:241133/148‑1 (MQ=255)
 tGGGCTGTAGAGCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACg                                                                                                                                          <  2:168475/149‑1 (MQ=255)
   ggCGGTAGACCTCTAATTCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTGATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGtt                                                                                                                                        <  3:9150/149‑1 (MQ=255)
   ggCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGt                                                                                                                                         <  1:169143/148‑1 (MQ=255)
     cGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGtt                                                                                                                                        <  3:69753/147‑1 (MQ=255)
     cGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGc                                                                                                                                      >  6:41829/1‑149 (MQ=255)
       gTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCTAAGCTTAGTTCGGACTGAATGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTTCGCCGGACGTTGCcg                                                                                                                                    <  7:301764/149‑1 (MQ=255)
          gACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCcgcgc                                                                                                                                 >  3:171878/1‑148 (MQ=255)
          gACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCcgcgc                                                                                                                                 <  8:164007/149‑2 (MQ=255)
           aCCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCcgcgct                                                                                                                                <  5:335750/149‑3 (MQ=255)
             ctctAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGTACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTaa                                                                                                                              >  4:197186/1‑149 (MQ=255)
             ctctAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTaa                                                                                                                              >  3:110377/1‑149 (MQ=255)
             ctctAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTaa                                                                                                                              >  1:312807/1‑149 (MQ=255)
             ctctAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTaa                                                                                                                              >  1:26587/1‑149 (MQ=255)
              tctAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTaa                                                                                                                              >  6:203257/1‑148 (MQ=255)
              tctAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTaa                                                                                                                              >  2:370498/1‑148 (MQ=255)
              tctAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTaa                                                                                                                              >  8:356647/1‑148 (MQ=255)
              tctAAATCGTGCACAGGCTCTGGCGAACAACGCAGAAATCAAAGCAAAGTTCGTAATGAAAGGAGCAAGAGTATGAAGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGATTCTGAGCGACGGACGTTGCCGCGCTaa                                                                                                                              >  8:320934/1‑148 (MQ=255)
               ctAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTaa                                                                                                                              >  8:185138/1‑147 (MQ=255)
                tAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAAc                                                                                                                             <  2:147741/147‑1 (MQ=255)
                tAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAg                                                                                                                           >  5:88496/1‑149 (MQ=255)
                tAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAg                                                                                                                           >  8:213665/1‑149 (MQ=255)
                tAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAg                                                                                                                           >  3:345325/1‑149 (MQ=255)
                  aaTCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACa                                                                                                                            <  1:87223/146‑1 (MQ=255)
                  aaTCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGgc                                                                                                                         <  8:26227/149‑1 (MQ=255)
                       tGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCaa                                                                                                                     >  1:102130/1‑148 (MQ=255)
                         cacaGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAAca                                                                                                                   <  1:351892/148‑2 (MQ=255)
                          acaGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTACCAGGCGCAAcag                                                                                                                  >  7:199725/1‑146 (MQ=255)
                           caGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGt                                                                                                                 >  5:126563/1‑148 (MQ=255)
                            aGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTaa                                                                                                                              <  3:145268/134‑1 (MQ=255)
                             ggCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGCACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAc                                                                                                              <  5:62560‑M1/149‑2 (MQ=255)
                              gCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTaccc                                                                                                             >  8:297909/1‑146 (MQ=255)
                               ctctGGCGAAGATAGCAGAAATCAAATCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCTCTAACAGGCGCAACAGTAAcca                                                                                                            <  2:243288‑M1/149‑4 (MQ=255)
                                    gCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGTGACGCATAATAGCTTCTGTGAGCCGTACGTTGCCGCGCTAGCAGGCGCAACAGTAAccagcat                                                                                                        >  6:143295‑M1/1‑141 (MQ=255)
                                    gCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcat                                                                                                        >  4:10871‑M1/1‑141 (MQ=255)
                                    gCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagc                                                                                                          <  4:31107‑M1/146‑6 (MQ=255)
                                      gAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaa                                                                                                     >  2:39238‑M1/1‑139 (MQ=255)
                                         gATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAc                                                                                                              >  7:141671‑M1/1‑136 (MQ=255)
                                          aTCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcag                                                                                                 <  1:38723‑M1/149‑15 (MQ=255)
                                           tcacagaAATCAAAGCTCAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcagg                                                                                                <  8:373937‑M1/146‑16 (MQ=255)
                                           tCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcagg                                                                                                >  2:90757‑M1/1‑134 (MQ=255)
                                           tCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcagg                                                                                                >  1:198821‑M1/1‑134 (MQ=255)
                                           gcgCAGAAATCACAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcagg                                                                                                <  5:310509‑M1/148‑16 (MQ=255)
                                           gcgCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcagg                                                                                                <  7:375194‑M1/148‑16 (MQ=255)
                                             gCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTGAccagcataaatcaggcc                                                                                              >  6:93102‑M1/1‑132 (MQ=255)
                                               aGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgc                                                                                            >  2:305945‑M1/1‑130 (MQ=255)
                                                      aaaGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaa                                                                                      >  5:277537‑M1/1‑123 (MQ=255)
                                                      aaaGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatca                                                                                                  >  8:6759‑M1/1‑123 (MQ=255)
                                                      aaaGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAcca                                                                                                            <  8:13348‑M1/126‑4 (MQ=255)
                                                       aaGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaa                                                                                    >  2:196997‑M1/1‑122 (MQ=255)
                                                        aGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaat                                                                                   <  2:348500‑M1/149‑29 (MQ=255)
                                                         gCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatg                                                                                  <  2:218421‑M1/149‑30 (MQ=255)
                                                         gCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggcagcggctaaaatg                                                                                  >  7:324004‑M1/1‑120 (MQ=255)
                                                          cTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAAca                                                                                                                   <  4:87369/115‑2 (MQ=255)
                                                          cTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaactgg                                                                                 <  6:350810‑M1/149‑31 (MQ=255)
                                                           tAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatgg                                                                                 >  6:323693‑M1/1‑118 (MQ=255)
                                                             aGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcac                                                                              >  1:235955‑M1/1‑116 (MQ=255)
                                                             aGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGCTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggca                                                                               <  3:323516‑M1/148‑33 (MQ=255)
                                                              gTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacg                                                                             <  6:230851‑M1/149‑35 (MQ=255)
                                                               ttCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAAc                                                                                                                             <  5:215792/100‑1 (MQ=255)
                                                               ttCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgg                                                                            <  4:337555‑M1/149‑36 (MQ=255)
                                                               ttCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgg                                                                            >  7:11972‑M1/1‑114 (MQ=255)
                                                                    aCTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggc                                                                          <  5:369089‑M1/146‑38 (MQ=255)
                                                                     cTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGgcgc                                                                                                                       <  1:148477/100‑1 (MQ=255)
                                                                      tGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtc                                                                     >  4:230756‑M1/1‑107 (MQ=255)
                                                                      tGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtc                                                                     >  5:157532‑M1/1‑107 (MQ=255)
                                                                        aaaGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACGTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtctt                                                                   >  4:139752‑M1/1‑105 (MQ=255)
                                                                        aaaGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtctt                                                                   >  2:270769‑M1/1‑105 (MQ=255)
                                                                        aaaGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtctt                                                                   >  7:232630‑M1/1‑105 (MQ=255)
                                                                        aaaGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacggg                                                                           >  2:115751‑M1/1‑105 (MQ=255)
                                                                         aaGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttc                                                                  <  1:57965‑M1/149‑46 (MQ=255)
                                                                         aaGGACCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAAAGCTTCTGTGCGACGGACGTTGCCGCGCTAACAGGCCCAACAGAAAccagcaaaaatcaggccgcggctaaaatggcacgggcaggtcttc                                                                  >  4:362626‑M1/1‑104 (MQ=255)
                                                                          aGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttct                                                                 <  8:103383‑M1/149‑47 (MQ=255)
                                                                              gCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggt                                                              >  6:271953‑M1/1‑99 (MQ=255)
                                                                              gCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttct                                                                 >  7:280457‑M1/1‑99 (MQ=255)
                                                                               cAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggtta                                                            >  1:291756‑M1/1‑98 (MQ=255)
                                                                                aaGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttat                                                           <  5:1886‑M1/149‑53 (MQ=255)
                                                                                  gtgtATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcag                                                                        <  6:173020‑M1/134‑40 (MQ=255)
                                                                                     tATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaa                                                       <  3:234465‑M1/148‑57 (MQ=255)
                                                                                         aCGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaa                                                  <  7:357306‑M1/149‑62 (MQ=255)
                                                                                         aCGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaa                                                  <  7:41525‑M1/149‑62 (MQ=255)
                                                                                           gggCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaagg                                                >  8:39832‑M1/1‑86 (MQ=255)
                                                                                            ggCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggt                                               >  3:294544‑M1/1‑85 (MQ=255)
                                                                                             gCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtc                                              >  6:112895‑M1/1‑84 (MQ=255)
                                                                                             gCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtc                                              >  2:63003‑M1/1‑84 (MQ=255)
                                                                                             gCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtc                                              <  7:72680‑M1/149‑66 (MQ=255)
                                                                                               aaaGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctg                                            >  4:157796‑M1/1‑82 (MQ=255)
                                                                                               aaaGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctg                                            >  4:265212‑M1/1‑82 (MQ=255)
                                                                                                aaGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctgg                                           <  5:48028‑M1/149‑69 (MQ=255)
                                                                                                aaGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctgg                                           >  8:234440‑M1/1‑81 (MQ=255)
                                                                                                  gagGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggc                                          >  6:294677‑M1/1‑79 (MQ=255)
                                                                                                  gagGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggc                                          >  6:68311‑M1/1‑79 (MQ=255)
                                                                                                  gagGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcag                                                                        >  7:249748‑M1/1‑79 (MQ=255)
                                                                                                   agGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcg                                         >  5:103010‑M1/1‑78 (MQ=255)
                                                                                                   agGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcg                                         >  4:68098‑M1/1‑78 (MQ=255)
                                                                                                      cAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCGACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaa                                     <  3:266373/149‑1 (MQ=255)
                                                                                                            ttCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACa                                <  7:173610/148‑1 (MQ=255)
                                                                                                                   ccTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacggtgtattaccggt                        <  1:131727‑M1/149‑88 (MQ=255)
                                                                                                                    cTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacggtgtattaccggtt                       <  3:198316‑M1/149‑89 (MQ=255)
                                                                                                                                      ttCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTaaa     >  3:369029/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                             
ATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAG  >  NC_000913/566598‑566883

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