Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,279,940 T→A 8.1% L26H (CTC→CAC)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,279,9400TA 8.1% 231.8 / 7.1 76L26H (CTC→CAC) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (28/40);  new base A (3/3);  total (32/44)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCGGCTCCTGCAGCTGGAAAATCTACGTCAAAAACGGTCTGGTCACCTGGGAAACCCA  >  NC_000913/1279793‑1280084
                                                                                                                                                   |                                                                                                                                                
gTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTtct                                                                                                                                                  >  8:42986/1‑148 (MQ=255)
    ccTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAAt                                                                                                                                              <  3:125115/148‑1 (MQ=255)
       tGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATAc                                                                                                                                            >  7:437201/1‑147 (MQ=255)
          ccTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAAc                                                                                                                                        <  4:31931/148‑1 (MQ=255)
               aaCATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTg                                                                                                                                    <  2:59094/147‑1 (MQ=255)
               aaCATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTg                                                                                                                                    <  4:120011/147‑1 (MQ=255)
                  aTTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTg                                                                                                                                <  3:291422/148‑1 (MQ=255)
                   ttAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTgg                                                                                                                               >  4:296087/1‑148 (MQ=255)
                       cAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGggagg                                                                                                                           <  3:401852/148‑1 (MQ=255)
                        aaTGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGggagga                                                                                                                          <  8:437201/148‑1 (MQ=255)
                         aTGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAggat                                                                                                                         <  7:214406/148‑1 (MQ=255)
                         aTGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAggat                                                                                                                         <  3:206068/148‑1 (MQ=255)
                           gTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGggagg                                                                                                                           <  4:232599/144‑1 (MQ=255)
                           gTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGggagg                                                                                                                           >  3:232599/1‑144 (MQ=255)
                             cGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCACAATACCAACCGTGACTGGGAggatgga                                                                                                                      >  7:76280/1‑147 (MQ=255)
                             cGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCACAATACCAACCGTGACTGGGAggatgga                                                                                                                      >  1:398512/1‑147 (MQ=255)
                             cGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGAAGGGCATGGCCAGCTTCACAATACCAACCGAGACTGGGAggatgca                                                                                                                      >  1:398508/1‑145 (MQ=255)
                                tttATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATc                                                                                                                  >  7:223265/1‑148 (MQ=255)
                                 ttATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAAACTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCg                                                                                                                 >  3:126574/1‑148 (MQ=255)
                                    tCAGAGGGCCGACAGGCTGCCACCGGAGAAAACCGATGTGTAAATTCCTGGACAGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCACAATACCAACCGTGACTGGGAGGATGGATATCGcc                                                                                                               <  2:398508/147‑1 (MQ=255)
                                    tCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCACAATACCAACCGTGACTGGGAGGATGGATATCGcc                                                                                                               <  8:76280/147‑1 (MQ=255)
                                    tCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCACAATACCAACCGTGACTGGGAGGATGGATATCGcc                                                                                                               <  2:398512/147‑1 (MQ=255)
                                     cAGAGGGCCGCCAGGCTCACACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAACCAGAAGGGTGAAACATTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGGGACTGGGAGGATGGATATCGcccg                                                                                                             >  4:49070/1‑146 (MQ=255)
                                         gggCCGAGAGGCTCCGACTGGAGAAAACCTATGAGTAAATTCGTGGACCGGTTTCGCTTCTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGtt                                                                                                         <  2:571/148‑1 (MQ=255)
                                          ggCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTg                                                                                                        <  3:378796/148‑1 (MQ=255)
                                               aCAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGgcagc                                                                                                   <  8:126410/148‑1 (MQ=255)
                                                cAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGgcagca                                                                                                  <  7:79713/148‑1 (MQ=255)
                                                  ggCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCGCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCAt                                                                                                 >  3:386154/1‑147 (MQ=255)
                                                       ccACAGGAGAAAACCGATGAGTAAATTCCTGGCCCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACaa                                                                                            <  7:100734/147‑1 (MQ=255)
                                                       ccACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACaa                                                                                            >  8:11849/1‑147 (MQ=255)
                                                       ccACAGGAGAAAACAGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACaa                                                                                            >  6:345288/1‑147 (MQ=255)
                                                         acaGTAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAAt                                                                                         <  3:296087/148‑1 (MQ=255)
                                                         acaGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAAt                                                                                         <  2:310140/148‑1 (MQ=255)
                                                         acaGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAAt                                                                                         >  2:234592/1‑148 (MQ=255)
                                                          caGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGatat                                                                                                                   <  7:70659/121‑1 (MQ=255)
                                                          caGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGatat                                                                                                                   >  8:70659/1‑121 (MQ=255)
                                                            ggAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGt                                                                                      <  2:5832/148‑1 (MQ=255)
                                                             gagaAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCg                                                                                                           <  6:232635/126‑1 (MQ=255)
                                                             gagaAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCg                                                                                                           >  5:232635/1‑126 (MQ=255)
                                                              agaAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTcc                                                                                    <  5:345288/148‑1 (MQ=255)
                                                              agaAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTcc                                                                                    <  2:263676/148‑1 (MQ=255)
                                                                 aaaCCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGc                                                                                  <  4:126574/147‑1 (MQ=255)
                                                                   aCCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGctc                                                                                >  2:408865/1‑147 (MQ=255)
                                                                      gATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTAcc                                                                            >  1:62027/1‑148 (MQ=255)
                                                                        tGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCa                                                                          >  7:309147/1‑148 (MQ=255)
                                                                          aGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCAc                                                                         <  5:136636/147‑1 (MQ=255)
                                                                                ttCGTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTaa                                                                  <  7:232712/148‑1 (MQ=255)
                                                                                  ccTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTaaa                                                                 <  1:408865/147‑1 (MQ=255)
                                                                                   cTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACt                                                               >  4:335378/1‑148 (MQ=255)
                                                                                   cTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACt                                                               <  7:42986/148‑1 (MQ=255)
                                                                                     ggACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTg                                                              >  6:61773/1‑147 (MQ=255)
                                                                                         cGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCg                                                         >  3:422103/1‑148 (MQ=255)
                                                                                           gTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCGGc                                                       <  7:11849/148‑1 (MQ=255)
                                                                                            tttCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGctc                                                                                <  5:418334/122‑1 (MQ=255)
                                                                                            tttCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGctc                                                                                >  6:418334/1‑122 (MQ=255)
                                                                                             ttCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCGGCTc                                                     >  7:101581/1‑148 (MQ=255)
                                                                                                 cTACTTCACACAGAAGGGTGAAACCTTTGCCGATGGGCAGGGCCAGCTTCTCACTACCAACCGTGACTGGGAGGAGGGATATCGCCAGCGTTGGCAGCATGACATAATCGTCCGCTCTACCCCCGGGGTAAACTGCACCGGCTCCTgg                                                 >  2:195362/1‑147 (MQ=255)
                                                                                                    cTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCGGCTCCTGCa                                                <  5:22598/146‑1 (MQ=255)
                                                                                                    cTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCGGCTCCTGCAGc                                              >  6:180386/1‑148 (MQ=255)
                                                                                                      tCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCGGCTCCTGCAGCTg                                            >  8:359080/1‑148 (MQ=255)
                                                                                                        aaaCAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCGGCTCCTGCAGCTGGa                                          <  5:61773/148‑1 (MQ=255)
                                                                                                              aaGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCCCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCGGCTCCTGCAGCTGGAAAATCt                                    <  4:386154/148‑1 (MQ=255)
                                                                                                                 ggTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCGGCTCCTGCAGCTGGAAAATCTACg                                 >  5:112631/1‑148 (MQ=255)
                                                                                                                 ggTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCGGCTCCTGCAGCTGGAAAATCTACg                                 <  3:404508/148‑1 (MQ=255)
                                                                                                                  gTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGt                                                                                      >  3:24584/1‑94 (MQ=255)
                                                                                                                  gTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGt                                                                                      <  4:24584/94‑1 (MQ=255)
                                                                                                                         cTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCGGCTCCTGCAGCTGGAAAATCTACGTCAAAAAc                         <  8:438263/148‑1 (MQ=255)
                                                                                                                               cGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCGGCTCCTGCAGCTGGAAAATCTACGTCAAAAACGGTCTg                   <  8:309147/148‑1 (MQ=255)
                                                                                                                                gATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCGGCTCCTGCAGCTGGAAAATCTACGTCAAAAACGGTCTgg                  <  4:422103/148‑1 (MQ=255)
                                                                                                                                    ggCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCGGCTCCTGCAGCTGGAAAATCTACGTCAAAAACGGTCTGGTCa               <  1:415932/147‑1 (MQ=255)
                                                                                                                                     gCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCGGCTCCTGCAGCTGGAAAATCTACGTCAAAAACGGTCTGGTCAcc             <  3:335378/148‑1 (MQ=255)
                                                                                                                                         ggCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTAcc                                                                            >  7:484732/1‑81 (MQ=255)
                                                                                                                                         ggCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTAcc                                                                            <  8:484732/81‑1 (MQ=255)
                                                                                                                                         ggCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCGGCTCCTGCAGCTGGAAAATCTACGTCAAAAACGGTCTGGTCACCTgg          >  3:383702/1‑147 (MQ=255)
                                                                                                                                         ggCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACAGGGTAAACTGCACCGGCTCCTGCAGCTGGAAAATCTACGTCAAAAACGGTCTGGTCACCTg           >  3:69761/1‑146 (MQ=255)
                                                                                                                                         ggCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACAGGGTAAACTGCACCGGCTCCTGCAGCTGGAAAATCTACGTCAAAAACGGTCTGGTCACCTg           <  4:69761/146‑1 (MQ=255)
                                                                                                                                            cAGCTTCTCAATACAAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCCCCGGCTCCTGCCGCGGGAACATCTACGTCACAAACGGTCTTGTACCCGGGGaac      >  8:335356/1‑147 (MQ=25)
                                                                                                                                                  ctcAATACCAACCGTGACTTGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCGGCTCCTGCAGCTGGAAAATCTACTTCAAAAACGGTCTGGTCACCTGGGAAACCca  <  3:36391/146‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCGGCTCCTGCAGCTGGAAAATCTACGTCAAAAACGGTCTGGTCACCTGGGAAACCCA  >  NC_000913/1279793‑1280084

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