Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,635,9740TG54.3% ‑1.5 / 37.4 35I132L (ATC→CTC) nohAputative prophage DNA‑packaging protein NohA
Reads supporting (aligned to +/- strand):  ref base T (8/8);  new base G (10/9);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.90e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTCCCTGGCTGGAGATCTGTCTCGCTGGCCTGCCGCAGTTCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCAT  >  NC_000913/1635731‑1636219
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
cATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTc                                                                                                                                                                                                                                                 <  1:25611/250‑1 (MQ=255)
          aGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTgg                                                                                                                                                                                                                                      <  1:4648/251‑1 (MQ=255)
              gAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAAt                                                                                                                                                                                                                                  <  1:104403/251‑1 (MQ=255)
                                                gttATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAg                                                                                                                                                                                                <  1:76972/251‑1 (MQ=255)
                                                          ccGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGttt                                                                                                                                                                                      <  2:78800/250‑1 (MQ=255)
                                                              ataACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTccac                                                                                                                                                                                  >  1:301394/1‑251 (MQ=255)
                                                                                                        cgGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTc                                                                                                                                        <  1:166894/251‑1 (MQ=32)
                                                                                                                     tCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTGCGCAcgc                                                                                                                           >  1:364661/1‑251 (MQ=31)
                                                                                                                         cAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCCCGCGTa                                                                                                                        >  2:274574/1‑250 (MQ=9)
                                                                                                                           aCCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCAACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGt                                                                                                                     >  2:1657/1‑251 (MQ=9)
                                                                                                                                   aTCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTGCGCAGGCGTAAGTCGATGGCg                                                                                                             >  2:126968/1‑251 (MQ=14)
                                                                                                                                                 cgcgGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGCCTCTGGCCTTCTTCAGTTCCTGTGCGTCGGCCTCCGCACGCGTCAGTCGCCGGCGTTCGTCCTccat                                                                                                 >  2:206981/1‑246 (MQ=11)
                                                                                                                                                        gctTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGTCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCctg                                                                                         <  1:244950/250‑1 (MQ=9)
                                                                                                                                                         ctTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAACACGCCGCTGCACCGACCGGGGGATCCCGTCGAGAATACTGGCCATTTCACCTGCGATCCGCGACAGCACGAACGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTCCctggc                                                                                       >  2:64879/1‑251 (MQ=9)
                                                                                                                                                            tGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTCCctggctg                                                                                     <  1:209319/250‑1 (MQ=31)
                                                                                                                                                                         gatgatATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGCAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAACGTACAGAATGCGGATTCAACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCATGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATCGTTCCTGGCTGGAGATCTGCCTCGc                                                                       >  1:20676/1‑251 (MQ=11)
                                                                                                                                                                          atgatATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTCCCTGGCTGGAGCTCTGTCTCGCt                                                                      >  1:103240/1‑251 (MQ=14)
                                                                                                                                                                          atgatATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGATATCTGCCTCGCt                                                                      <  1:349592/251‑1 (MQ=255)
                                                                                                                                                                                   cGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTACGGAAAACGCCGCTGCACCGACAGAGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTCCCTGGCTGGAGATCTGTCTCGCTGGCCTGCCg                                                             >  2:157655/1‑251 (MQ=17)
                                                                                                                                                                                    gTTTCAGGAAATCAACATGTCGGATTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGc                                                            >  2:304900/1‑251 (MQ=255)
                                                                                                                                                                                        cAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGTCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCTTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCCGCCGCAGtt                                                        <  1:304900/251‑1 (MQ=11)
                                                                                                                                                                                                      gTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTCCCTGGCTGGAGATCTGTCTCGCTGGCCTGCCGCAGTTCTTCAACTTCCCgg                                          >  1:186550/1‑251 (MQ=34)
                                                                                                                                                                                                         ggTTTTACAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGgcgc                                       >  1:120726/1‑251 (MQ=11)
                                                                                                                                                                                                            tttCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACATAATGCGGTTTCCACCACTTCATCGGAGGCTCTGGCATTCTTCAGTTCCTGTGCGTCGTCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCGGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAg                                     <  1:61614/250‑1 (MQ=11)
                                                                                                                                                                                                            tttCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAg                                     <  1:159384/250‑1 (MQ=255)
                                                                                                                                                                                                               ccAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGGTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCttt                                 >  2:290134/1‑251 (MQ=11)
                                                                                                                                                                                                                    tCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGGCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCAGTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGt                             <  2:69527/250‑1 (MQ=11)
                                                                                                                                                                                                                     ccGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGtt                            <  1:206981/250‑1 (MQ=255)
                                                                                                                                                                                                                      cGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTtct                          >  1:84392/1‑251 (MQ=255)
                                                                                                                                                                                                                         aacaCGCCGCTGCCCCGACAGGGGGAGCCCGCCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCCCTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCCTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTTGCTTGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTtcgcta                       <  2:312985/248‑5 (MQ=11)
                                                                                                                                                                                                                            aCGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTGCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGGTCTCAAttt                    <  2:60311/251‑1 (MQ=11)
                                                                                                                                                                                                                               ccgcTGCGCCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCGCCGGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCGCTTCAGCGTAGTCTCTGGCATTTTTCAGCTCCGGTGCGTCGGCTTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTCCCTGGCTGGATATCTGGCTCGCTGGCTTGCCGCAGTTCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAg                 <  1:63086/251‑1 (MQ=11)
                                                                                                                                                                                                                                   tGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTCCCTGGCTGGAGATCTGTCTCGCTGGCCTGCCGCAGTTCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATc             >  2:369016/1‑251 (MQ=34)
                                                                                                                                                                                                                                            aGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGc    >  1:138619/1‑251 (MQ=255)
                                                                                                                                                                                                                                              ggggAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCAt  >  2:352882/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
CATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTCCCTGGCTGGAGATCTGTCTCGCTGGCCTGCCGCAGTTCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCAT  >  NC_000913/1635731‑1636219

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