Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,214,497 T→A 5.3% I327F (ATC→TTC)  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,4970TA 5.3% 246.6 / 2.1 85I327F (ATC→TTC) 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/44);  new base A (2/3);  total (38/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCA  >  NC_000913/1214348‑1214645
                                                                                                                                                     |                                                                                                                                                    
tGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGAt                                                                                                                                                      <  1:20769/150‑1 (MQ=255)
    gATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTa                                                                                                                                                  <  1:716632/150‑1 (MQ=255)
        cAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTAGAATTTCAAGTGATGTACAGCTTTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTc                                                                                                                                              >  1:189069/1‑150 (MQ=255)
            ccACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTc                                                                                                                                          <  2:837845/150‑1 (MQ=255)
              acaTCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATTATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAg                                                                                                                                        <  2:364191/150‑1 (MQ=255)
              acaTCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAg                                                                                                                                        >  2:755024/1‑150 (MQ=255)
                aTCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACAGAATTAAAATCAGTCAgga                                                                                                                                      <  1:1115064/150‑2 (MQ=255)
                aTCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACAGAATTAAAATCAGTCAgga                                                                                                                                      <  1:1120255/150‑2 (MQ=255)
                    aCTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGa                                                                                                                                  <  1:1019598/150‑1 (MQ=255)
                     ctctTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACAGAATTAAAATCAGTCAggaattga                                                                                                                                 <  2:888208/150‑7 (MQ=255)
                      tcttctGGTGTTGCCACGCCCATAGCACTGACTTGAATTTGAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAg                                                                                                                                <  2:113119/150‑1 (MQ=255)
                       cttctGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGc                                                                                                                               <  1:863033/150‑1 (MQ=255)
                          ctGGTGTTGCCTCGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATATTTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTg                                                                                                                            <  1:13828/150‑1 (MQ=255)
                          atGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGACTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTg                                                                                                                            <  1:124770/149‑1 (MQ=255)
                           tGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGa                                                                                                                           >  1:139399/1‑150 (MQ=255)
                            ggTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGAc                                                                                                                          >  2:478296/1‑150 (MQ=255)
                              tgtTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAg                                                                                                                        >  1:1055292/1‑150 (MQ=255)
                               gtTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAgg                                                                                                                       <  1:389498/150‑1 (MQ=255)
                                ttGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGa                                                                                                                      >  1:956236/1‑150 (MQ=255)
                                 tGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAg                                                                                                                     <  1:867521/150‑1 (MQ=255)
                                  gCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGt                                                                                                                    >  1:1211096/1‑150 (MQ=255)
                                   ccACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGACCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTa                                                                                                                   <  1:675056/150‑1 (MQ=255)
                                    cACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTaa                                                                                                                  >  2:1074808/1‑150 (MQ=255)
                                     actcacATACCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATTATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAAc                                                                                                                 <  2:181675/145‑1 (MQ=255)
                                     aCGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAAc                                                                                                                 <  2:13336/150‑1 (MQ=255)
                                      cGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAAcn                                                                                                                <  2:751298/150‑2 (MQ=255)
                                      cGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACa                                                                                                                <  2:746068/150‑1 (MQ=255)
                                         ccATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAAc                                                                                                             >  2:677648/1‑150 (MQ=255)
                                            tAGCACTGACTTGAATTTCCAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAg                                                                                                          <  2:860773/150‑1 (MQ=255)
                                             aGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAAGCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCacc                                                                                                         <  2:370213/150‑3 (MQ=255)
                                              gCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGATCATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACCAACCAGCa                                                                                                        >  2:1301905/1‑150 (MQ=255)
                                                 cTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCaaaa                                                                                                     >  2:1055982/1‑150 (MQ=255)
                                                   gACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTCATTTTAATTCTGTCAGGCTGGAAGCGTGCCAGGAGTAACAAACCAGCAAAAcc                                                                                                   >  2:1276215/1‑150 (MQ=255)
                                                       tGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCATGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTgcg                                                                                               >  1:811783/1‑150 (MQ=255)
                                                       tGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTgcg                                                                                               >  2:275174/1‑150 (MQ=255)
                                                       tGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTgcg                                                                                               <  2:685233/150‑1 (MQ=255)
                                                         aaTTTCATGTGATGAACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGcc                                                                                             <  2:267812/150‑1 (MQ=255)
                                                              cAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACGTGCGCCAAAAt                                                                                        >  1:675806/1‑150 (MQ=255)
                                                                tgTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGAGTCCCACTTATATGAAAATTGGTAATCAAATACTGAATGATTTTAATTCTGTCAGGCTGTAAGAGTGACAGGAGTAACAAACCAGCAAAAAATGAGCCAACATGa                                                                                      >  1:908141/2‑150 (MQ=255)
                                                                attGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGACTTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGa                                                                                      <  1:935490/148‑1 (MQ=255)
                                                                      gTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAAtt                                                                                >  1:422858/1‑150 (MQ=255)
                                                                       tACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGTAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTg                                                                               >  1:114168/1‑150 (MQ=255)
                                                                       tACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTg                                                                               >  1:874374/1‑150 (MQ=255)
                                                                        aCAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGc                                                                              <  2:189069/150‑1 (MQ=255)
                                                                          aGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTa                                                                            >  1:175022/1‑150 (MQ=255)
                                                                          aGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTa                                                                            <  1:1043095/150‑1 (MQ=255)
                                                                                ttATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAAAAAACCAGCAAACCCTGCGCCAAAATGATCAATTGCTACACTGa                                                                      >  1:852746/1‑150 (MQ=255)
                                                                                   tGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATat                                                                   <  1:226313/150‑2 (MQ=255)
                                                                                    gATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTATTCAATACCTGACTGAATTAAATTCTGTCAGGCTGCAAGCGTGACAGTAGTAACAATCCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATAcc                                                                  >  2:1227781/1‑150 (MQ=255)
                                                                                         ccTGAATTATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACcagcag                                                             >  1:540998/1‑150 (MQ=255)
                                                                                         ccTGAATAATTGCCTGCCGTGGCCAACTTTTATGAACATTGGTAATCAATTCCTGAATGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATAAATTGCTAAAATGATACcagcag                                                             >  2:238638/1‑150 (MQ=255)
                                                                                          cTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCagcagc                                                            <  2:883359/150‑1 (MQ=255)
                                                                                           tGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGcc                                                           <  2:667521/150‑1 (MQ=255)
                                                                                               taatTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTa                                                       >  2:157878/1‑150 (MQ=255)
                                                                                                aatTGCCTTCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAg                                                      <  1:1144663/150‑1 (MQ=255)
                                                                                                  ttGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCg                                                    <  1:493169/150‑1 (MQ=255)
                                                                                                   tgcctgccGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGGGa                                                   >  1:356620/1‑150 (MQ=255)
                                                                                                    gcctgccGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGAt                                                  <  1:554234/150‑1 (MQ=255)
                                                                                                     cctgccGTGGCCCACTTTTATGAACATTGGTAATCTATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAACCCCGCTAGCCCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCGTTGAGCgctg                                                 >  2:1204701/1‑147 (MQ=255)
                                                                                                     cctgccGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGAtt                                                 >  2:1101515/1‑150 (MQ=255)
                                                                                                     cctgccGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGAtt                                                 <  1:1263379/150‑1 (MQ=255)
                                                                                                          cGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATc                                            <  1:397825/150‑1 (MQ=255)
                                                                                                           gTGGCCCACTTTTATGAACATTTGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCg                                           <  2:643731/150‑1 (MQ=255)
                                                                                                           gTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCg                                           <  1:1265111/150‑1 (MQ=255)
                                                                                                             ggCCCACTTTTATGAACATTTGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAAACTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCt                                         >  2:503584/1‑150 (MQ=255)
                                                                                                             ggCCCACTTTTATGAACATTGGTTATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCt                                         >  2:355218/1‑150 (MQ=255)
                                                                                                               cccACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTc                                       <  2:283662/150‑1 (MQ=255)
                                                                                                               cccACTTTTATGAACATTGGTAATCAATTCCTGACTCATTTTAATTCTGTCAGGCTGGAAGATTGACAGTAGTAACAAACCAGCAAAACATGAGCATAAATGAAAAATTGATACACTGATACCAGCAGCCTTTATCGATTTAATCgaatc                                       >  1:153809/1‑146 (MQ=255)
                                                                                                                       tATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACt                               <  1:677648/150‑1 (MQ=255)
                                                                                                                       tATGAACATTGGTAATCAATTACTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAACCCAGCATAACCTGCGCCAAAATGATCTATTGCTTCACTGATACCAGCAGCCTTTTGCGTTTTAATCGCTTCGGCAAACt                               <  2:638492/150‑1 (MQ=255)
                                                                                                                            aCATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCTAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACtcatca                          >  1:625083/1‑150 (MQ=255)
                                                                                                                                tGGTAATCAAATCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGATACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCATACTCATCAAAcc                      >  1:265063/1‑150 (MQ=255)
                                                                                                                                   tAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCgag                   <  2:114168/150‑1 (MQ=255)
                                                                                                                                     aTCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGat                 <  2:608461/150‑1 (MQ=255)
                                                                                                                                      taaaTTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGata                <  1:192930/148‑1 (MQ=255)
                                                                                                                                      tCAATTCCTGAATTATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGata                <  2:83137/150‑1 (MQ=255)
                                                                                                                                         cagccTGACAGAATTAAAATCAGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGCTCAATTGCTACACTGATACCAGCAGCCTTTAGAGATTTAATCGCTTCGGCAAACTCATCTACCCGAGATATGa             >  1:179577/4‑150 (MQ=255)
                                                                                                                                         aTTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGa             <  1:1045478/150‑1 (MQ=255)
                                                                                                                                           tCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACt           <  1:580577/150‑1 (MQ=255)
                                                                                                                                            ccTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACtt          >  1:1236890/1‑150 (MQ=255)
                                                                                                                                               gactgaTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCAc       <  2:385059/150‑1 (MQ=255)
                                                                                                                                               gactgaTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCAc       >  1:1022913/1‑150 (MQ=255)
                                                                                                                                                actgaTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACt      <  2:151269/150‑1 (MQ=255)
                                                                                                                                                actgaTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACt      >  1:1093297/1‑150 (MQ=255)
                                                                                                                                                 ctgaTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACtt     <  1:1002793/150‑1 (MQ=255)
                                                                                                                                                    aTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTAATCAAACCGAGATATGACTTCACTTTCa  >  2:589783/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCA  >  NC_000913/1214348‑1214645

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: