Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,214,861 C→A 5.3% A205A (GCG→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,8610CA 5.3% 352.5 / 4.7 120A205A (GCG→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 C (51/61);  new base A (5/3);  total (56/64)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

TCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAGTTCTTTATCAGA  >  NC_000913/1214718‑1215007
                                                                                                                                               |                                                                                                                                                  
gcAGGTTCGTTAACCAGGGTCATAGGGAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTAt                                                                                                                                               <  1:1012953/148‑1 (MQ=255)
  aGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATgg                                                                                                                                             <  1:581870/149‑1 (MQ=255)
  aGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATgg                                                                                                                                             <  1:506777/149‑1 (MQ=255)
     ttCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTg                                                                                                                                          >  1:178296/1‑149 (MQ=255)
        gTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCATTACGCTGCCCAACCGCTATGCCTGATg                                                                                                                                       >  2:194229/1‑149 (MQ=255)
         ttAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATgg                                                                                                                                      >  2:738186/1‑149 (MQ=255)
          tAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATggg                                                                                                                                     <  1:461120/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)
                   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)
                                              ttttATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCg                                                                                                 >  1:840425/1‑149 (MQ=255)
                                               tttATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGGCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGa                                                                                                <  2:777933/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)
                                                                gTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCAAAATGGCTTAAAAAGATATTATCAGACGCGATaat                                                                               >  2:1278509/1‑148 (MQ=255)
                                                                 tGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAgg                                                                              <  2:784418/149‑1 (MQ=255)
                                                                      gtgCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGa                                                                                                >  1:476005/1‑126 (MQ=255)
                                                                      gtgCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGa                                                                                                <  2:476005/126‑1 (MQ=255)
                                                                       tgCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGa                                                                        <  1:738186/149‑1 (MQ=255)
                                                                       ggCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGGAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATCGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGa                                                                        <  2:202196/148‑1 (MQ=255)
                                                                        gCCATCGAGAATGCAAGCGCCTTACTTTTGAGATCCGCGGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGAc                                                                       <  2:1010395/149‑1 (MQ=255)
                                                                          cATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATATCCCAGTCATAACTCTCAAAAACCACTATTGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACaa                                                                     >  2:1418818/1‑149 (MQ=255)
                                                                           aTCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCAGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAAt                                                                    >  1:390581/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)
                                                                              gTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATATCCCAGTAATAACTATCACAAAAAAAAATAGATGATGGGCTATCTTCATTTTTTTGCACAATTGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAgg                                                                 >  2:746759/1‑149 (MQ=255)
                                                                               tGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGa                                                                <  1:1242472/149‑1 (MQ=255)
                                                                               ggAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGa                                                                <  1:937065/148‑1 (MQ=255)
                                                                                  aTGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTAACTCTGCCAAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGa                                                             >  2:1083134/1‑149 (MQ=255)
                                                                                      aaGCGCCTTTCTTTTGAGATCCGCTGTGTAGATTTCCCCGTGTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAgg                                                         <  1:194229/149‑1 (MQ=255)
                                                                                          gcCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCaaa                                                     >  1:1578175/1‑149 (MQ=255)
                                                                                           ccTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAg                                                    <  2:429097/149‑1 (MQ=255)
                                                                                            cTTACTTTTGAGATCCGCTGTGTAGATATCCCAGTAATAACTATCAAAAAAAAAAAATACTGATGGGCTATCTTCATTTTTTTGCACAATGGGTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGATAGGCaatgg                                                   >  2:143869/1‑146 (MQ=255)
                                                                                             ttACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAAt                                                                    <  2:138156/131‑1 (MQ=255)
                                                                                             ttACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAAt                                                                    >  1:138156/1‑131 (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)
                                                                                                   ttGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGtt                                            <  1:324313/149‑1 (MQ=255)
                                                                                                    tGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGttt                                           <  1:53630/149‑1 (MQ=255)
                                                                                                     gagaTCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTta                                          >  1:840990/1‑149 (MQ=255)
                                                                                                        aTCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAg                                       <  1:1266274/149‑1 (MQ=255)
                                                                                                        aTCCGCTGTGTAGATTTCCCCGTCATTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCTCAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGGTTAAAg                                       >  1:277453/1‑149 (MQ=255)
                                                                                                         tCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGAGGTCGtt                                            <  2:1415822/143‑1 (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)
                                                                                                                 gtAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGaaa                              <  1:1473636/149‑1 (MQ=255)
                                                                                                                 gtAGATTTCCCCGTCTTTACGCAGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGaaa                              <  2:390581/149‑1 (MQ=255)
                                                                                                                  tAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAg                             >  1:822707/1‑149 (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)
                                                                                                                          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)
                                                                                                                                 ttACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAg              >  1:957785/1‑149 (MQ=255)
                                                                                                                                  aagatagCCCATCAGCCATAGCGGTTGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAg              <  2:28185/142‑1 (MQ=255)
                                                                                                                                   aCGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAGtt            >  2:1091611/1‑149 (MQ=255)
                                                                                                                                    cGCTGCCCAACCGCTATGGCTGCTGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAGtta           >  2:581922/1‑148 (MQ=255)
                                                                                                                                      tagCCCATCAGCCATAGCGGTTGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAGTTCtt         >  1:1515669/3‑149 (MQ=255)
                                                                                                                                       tGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAGTTCttt        <  2:1046675/149‑1 (MQ=255)
                                                                                                                                           cAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGATTCAAAAGCGATTATCCGCCGTGATAAGTTATCGACAATATGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAGTTCTTTATCa    <  1:1278509/149‑1 (MQ=255)
                                                                                                                                             aCCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAGTTCTTTATCAGa  >  2:853840/1‑149 (MQ=255)
                                                                                                                                             aCCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAGTTCTTTATCAGa  >  1:117187/1‑149 (MQ=255)
                                                                                                                                             aCCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAGTTCTTTATCAGa  <  2:1578175/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGATAAGGGATCGACAATAGGATGAAAGGCAAAGTGGTCGTTTATAGTTGGTGAAAGGGCACAGGAATCAAGTTCTTTATCAGA  >  NC_000913/1214718‑1215007

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