Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 917,078 T→G D202E (GAT→GAG ybjD → putative OLD family ATP‑dependent endonuclease; DUF2813 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913917,0780TG98.4% 222.0 / NA 63D202E (GAT→GAGybjDputative OLD family ATP‑dependent endonuclease; DUF2813 family protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (34/28);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGTTTAATGCCGGTGCTGCGCTTGCGTGATGCCCGTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGATGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTTTGTTGCAGGCAAAAGGCAC  >  NC_000913/916847‑917308
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                      
gaCACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGTTTAATGCCGGTGCTGCGCTTGCGTGATGCCCGTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAgggg                                                                                                                                                                                                                                     <  1:66351/235‑1 (MQ=255)
  ccttGCTCAGTTGTCTCGATAAAGGCGGACAACCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGTTTAATGCCGGTGCTGCGCTTGCGTGATGCCCGTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCa                                                                                                                                                                                                                                   <  2:244659/232‑1 (MQ=255)
    cTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGTTTAATGCCGGTGCTGCGCTTGCGTGATGCCCGTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGATGGGCAGAt                                                                                                                                                                                                                                >  1:487865/1‑236 (MQ=255)
     tGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGTTTAATGCCGGTGCTGCGCTTGCGTGATGCCCGTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGAt                                                                                                                                                                                                                                >  2:268031/1‑235 (MQ=255)
     tGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGTTTAATGCCGGTGCTGCGCTTGCGTGATGCCCGTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGAt                                                                                                                                                                                                                                <  2:599585/235‑1 (MQ=255)
       cgcAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGTTTAATGCCGGTGCTGCGCTTGCGTGATGCCCGTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGCCAGATTCg                                                                                                                                                                                                                             >  1:198543/1‑236 (MQ=255)
       cgcAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGTTTAATGCCGGTGCTGCGCTTGCGTGATGCCCGTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCACCCCGATTTCTTCGCCCGTGGGTTTTCTTCCCATCCGCAACATCTTCCTGAGGGGCAGATTCg                                                                                                                                                                                                                             >  2:409564/1‑236 (MQ=255)
           gTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGTTTAATGCCGGTGCTGCGCTTGCGTGATGCCCGTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGAt                                                                                                                                                                                                                                <  1:409564/229‑1 (MQ=255)
                 tcGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGTTTAATGCCGGTGCTGCGCTTGCGTGATGCCCTTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACtt                                                                                                                                                                                                                   >  2:174664/1‑236 (MQ=255)
                   gATAAAGACGGACATCCGATGTATGTCGAGGATATGAACGATCAGGCACGCCATCTGGTGCGTTTAATGCCGGTGCTGCGCTTGCGTGATGCCCGTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTc                                                                                                                                                                                                                 <  1:590915/236‑1 (MQ=255)
                      aaaGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGTTTAATGCCGGTGCTGCGCTTGCGTGATGCCCGTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGc                                                                                                                                                                                                              <  1:139264/236‑1 (MQ=255)
                                                        aCGATCAGGCACGCCATCTGGTGCGTTTAATGCCGGTGCTGCGCTTGCGTGCTGCCCGTTTTATTCGCCGTATTCGTAACGGCACGGTGCCAAACGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGCTTCCCTCGCCCGGCAGTTATCCTCCCATCAGCAAAATACACATGAGGGGCAGAATCGTCATGGACCTTCTGCCCTGCTCAAGCTGCTGTAGCCTTATGtctataag                                                                                                                                                                            >  2:312682/1‑231 (MQ=255)
                                                          gATCAGGCACGCGATATGGATCGTTTAACGCCGGTGCTGCGCTAGCGTGATGCCCGTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCa                                                                                                                                                                          <  2:629207/236‑1 (MQ=255)
                                                          gATCAGGCACGCCATCTGGTGCGTTTAATGCCGGTGCTGCGCTTGCGTGATGCCCGTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCa                                                                                                                                                                          <  1:268031/236‑1 (MQ=255)
                                                               ggCACGCCATCTGGTGCGTTTAATGCCGGTGCTGCGCTTGCGTGATGCCCGTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAggggg                                                                                                                                                                     >  1:716861/1‑236 (MQ=255)
                                                                                gTTTAATGCCGGTGCTGCGCTTGCGTGATGCCCGTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGata                                                                                                                                                     <  1:299233/235‑1 (MQ=255)
                                                                                 tttAATGCCGGTGCTGCGCTTGCGTGATGCCCGTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGTGGCCGGACAGGCGCGATATc                                                                                                                                                   >  1:578399/1‑236 (MQ=255)
                                                                                  ttAATGCCGATGCTGCGCTTGCGTGATGCCCGTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGata                                                                                                                                                     >  2:21230/1‑233 (MQ=255)
                                                                                  ttAATGCCGATGCTGCGCTTGCGTGATGCCCGTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGata                                                                                                                                                     <  1:21230/233‑1 (MQ=255)
                                                                                                gcgcTTGCGTTGTGCCGGTTGTAGGCGCCGTATTCGTAACGGCCCGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATgcggcgg                                                                                                                                     <  2:198543/235‑1 (MQ=255)
                                                                                                gcgcTTGCGTGATGCCCGTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCGCCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATgcggcgg                                                                                                                                     <  1:751693/235‑1 (MQ=255)
                                                                                                     tGCGTGATGCCCTTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGc                                                                                                                                <  1:174664/235‑1 (MQ=255)
                                                                                                         tGATGCCCGTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCa                                                                                                                           <  2:716861/236‑1 (MQ=255)
                                                                                                                       tGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAAc                                                                                                                  >  1:432207/1‑231 (MQ=255)
                                                                                                                            cGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGAc                                                                                                                                                               <  1:438989/181‑1 (MQ=255)
                                                                                                                            cGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGAc                                                                                                                                                               >  2:438989/1‑181 (MQ=255)
                                                                                                                                tgcggaaCGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCGCGATTTCCTCGCCGGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTAt                                                                                                    <  2:611109/231‑1 (MQ=255)
                                                                                                                                     aacgggctggTGCCAAATGTCCCTAATGATGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTCTCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTgg                                                                                                <  2:1052/227‑1 (MQ=255)
                                                                                                                                      aCGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGat                                                                                              >  2:757168/1‑236 (MQ=255)
                                                                                                                                        ggCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGcta                                                                                             >  1:375966/1‑235 (MQ=255)
                                                                                                                                          cACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGAt                                                                                                                                                                                                                                >  2:202902/1‑102 (MQ=255)
                                                                                                                                          cACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGAt                                                                                                                                                                                                                                <  1:202902/102‑1 (MQ=255)
                                                                                                                                          cACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATatca                                                                                          >  1:673035/1‑236 (MQ=255)
                                                                                                                                          cACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCCCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGCTatc                                                                                           >  1:62724/1‑235 (MQ=255)
                                                                                                                                              gTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATatca                                                                                          <  1:520581/232‑1 (MQ=255)
                                                                                                                                               tGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTg                                                                                                           >  1:480226/1‑214 (MQ=255)
                                                                                                                                               tGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTg                                                                                                           <  2:480226/214‑1 (MQ=255)
                                                                                                                                               tGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGgcg                                                                                                        <  1:107056/217‑1 (MQ=255)
                                                                                                                                               tGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGgcg                                                                                                        >  2:107056/1‑217 (MQ=255)
                                                                                                                                                 ccAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGGAATCATCAACCg                                                                                   >  1:320572/1‑236 (MQ=255)
                                                                                                                                                 ccAAATGTCCCTAATGGGGAAGTCACCGCGCGCCAGCTCGATTGCCTCGCCCGTGAGGTATCCTCACATCCGAAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATatcatc                                                                                        <  1:27670/231‑1 (MQ=255)
                                                                                                                                                         cccTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATatc                                                                                           >  1:495059/1‑220 (MQ=255)
                                                                                                                                                         cccTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATATCATCAACCGGATGATTg                                                                           >  1:106808/1‑236 (MQ=255)
                                                                                                                                                         cccTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATATCATCAACCGGATGATTg                                                                           >  1:238328/1‑236 (MQ=255)
                                                                                                                                                         cccTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCAGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATatc                                                                                           <  2:495059/220‑1 (MQ=255)
                                                                                                                                                          cataagtggaaaGTCACCGCGGGCCAGCTAGATTTCCTGGCCCGTGAGTTAACCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATATCATCAACCGGATGATTga                                                                          <  2:487865/226‑1 (MQ=255)
                                                                                                                                                                gtgGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATATCATCAACCGGATGATTgaccgacc                                                                    >  2:364102/1‑236 (MQ=255)
                                                                                                                                                                         aCCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATATCATCAACCGGATGATTGACCGACCTGGTGGgcg                                                           >  1:588604/1‑236 (MQ=255)
                                                                                                                                                                         aCCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATATCATCAACCGGATGATTGACCGACCTGGTGGgcg                                                           >  1:371114/1‑236 (MQ=255)
                                                                                                                                                                           cgcgcgCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATATCATCAACCGGATGATTGACCGACCTGGCGGgcgc                                                          >  2:248702/1‑235 (MQ=255)
                                                                                                                                                                               cgcCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTa                                                     <  1:757168/236‑1 (MQ=255)
                                                                                                                                                                                 ccAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATATCATCAACCGGATGATTga                                                                          <  2:367397/213‑1 (MQ=255)
                                                                                                                                                                                 ccAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATATCATCAACCGGATGATTga                                                                          >  1:367397/1‑213 (MQ=255)
                                                                                                                                                                                              cTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCt                                      <  2:238328/236‑1 (MQ=255)
                                                                                                                                                                                               tCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGGCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGGATCGGGTTATTTTGc                                       >  2:122132/1‑234 (MQ=255)
                                                                                                                                                                                                cGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGGTATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGGTATTTTGCt                                      >  1:428567/1‑234 (MQ=255)
                                                                                                                                                                                                  cccGTGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGc                                  <  2:578399/236‑1 (MQ=255)
                                                                                                                                                                                                      tGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCaa                                                                                                                                                                                                            <  1:84299/62‑1 (MQ=255)
                                                                                                                                                                                                      tGAGTTATCCTCACATCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCaa                                                                                                                                                                                                            >  2:84299/1‑62 (MQ=255)
                                                                                                                                                                                                                  acaTCCGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTttgt                  >  1:715236/1‑236 (MQ=255)
                                                                                                                                                                                                                      ccGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTTTGTTGCa              >  1:287033/1‑236 (MQ=255)
                                                                                                                                                                                                                       cGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTTTGTTGCAg             >  1:510699/1‑236 (MQ=255)
                                                                                                                                                                                                                       cGCAAAATCTCTCTGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATATCATCAACCGGATGATTGACCGAACTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTTTGTTGCAg             >  2:330524/1‑236 (MQ=255)
                                                                                                                                                                                                                                  tctGAGGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTTTGTTGCAGGCAAAAGGcac  >  2:24787/1‑236 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                      
GACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGTTTAATGCCGGTGCTGCGCTTGCGTGATGCCCGTTTTATGCGCCGTATTCGTAACGGCACGGTGCCAAATGTCCCTAATGTGGAAGTCACCGCGCGCCAGCTCGATTTCCTCGCCCGTGAGTTATCCTCACATCCGCAAAATCTCTCTGATGGGCAGATTCGTCAGGGACTTTCCGCAATGGTACAGCTGCTTGAGCATTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTTTGTTGCAGGCAAAAGGCAC  >  NC_000913/916847‑917308

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