Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,536,304 C→T G105D (GGC→GAC)  narV ← nitrate reductase Z subunit gamma

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,536,3040CT98.5% 216.1 / NA 68G105D (GGC→GAC) narVnitrate reductase Z subunit gamma
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (32/35);  minor base G (0/1);  total (32/36)
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.96e-01

CGGAAGACGAACGCTACGCCGTTGAGCATCTCTGACGATCCACCACGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGCCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGGTGTACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGA  >  NC_000913/1536061‑1536550
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
cGGAAGACGAACGCTACGCCGTTGAGCATCTCTGACGATCCACCACGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCGGTCCTGCGCCGTCAATCa                                                                                                                                                                                                                                                  >  1:506128/1‑250 (MQ=255)
 ggAAGACGAACGCTACGCCGTTGAGCATCTCTGACGATCCACCACGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAAc                                                                                                                                                                                                                                                >  1:56670/1‑251 (MQ=255)
  gAAGACGAACGCTACGCCGTTGAGCATCTCTGACTATCCACCACGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACTCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACg                                                                                                                                                                                                                                               <  1:353538/251‑1 (MQ=255)
      aCGAACGCTACGCCGTTGAGCATCTCTGACGATCCACCACGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCaa                                                                                                                                                                                                                                           <  2:266830/251‑1 (MQ=255)
       cGAACGCTACGCCGTTGAGCATCTCTGACGATCCACCACGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAac                                                                                                                                                                                                                                          <  1:323022/251‑1 (MQ=255)
         aaCGCTACGCCGTTGAGCATCTCTGACGATCCACCACGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACTTCAacac                                                                                                                                                                                                                                        <  2:56670/251‑1 (MQ=255)
                cGCCGTTGAGCATCTCTGACGATCCACCACGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCa                                                                                                                                                                                                                                  >  2:101648/1‑250 (MQ=255)
                  ccGTTGAGCATCTCTGACGATCCACCACGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAAt                                                                                                                                                                                                                               >  2:140223/1‑251 (MQ=255)
                         gCATCTCTGACGATCCACCACGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATAccgccg                                                                                                                                                                                                                        <  1:52640/251‑1 (MQ=255)
                                  aCGATCCACCACGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATc                                                                                                                                                                                                               <  2:209833/251‑1 (MQ=255)
                                  aCGATCCACCACGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATc                                                                                                                                                                                                               <  1:396691/251‑1 (MQ=255)
                                  aCGATCCACCACGGAAAGTGACAATGCTTGGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATc                                                                                                                                                                                                               <  2:506128/251‑1 (MQ=255)
                                                aaaGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATAAGGATACTGTGCTGACAACGGTATGGTACTTAGCCACAACAGACACTGGATCAGCAAACTGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTAATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGttt                                                                                                                                                                                                 >  1:373392/1‑251 (MQ=255)
                                                 aaGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTc                                                                                                                                                                                                <  1:556865/251‑1 (MQ=255)
                                                 aaGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTc                                                                                                                                                                                                >  1:83098/1‑251 (MQ=255)
                                                      aCAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGc                                                                                                                                                                                           >  1:544386/1‑251 (MQ=255)
                                                           gCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTGGTTTCGCTGCCACt                                                                                                                                                                                       <  1:449277/250‑1 (MQ=255)
                                                              gtGCGCCCAACCCACCAGATTCATCATTTCGCTACCAGCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGGGCCACAGCAGTCCTGCGCCGGCAAGCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCAGGAGTTGTTGCGCTGCCACTGGCa                                                                                                                                                                                   <  2:247098/250‑1 (MQ=255)
                                                                 cgcCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCaaa                                                                                                                                                                                 >  1:475114/1‑250 (MQ=255)
                                                                    ccAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAAcc                                                                                                                                                                             >  1:255505/1‑251 (MQ=255)
                                                                         ccaccaGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGCCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCAcgcg                                                                                                                                                                        <  1:311485/251‑1 (MQ=255)
                                                                         ccaccaGCTTCATCATGTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCAGGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCACCGGCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCAcgcg                                                                                                                                                                        <  1:92477/251‑1 (MQ=255)
                                                                          caccaGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGt                                                                                                                                                                       >  1:34582/1‑251 (MQ=255)
                                                                                    atcatTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTg                                                                                                                                                             >  2:507050/1‑251 (MQ=255)
                                                                                       atTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCAAGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCCTCGCCATGAGTTGTTTCGCTGCCCCTGACAAAACCCACGCGTACATCCAGTgcgg                                                                                                                                                          >  1:329622/1‑251 (MQ=255)
                                                                                          tCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTTCCACTGGCAAAAACCACGCGTACATCCAGTgcggcg                                                                                                                                                        <  2:54195/250‑1 (MQ=255)
                                                                                              tACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAAc                                                                                                                                                   >  1:574922/1‑251 (MQ=255)
                                                                                                    cAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGcc                                                                                                                                              <  2:424988/250‑1 (MQ=255)
                                                                                                          aCTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAg                                                                                                                                       <  1:528058/251‑1 (MQ=255)
                                                                                                          aCTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAg                                                                                                                                       <  1:317602/251‑1 (MQ=255)
                                                                                                                    aaaCGGTATGGTACTTAGCCCAAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCgaaga                                                                                                                             <  2:155862/251‑1 (MQ=255)
                                                                                                                       cGGTAAGGTACTTCGCACCAACAGACACTGGATCAGAACAATGCTCATGATGATGATATCCGGCGTAGTGGAAAAGGCACGAACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGACGTCAATCAACGTCACCACACCGCAAATACCGCCGAGCACCATCGCAATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAAAAAGTGCCCGAAGaaa                                                                                                                           >  2:280428/1‑250 (MQ=255)
                                                                                                                         gTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATc                                                                                                                        <  2:493509/251‑1 (MQ=255)
                                                                                                                          tATGGTACTTAGCCCCCACAGACACTGGATCTGCAAAATGCTCATGATGATGATATCCGGCGTCGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGCGCGGAGTTAACATGCAGAACAAGTGCCCGAAGAAAATcc                                                                                                                       >  1:422663/1‑251 (MQ=255)
                                                                                                                                     gccacAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACATCAGTCCTGCGCCGTCAAGCAACTTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCg                                                                                                            <  2:94880/247‑1 (MQ=255)
                                                                                                                                       cccAACAGACACGGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGat                                                                                                          <  2:416662/251‑1 (MQ=255)
                                                                                                                                              gACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAAc                                                                                                   <  2:116561/251‑1 (MQ=255)
                                                                                                                                               acacTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAAc                                                                                                   >  2:503806/1‑250 (MQ=255)
                                                                                                                                               acacTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGAGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAAc                                                                                                   >  2:32378/1‑250 (MQ=255)
                                                                                                                                                     gATCAGCAAAATGCTCATGATGATGATATCCGTCGTCGTGGAAGTGTCCCGCACGCGCTTATTTGTCAGCCTGCGCCTCAGCAGTCCTGCGCCGTCAATCATCGTCAACTCACCGCAAATACCGCCGAGCACCTTCGCCATGAGTTGGTTCGCTTCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAACCCCCAACACGCCGATATGGTACAGATTcc                                                                                            <  2:481702/251‑2 (MQ=255)
                                                                                                                                                         aGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGAcc                                                                                         >  1:304521/1‑250 (MQ=255)
                                                                                                                                                          gCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCat                                                                                       <  1:507050/251‑1 (MQ=255)
                                                                                                                                                           cAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCat                                                                                       >  1:134373/1‑250 (MQ=255)
                                                                                                                                                                 gCTCATGAGGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAAGGCCGATATGGAACAGATTCGACCATATCAcc                                                                                 <  2:31911/250‑1 (MQ=255)
                                                                                                                                                                    cATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATcc                                                                             <  1:140223/251‑1 (MQ=255)
                                                                                                                                                                       gatgatgatATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCcgc                                                                          >  1:246041/1‑251 (MQ=255)
                                                                                                                                                                         tgatgatATCCGGCGTAGTGGACGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAACTACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCACAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAACATGCCGCTATGGAACAGATTCGACCATATCACCATCCCGCGt                                                                        >  2:506234/1‑251 (MQ=255)
                                                                                                                                                                          aatgatATCCGGCGTAGTGGAAGGGGCACGCACGCGCTGATTTGTCAGACTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACCCACCGCAAATCCCGCCGAGCACCAGCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTTCAGCCAGTGCGGCGGTAACATTCCGAGCAAGTGCCCGAAGAAAATCCCCAAAAAGCCGAGAATGAACTGACGCGGCCATATCAGCATCCCGCGt                                                                        <  2:373392/249‑1 (MQ=255)
                                                                                                                                                                               tatCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAACATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCg                                                                  >  2:514427/1‑251 (MQ=255)
                                                                                                                                                                                   cGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACAGCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCAGATCACCATCCCGCGTTTATCGAGCa                                                              <  2:34582/251‑1 (MQ=255)
                                                                                                                                                                                      ggTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCAttt                                                           <  2:232069/250‑1 (MQ=255)
                                                                                                                                                                                         aGTGGAAGCGGCACGCACGCGCTGATTTGTCAGCCTGCGACACCGCAGTCCTGCGCCGTCAATCACCGTCACCACAACGCAAATACCGCCGCGCACCATCGCCCTCAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGCACATCCCCAAAATGCCGTTATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGAc                                                        >  2:341796/1‑251 (MQ=255)
                                                                                                                                                                                                  ggCACGCACGCGCTGATTTGTCAGCATGCGCCACAGCAGTCCTGCGACGTCAATCAACGTCAACACACCGCAAATACCGCCGCGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCATTGCGGCGTTAACATGCCGCACAGGGGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATAACACCATCCCGCGTTTATCTAGCATTTGACTTGgggagg                                               >  2:470704/1‑251 (MQ=255)
                                                                                                                                                                                                     acgcacgcGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGgcgc                                            >  1:201533/1‑251 (MQ=255)
                                                                                                                                                                                                       gcacgcGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCTACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAACTCCCCAAAATGCCGATATGGCACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGgcgcgc                                          >  2:135502/1‑251 (MQ=255)
                                                                                                                                                                                                           gcgcTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAggg                                     >  1:324869/1‑250 (MQ=255)
                                                                                                                                                                                                                 aTTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGGTGTACt                                >  2:480550/1‑251 (MQ=255)
                                                                                                                                                                                                                      gcAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACCCACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGGTGTACTGCCCg                           <  2:513168/250‑1 (MQ=255)
                                                                                                                                                                                                                                cctcagcagTCCTGCGCCTTCAATCAACGTCACCACACCGCAAATACCGCCGAGCACCAGCGCCATGTGTTGTGTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATGCGACCATATCTCCATCCCGCGGTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGGTGTACGGCCCGTATGCGTAgc                 <  2:234949/248‑1 (MQ=255)
                                                                                                                                                                                                                                ccgcagcagTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCAGCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGGTGTACTGCCCGTAGTCGTAgc                 <  2:134373/248‑1 (MQ=255)
                                                                                                                                                                                                                                  acaGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACAGCCGAGAACCATCGCCATGAGTTGTTTCGCTGCCACCGGCAAAAACCACGAGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGGTGTACTGCCCGTAGTCGTAgcgc               >  1:559150/1‑251 (MQ=255)
                                                                                                                                                                                                                                   cagcagTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGGTGTACTGCCCGTAGTCGTAgcgc               >  2:4217/1‑250 (MQ=255)
                                                                                                                                                                                                                                   cagcagTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGGTGTACTGCCCGTAGTCGTAGCGca              <  2:246041/251‑1 (MQ=255)
                                                                                                                                                                                                                                   cagcagTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATGCGACCATATCTCCATCCCGCGTTTATCGAGCATTTGACTGGAGGAGGCGCGCCAGGTGTACTGCCCGTAGTCGTAGCGca              <  2:574922/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cagTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGGTGTCCTGCCCGTAGTCGTAGCGcagcc           >  2:201233/1‑251 (MQ=255)
                                                                                                                                                                                                                                       agTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGGTGTACTGCCCGTAGTCGTAGCGcagcc           >  2:105618/1‑250 (MQ=255)
                                                                                                                                                                                                                                          ccTGCTCCGTCAAGCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTGCGCTGCCACTGGCAAAAACCCCGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTAGCGAGCATGTGACTTGAGGAGGCGCGCCAGGTGTACTGCCCGTATTCGTAGCGCAGCCAGCt       <  2:380006/251‑1 (MQ=255)
                                                                                                                                                                                                                                          ccTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGGTGTACTGCCCGTAGTCGTAGCGCAGCCAGCt       <  2:201533/251‑1 (MQ=255)
                                                                                                                                                                                                                                               gccgtcAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGGTGTACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCgc  >  1:227549/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
CGGAAGACGAACGCTACGCCGTTGAGCATCTCTGACGATCCACCACGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGCCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGGTGTACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGA  >  NC_000913/1536061‑1536550

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