Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,214,599 T→A 9.9% K293* (AAA→TAA)  bluF ← anti‑repressor for YcgE, blue light‑responsive; FAD‑binding; inactive c‑di‑GMP phosphodiesterase‑like EAL domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,214,5990TA 9.9% 267.5 / 17.8 81K293* (AAA→TAA) bluFanti‑repressor for YcgE, blue light‑responsive; FAD‑binding; inactive c‑di‑GMP phosphodiesterase‑like EAL domain protein
Reads supporting (aligned to +/- strand):  ref base T (36/37);  new base A (4/4);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

TTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTC  >  NC_000913/1214466‑1214725
                                                                                                                                     |                                                                                                                              
ttATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCaa                                                                                                                  <  1:840349/148‑1 (MQ=255)
             tAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGGTCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCg                                                                                                     >  2:630013/1‑148 (MQ=255)
               aTCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCtt                                                                                                                        >  1:409699/1‑127 (MQ=255)
               aTCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCtt                                                                                                                        <  2:409699/127‑1 (MQ=255)
                   aTTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTAAATCGCtaa                                                                                                                       <  3:614776/124‑3 (MQ=255)
                   aTTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTAAATCGCtaa                                                                                                                       >  4:614776/1‑122 (MQ=255)
                               ttttAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTc                                                                                   >  2:921035/1‑148 (MQ=255)
                                    aTTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTAAATCGCtaaa                                                                                                                      >  4:877305/1‑105 (MQ=255)
                                    aTTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTAAATCGCtaaa                                                                                                                      <  3:877305/108‑4 (MQ=255)
                                          tCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACtt                                                                                     >  1:476566/1‑135 (MQ=255)
                                          tCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACtt                                                                                     <  2:476566/135‑1 (MQ=255)
                                           cAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCaaa                                                                                                                 >  3:1013214/1‑106 (MQ=255)
                                           cAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCaaa                                                                                                                 <  4:1013214/106‑1 (MQ=255)
                                             ggCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTaa                                                                                                                              <  3:314433/91‑1 (MQ=255)
                                             ggCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTaa                                                                                                                              >  4:314433/1‑91 (MQ=255)
                                             ggCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACga                                                                     >  4:37503/1‑148 (MQ=255)
                                                 ggAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGGTCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACgatgat                                                                 <  1:630013/148‑1 (MQ=255)
                                                 ggAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACgatgat                                                                 >  4:644923/1‑148 (MQ=255)
                                                   aaGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTc                                                                                   <  1:1084074/128‑1 (MQ=255)
                                                   aaGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTc                                                                                   >  2:1084074/1‑128 (MQ=255)
                                                   aaGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGAttt                                                               <  1:921035/148‑1 (MQ=255)
                                                     gCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGt                                                             <  3:37503/148‑1 (MQ=255)
                                                             ggAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAAt                                                                            >  1:333222/1‑125 (MQ=255)
                                                             ggAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAACCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAAt                                                                            <  2:333222/125‑1 (MQ=255)
                                                               aGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCa                                                   >  3:692804/1‑148 (MQ=255)
                                                                   aCAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTaa                                                                                                                              >  1:835142/1‑69 (MQ=255)
                                                                   aCAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTaa                                                                                                                              <  2:835142/69‑1 (MQ=255)
                                                                    cAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACgatgat                                                                 <  3:953179/129‑1 (MQ=255)
                                                                    cAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACgatgat                                                                 >  4:953179/1‑129 (MQ=255)
                                                                     aaaCCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCAt                                             <  4:692804/148‑1 (MQ=255)
                                                                      aaCCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTc                                                                                   >  4:489522/1‑109 (MQ=255)
                                                                      aaCCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTc                                                                                   <  3:489522/109‑1 (MQ=255)
                                                                       aCCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTg                                           >  3:1028611/1‑148 (MQ=255)
                                                                        ccAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTgg                                          >  4:316073/1‑148 (MQ=255)
                                                                        ccAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTgg                                          <  3:316073/148‑1 (MQ=255)
                                                                          aGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGcc                                        >  1:752457/1‑148 (MQ=255)
                                                                           gCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCt                                       <  3:644923/148‑1 (MQ=255)
                                                                             aaaaCCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCAt                                             <  1:609630/140‑1 (MQ=255)
                                                                             aaaaCCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCttt                                     <  2:752457/148‑1 (MQ=255)
                                                                             aaaaCCGGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCAt                                             >  2:609630/1‑140 (MQ=255)
                                                                               aaCCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTaa                                                                                                                              >  3:383146/1‑57 (MQ=255)
                                                                               aaCCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTaa                                                                                                                              <  4:383146/57‑1 (MQ=255)
                                                                                aCCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTAtt                                  >  4:577806/1‑148 (MQ=255)
                                                                                  cTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTc                                >  1:207239/1‑148 (MQ=255)
                                                                                   tGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCgag                                                                                                   >  1:585303/1‑80 (MQ=255)
                                                                                   tGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCgag                                                                                                   <  2:585303/80‑1 (MQ=255)
                                                                                   tGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCa                               <  4:1028611/148‑1 (MQ=255)
                                                                                    gcgcCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGc                                                                                                                          <  3:1126009/56‑1 (MQ=255)
                                                                                    gcgcCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGc                                                                                                                          >  4:1126009/1‑56 (MQ=255)
                                                                                    gcgcCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCAt                              <  3:577806/148‑1 (MQ=255)
                                                                                      gcCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTa                            <  2:207239/148‑1 (MQ=255)
                                                                                          aaaTGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGAc                                                                                            <  3:1037474/80‑1 (MQ=255)
                                                                                          aaaTGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGAc                                                                                            >  4:1037474/1‑80 (MQ=255)
                                                                                          aaaTGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCAt                                             >  4:754227/1‑127 (MQ=255)
                                                                                          aaaTGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCAt                                             <  3:754227/127‑1 (MQ=255)
                                                                                                 cAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCa                               <  3:359088/134‑1 (MQ=255)
                                                                                                 cAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCa                               >  4:359088/1‑134 (MQ=255)
                                                                                                     tGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTc                                                                                   <  3:131782/78‑1 (MQ=255)
                                                                                                     tGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTc                                                                                   >  4:131782/1‑78 (MQ=255)
                                                                                                          cacTGATACCAGCAGCCTTTAGCGATTTAATTGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCAtt                                            <  3:291383/112‑1 (MQ=255)
                                                                                                          cacTGATACCAGCAGCCTTTAGCGATTTAATTGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCAtt                                            >  4:291383/1‑112 (MQ=255)
                                                                                                          cacTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAAtt                                                                           <  1:340598/81‑1 (MQ=255)
                                                                                                          cacTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAAtt                                                                           >  2:340598/1‑81 (MQ=255)
                                                                                                                     gcagcCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAagag                <  2:618363/129‑1 (MQ=255)
                                                                                                                     gcagcCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAagag                >  1:618363/1‑129 (MQ=255)
                                                                                                                      cagcCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTAttt                                 <  3:82505/111‑1 (MQ=255)
                                                                                                                      cagcCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTAttt                                 >  4:82505/1‑111 (MQ=255)
                                                                                                                      cagcCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTc  >  4:890793/1‑142 (MQ=255)
                                                                                                                      cagcCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTc  <  3:890793/142‑1 (MQ=255)
                                                                                                                          cTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTc                                <  2:888891/108‑1 (MQ=255)
                                                                                                                          cTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTc                                >  1:888891/1‑108 (MQ=255)
                                                                                                                          cTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAagag                <  1:665642/124‑1 (MQ=255)
                                                                                                                          cTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAagag                >  2:665642/1‑124 (MQ=255)
                                                                                                                           tttAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACgat                                                                    <  3:563580/71‑1 (MQ=255)
                                                                                                                           tttAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACgat                                                                    >  4:563580/1‑71 (MQ=255)
                                                                                                                             aaGCGATTAAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTc                                                                                   >  1:436053/2‑54 (MQ=255)
                                                                                                                             aaGCGATTAAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTc                                                                                   <  2:436053/53‑1 (MQ=255)
                                                                                                                             aaGCGATTAAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGc                                         <  4:370495/95‑1 (MQ=255)
                                                                                                                             aaGCGATTAAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGc                                         >  3:370495/2‑96 (MQ=255)
                                                                                                                                  aTTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACg                                                                      >  2:367207/1‑62 (MQ=255)
                                                                                                                                  aTTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACg                                                                      <  1:367207/62‑1 (MQ=255)
                                                                                                                                   tttAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAag                  >  1:621353/1‑113 (MQ=255)
                                                                                                                                   tttAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAag                  <  2:621353/113‑1 (MQ=255)
                                                                                                                                    ttAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGt                                                             >  3:810339/1‑69 (MQ=255)
                                                                                                                                    ttAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGt                                                             <  4:810339/69‑1 (MQ=255)
                                                                                                                                     |                                                                                                                              
TTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTC  >  NC_000913/1214466‑1214725

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