Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,214,867 G→A 6.2% A203A (GCC→GCT)  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,8670GA 6.2% 361.3 / 6.7 119A203A (GCC→GCT) 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 G (51/58);  new base A (4/3);  total (57/62)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

AGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAGTTCTTTATCAGATCCGTCAG  >  NC_000913/1214720‑1215015
                                                                                                                                                   |                                                                                                                                                    
aGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATgg                                                                                                                                                     <  1:581870/149‑1 (MQ=255)
aGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATgg                                                                                                                                                     <  1:506777/149‑1 (MQ=255)
  gTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCt                                                                                                                                                   >  2:1266274/1‑149 (MQ=255)
   ttCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTg                                                                                                                                                  >  1:178296/1‑149 (MQ=255)
      gTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCATTACGCTGCCCAACCGCTATGCCTGATg                                                                                                                                               >  2:194229/1‑149 (MQ=255)
       ttAACCAGGGTCATTGGTAATAGATTGATTGAAATCATTTTATCACCGATCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATgg                                                                                                                                              <  1:860391/149‑1 (MQ=255)
       ttAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATgg                                                                                                                                              >  2:738186/1‑149 (MQ=255)
        tAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATggg                                                                                                                                             <  1:461120/149‑1 (MQ=255)
        tAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATggg                                                                                                                                             <  1:1442652/149‑1 (MQ=255)
          aCCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCt                                                                                                                                           >  2:324313/1‑149 (MQ=255)
           ccAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCATCAGCCATAGCGGTTGGGcag                                                                                                                                          <  1:1299407/149‑3 (MQ=255)
           ccAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTa                                                                                                                                          >  1:1490663/1‑149 (MQ=255)
           ccAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTACCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTa                                                                                                                                          >  1:562420/1‑149 (MQ=255)
            cAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCATCAGCCATAGCGGTTGGGcagc                                                                                                                                         <  1:133689/149‑4 (MQ=255)
             agtgtaaTAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGGGCCATCGTTAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATc                                                                                                                                        <  2:976039/143‑1 (MQ=255)
             aGGGTCATAGTTAATATATTGATTGAAATCATTTTATCACCTAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTATTTTTTAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATc                                                                                                                                        <  2:278343/149‑1 (MQ=255)
             aGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATc                                                                                                                                        <  2:178296/149‑1 (MQ=255)
               ggTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCtt                                                                                                                                      >  1:784418/1‑149 (MQ=255)
                 tCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCa                                                                                                                                    >  2:1465687/1‑149 (MQ=255)
                  cATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCAt                                                                                                                                   <  1:1423055/149‑1 (MQ=255)
                         aaTAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTg                                                                                                                            <  2:70470/149‑1 (MQ=255)
                          aTAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGAGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGc                                                                                                                           <  1:119024/149‑1 (MQ=255)
                           tAGATTGATTGAAATCATTTTATAAACGAGCTCAAGTGCGTTTGACATCGTGAATGCAAGCGCATTACTTTTGAGATACTCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTgaa                                                                                                                          <  2:592233/149‑3 (MQ=255)
                              attgattgAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACaa                                                                                                                       >  2:1306973/1‑149 (MQ=255)
                               ttgattgaAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAAt                                                                                                                      <  1:830283/149‑1 (MQ=255)
                                   ttgaAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCt                                                                                                                  >  1:424463/1‑149 (MQ=255)
                                   ttgaAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGCAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCt                                                                                                                  >  2:61762/1‑149 (MQ=255)
                                    tgaAATCATTTTATCACCGATCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCtt                                                                                                                 <  2:522935/149‑1 (MQ=255)
                                    tgaAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCtt                                                                                                                 <  2:1569238/149‑1 (MQ=255)
                                      aaaTCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCa                                                                                                               >  2:44675/1‑149 (MQ=255)
                                       aaTCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCaa                                                                                                              >  1:1067040/1‑149 (MQ=255)
                                        aTCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCaaa                                                                                                             <  1:1352403/149‑1 (MQ=255)
                                           aTTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGc                                                                                                          >  1:1352765/1‑149 (MQ=255)
                                            ttttATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACaa                                                                                                                       >  2:926099/1‑135 (MQ=255)
                                            ttttATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACaa                                                                                                                       <  1:926099/135‑1 (MQ=255)
                                             tttATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGGCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGa                                                                                                        <  2:777933/149‑1 (MQ=255)
                                             tttATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGa                                                                                                        <  2:562420/149‑1 (MQ=255)
                                              ttATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGAt                                                                                                       <  2:752490/149‑1 (MQ=255)
                                                aTCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTa                                                                                                     <  1:588407/149‑1 (MQ=255)
                                                    ccGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATccg                                                                                                 >  2:579650/1‑149 (MQ=255)
                                                     cGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCTCTGTGTAGATTTCCCCGTCTTAACGATGCACAACAGCTATGTCTTAAGGGCTATCTTCATTTTTTTTCACAATGGCTTCAAAAGCGATTATacgc                                                                                                >  2:1176353/1‑149 (MQ=255)
                                                       tgCTCAAGTGCGTGTGCCATCGTGATTGAAAGCGGCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATccgccg                                                                                              <  1:123195/148‑1 (MQ=255)
                                                        gCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCcgccgc                                                                                             <  2:1471315/149‑1 (MQ=255)
                                                        gCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCAAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCcgccgc                                                                                             >  2:986509/1‑149 (MQ=255)
                                                          gcAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATACGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCTATTATCCGCCGCGa                                                                                           <  1:1167639/148‑1 (MQ=255)
                                                               tGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAgg                                                                                      <  2:784418/149‑1 (MQ=255)
                                                                    gtgCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGa                                                                                                        >  1:476005/1‑126 (MQ=255)
                                                                      gCCATCGAGAATGCAAGCGCCTTACTTTTGAGATCCGCGGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGAc                                                                               <  2:1010395/149‑1 (MQ=255)
                                                                        cATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATATCCCAGTCATAACTCTCAAAAACCACTATTGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACaa                                                                             >  2:1418818/1‑149 (MQ=255)
                                                                         aTCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCAGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAAt                                                                            >  1:390581/1‑149 (MQ=255)
                                                                         aTCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATATCCCAGTAATAACTATCAAAAACCAATATAGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAAt                                                                            >  2:1408222/1‑149 (MQ=255)
                                                                          tCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTTCCCAACCGATATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATa                                                                           <  1:820965/149‑1 (MQ=255)
                                                                          tCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATa                                                                           >  1:1046675/1‑149 (MQ=255)
                                                                          tCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATa                                                                           <  1:950200/149‑1 (MQ=255)
                                                                            gTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACTATCCCAAAACGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAgg                                                                         >  2:54338/1‑149 (MQ=255)
                                                                             tGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTTATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGa                                                                        <  1:1248583/149‑1 (MQ=255)
                                                                             ggAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGa                                                                        <  1:937065/148‑1 (MQ=255)
                                                                                  gCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGaaa                                                                   <  2:840425/149‑1 (MQ=255)
                                                                                    aaGCGCCTTTCTTTTGAGATCCGCTGTGTAGATTTCCCCGTGTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAgg                                                                 <  1:194229/149‑1 (MQ=255)
                                                                                         ccTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAg                                                            <  2:429097/149‑1 (MQ=255)
                                                                                          cTTACTTTTGAGATCCGCTGTGTAGATATCCCAGTAATAACTATCAAAAAAAAAAAATACTGATGGGCTATCTTCATTTTTTTGCACAATGGGTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGATAGGCaatgg                                                           >  2:143869/1‑146 (MQ=255)
                                                                                           ttACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAAt                                                                            <  2:138156/131‑1 (MQ=255)
                                                                                           ttACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTg                                                          >  2:1335120/1‑149 (MQ=255)
                                                                                             aCTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGt                                                        >  1:1516981/1‑149 (MQ=255)
                                                                                               ttttGAGATACGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTATTTTGCACAATGGCTTCAAAAGCTATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGAAAATTGGTCg                                                      >  1:435839/1‑149 (MQ=255)
                                                                                                tttGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGt                                                     >  1:610369/1‑149 (MQ=255)
                                                                                                  tGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGttt                                                   <  1:53630/149‑1 (MQ=255)
                                                                                                   gagaTCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTta                                                  >  1:840990/1‑149 (MQ=255)
                                                                                                   gagaTCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTta                                                  <  2:1027176/149‑1 (MQ=255)
                                                                                                      aTCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAg                                               <  1:1266274/149‑1 (MQ=255)
                                                                                                      aTCCGCTGTGTAGATTTCCCCGTCATTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCTCAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGGTTAAAg                                               >  1:277453/1‑149 (MQ=255)
                                                                                                       tCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGAGGTCGtt                                                    >  1:1415822/1‑143 (MQ=255)
                                                                                                        ccGCTGTGTAGATTTCCACGTCTTTACGGTGCCGAACCGCTATGGCTGATGGGCTATCTTGATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGtt                                             <  2:1100940/149‑1 (MQ=255)
                                                                                                         cGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTg                                            >  1:1494164/1‑149 (MQ=255)
                                                                                                             gtgtAGATTTCCCAGTCTTTACGCTGCCAAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGa                                        >  2:366141/1‑149 (MQ=255)
                                                                                                                tAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAg                                     >  1:822707/1‑149 (MQ=255)
                                                                                                                tAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAg                                     >  2:715236/1‑149 (MQ=255)
                                                                                                                tAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAg                                     <  1:1516128/149‑1 (MQ=255)
                                                                                                                tAGATTACGCCGTCTTTTCGCTGGCCAACCGCTATGGCTGATGGGCTATCTTCTTTTTTTTGCACAATGTCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGTAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAg                                     <  1:1176353/149‑1 (MQ=255)
                                                                                                                   aTTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGc                                  >  2:1313998/1‑149 (MQ=255)
                                                                                                                   aTTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGc                                  <  2:891668/149‑1 (MQ=255)
                                                                                                                    tttCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGca                                 <  1:1306973/149‑1 (MQ=255)
                                                                                                                    tttCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGca                                 >  1:1343868/1‑149 (MQ=255)
                                                                                                                        cccGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAgg                             >  1:398871/1‑149 (MQ=255)
                                                                                                                         ccGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGa                            >  1:1204228/1‑149 (MQ=255)
                                                                                                                          cGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGaa                           <  1:61762/149‑1 (MQ=255)
                                                                                                                            tCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATc                         <  1:1465687/149‑1 (MQ=255)
                                                                                                                                aagatagCCCATCAGCCATAGCGGTTGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAg                      <  2:28185/142‑1 (MQ=255)
                                                                                                                                aagatagCCCATCAGCCATAGCGGTTGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAg                      >  1:28185/7‑148 (MQ=255)
                                                                                                                                 aCGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAGtt                    >  2:1091611/1‑149 (MQ=255)
                                                                                                                                    tagCCCATCAGCCATAGCGGTTGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAGTTCtt                 >  1:1515669/3‑149 (MQ=255)
                                                                                                                                     tGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAGTTCttt                <  2:1046675/149‑1 (MQ=255)
                                                                                                                                         cAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAGTTCTTTATCa            <  2:424463/149‑1 (MQ=255)
                                                                                                                                         cAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGATTCAAAAGCGATTATCCGCCGTGATAAGTTATCGACAATATGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAGTTCTTTATCa            <  1:1278509/149‑1 (MQ=255)
                                                                                                                                           aCCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAGTTCTTTATCAGa          <  2:1578175/149‑1 (MQ=255)
                                                                                                                                           aCCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAGTTCTTTATCAGa          >  1:117187/1‑149 (MQ=255)
                                                                                                                                           aCCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGAACGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAGGTCTTTATCAGa          >  2:239296/1‑149 (MQ=255)
                                                                                                                                                 aTGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAGTTCTTTATCAGATCCGTc    >  2:1082829/1‑149 (MQ=255)
                                                                                                                                                 aTGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAGTTCTTTATCAGATCCGTc    <  1:1162155/149‑1 (MQ=255)
                                                                                                                                                   ggCTGATGGGCTATATTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGAGAATAGGATGAAAGGCAAAGTGGTCGTTTTTAGTTGGTGAAAGGGCACAGGAATCAAGTTCTTTATCAGATCCGTCAg  <  2:435839/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAGTTCTTTATCAGATCCGTCAG  >  NC_000913/1214720‑1215015

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