Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,214,515 G→T 5.0% Q321K (CAG→AAG)  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,5150GT 5.0% 357.8 / 7.6 120Q321K (CAG→AAG) 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 (64/50);  new base T (4/2);  total (68/52)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.29e-01

Only 100 of 124 total aligned reads displayed.

CCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATG  >  NC_000913/1214366‑1214660
                                                                                                                                                     |                                                                                                                                                  
ccACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCATTTCCTGACTGATTTTAATTCTGTCAGGCTg                                                                                                                                                    <  2:144411/150‑1 (MQ=255)
 cACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAATATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTgg                                                                                                                                                   >  1:1494052/1‑150 (MQ=255)
 cACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTgg                                                                                                                                                   >  1:402709/1‑150 (MQ=255)
  aCTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGa                                                                                                                                                  <  1:1520889/150‑1 (MQ=255)
   ctctTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGaa                                                                                                                                                 >  2:1242916/1‑150 (MQ=255)
     cttctGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGc                                                                                                                                               >  1:560076/1‑150 (MQ=255)
       tctGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGt                                                                                                                                             >  2:30226/1‑150 (MQ=255)
        ctGTTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGAGTTTAATTCTGTCAGGCTGGAAGCGTg                                                                                                                                            <  2:719378/150‑1 (MQ=255)
        ctGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTg                                                                                                                                            >  2:124530/1‑150 (MQ=255)
        ctGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTg                                                                                                                                            <  1:520639/150‑1 (MQ=255)
        ctGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTg                                                                                                                                            >  2:578281/1‑150 (MQ=255)
        ctGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTg                                                                                                                                            <  1:1013868/150‑1 (MQ=255)
          ggTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGAc                                                                                                                                          >  2:1579721/1‑150 (MQ=255)
           gtgtTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACa                                                                                                                                         <  1:794172/150‑1 (MQ=255)
            tgtTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTTCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAg                                                                                                                                        <  2:211423/150‑1 (MQ=255)
            tgtTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAg                                                                                                                                        >  2:1405030/1‑150 (MQ=255)
             gtTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAgg                                                                                                                                       >  2:950225/1‑150 (MQ=255)
             gtTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAgg                                                                                                                                       <  2:327542/150‑1 (MQ=255)
               tGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCACGCTTCCAGCCTGACAgaat                                                                                                                                     <  2:634544/150‑4 (MQ=255)
                gCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGt                                                                                                                                    >  1:1469293/1‑150 (MQ=255)
                  cACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTaa                                                                                                                                  >  2:1604547/1‑150 (MQ=255)
                  cACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTaa                                                                                                                                  <  1:1719126/150‑1 (MQ=255)
                  cACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCATGACTGATTTTAATTCTGTCAGGCTGGAAGAGTGACAGGAGTaa                                                                                                                                  >  1:862926/1‑150 (MQ=255)
                   aCGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCACGCTTCCAGCCTGACAGAATTaaa                                                                                                                                 <  1:638465/150‑2 (MQ=255)
                    cGCCCATAGAACTGACTTGAATTTCAAGTTATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACa                                                                                                                                <  2:1153248/150‑1 (MQ=255)
                       ccATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAAc                                                                                                                             >  1:807092/1‑150 (MQ=255)
                        cATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTTACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAAcc                                                                                                                            <  1:249673/150‑1 (MQ=255)
                          tagaaatGACTTGAATTTCAAGTGATGTACATCATTTTATGATCTCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTACTGAATGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAg                                                                                                                          <  2:444560/144‑1 (MQ=255)
                          tAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAg                                                                                                                          <  2:1341389/150‑1 (MQ=255)
                               cTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCaaaa                                                                                                                     >  1:1203523/1‑150 (MQ=255)
                                 gACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAAcc                                                                                                                   <  1:1644492/150‑1 (MQ=255)
                                 gACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAAcc                                                                                                                   <  1:1579022/150‑1 (MQ=255)
                                        aTTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCa                                                                                                            >  1:1300696/1‑150 (MQ=255)
                                           tCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCaaaa                                                                                                         >  1:1680875/1‑150 (MQ=255)
                                                 gATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCa                                                                                                   <  1:1268640/150‑1 (MQ=255)
                                                   tGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAAt                                                                                                 <  1:574507/150‑1 (MQ=255)
                                                      aCAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGc                                                                                              <  1:1125274/150‑1 (MQ=255)
                                                      aCAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGc                                                                                              >  2:521742/1‑150 (MQ=255)
                                                        aGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTa                                                                                            >  2:748114/1‑150 (MQ=255)
                                                         gCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTac                                                                                           <  1:1405030/150‑1 (MQ=255)
                                                            ttttATGATCGCCTGAATAATTGCCTGCCGTGGCCTACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACt                                                                                        >  2:559391/1‑150 (MQ=255)
                                                            ttttATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACt                                                                                        <  1:1584637/150‑1 (MQ=255)
                                                                  gATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGACGCGTGACAGGAGTAACATACCAGCAAATCCTGCGCCAAAATGATCAATTGCTACACTGATAcc                                                                                  >  1:272616/1‑150 (MQ=255)
                                                                  gATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAATAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATAcc                                                                                  <  2:1223922/150‑1 (MQ=255)
                                                                  gATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATAcc                                                                                  >  2:1709042/1‑150 (MQ=255)
                                                                  gATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATAcc                                                                                  >  1:1228442/1‑150 (MQ=255)
                                                                     cGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAGAACCTGCGCCAAAATGATCAATTGCTACACTGATACcagc                                                                               >  2:1677822/1‑150 (MQ=255)
                                                                      gCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGTAGTAACAAACCAGCAAAACCTGCGCCAAAATTATCAATTTCTACACTGATACaagcc                                                                              >  2:758594/1‑149 (MQ=255)
                                                                       ccTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACcagcag                                                                             <  2:1715926/150‑1 (MQ=255)
                                                                          gAATAATTGCCTGCCGTGTCCCACTTTTATGAACATTGGTAATCAAATCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑t                                                                         <  2:944354/150‑1 (MQ=255)
                                                                           aataatTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑tt                                                                        <  2:372873/150‑1 (MQ=255)
                                                                               atTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAGCCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑ATTAGc                                                                    >  2:471390/1‑150 (MQ=255)
                                                                               atTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATATTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGACAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGc                                                                    >  2:1227981/1‑150 (MQ=255)
                                                                                 tgcctgccGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGa                                                                  >  2:1452057/1‑150 (MQ=255)
                                                                                  gcctgccGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGAt                                                                 <  2:1114882/150‑1 (MQ=255)
                                                                                   cctgccGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGAtt                                                                <  2:1703439/150‑1 (MQ=255)
                                                                                       ccGTGCCCCACTTTTATGAACAATGGTTATCAATACCTCACAGATTTTAAATCAGAAAGGCAGGAAGCGAGACAGGACAAACAAACCAGCAAAACCTGCCCCAAAAAGATCAATTCCTACACTGATACCAGCAGCC‑TATACCGATTTaaa                                                            >  2:960259/1‑149 (MQ=255)
                                                                                           ggCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCt                                                        <  2:446996/150‑1 (MQ=255)
                                                                                             cccACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTc                                                      <  1:1212656/150‑1 (MQ=255)
                                                                                             cccACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGAGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTTAGCGATTTAATCGCtt                                                       >  1:58029/1‑150 (MQ=255)
                                                                                              ccACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCg                                                     >  2:790850/1‑150 (MQ=255)
                                                                                                 cTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCa                                                  >  1:149114/1‑150 (MQ=255)
                                                                                                      aTGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACtc                                             >  2:1183955/1‑150 (MQ=255)
                                                                                                         aaCATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACtcatc                                          >  1:1324947/1‑150 (MQ=255)
                                                                                                          accttGGTAATCAATTCCTGACTGAATTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAACCCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACtcatca                                         <  2:1187433/147‑1 (MQ=255)
                                                                                                          aCATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACtcatca                                         <  1:559391/150‑1 (MQ=255)
                                                                                                          aCATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACtcatca                                         <  2:752895/150‑1 (MQ=255)
                                                                                                          aCATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACATGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACtcatca                                         >  1:224746/1‑150 (MQ=255)
                                                                                                           cATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCaa                                        >  1:1592867/1‑150 (MQ=255)
                                                                                                             ttGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACATGAGTAACAAACCAGCAAAACCTGCGCCAAAATGAACAATTGCTACACTGATACCAGCAGCC‑TTTATCGATTTAATCGCTTCTGCAAACTCATCAAAc                                      >  2:1130432/1‑150 (MQ=255)
                                                                                                             ttGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAAc                                      >  2:437273/1‑150 (MQ=255)
                                                                                                               ggTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGTGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCg                                    >  1:915615/1‑150 (MQ=255)
                                                                                                               ggTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCg                                    >  1:739388/1‑150 (MQ=255)
                                                                                                               ggTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCg                                    <  1:1131181/150‑1 (MQ=255)
                                                                                                               ggTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCg                                    <  2:1029037/150‑1 (MQ=255)
                                                                                                                  aaTCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCgaaa                                 <  1:664516/150‑3 (MQ=255)
                                                                                                                   aTCAATTCCTGACTGATTTTAATTCTGTCAGTCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAAATGCTACACAGATACCAGCAGCC‑TTTAGCGATTTATTCGCTTCGGCAAACTCATCAAACCGAGat                                <  1:30226/150‑1 (MQ=255)
                                                                                                                   aTCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGat                                >  1:532921/1‑150 (MQ=255)
                                                                                                                   aTCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGat                                <  1:1579721/150‑1 (MQ=255)
                                                                                                                     cAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGatat                              >  1:1386397/1‑150 (MQ=255)
                                                                                                                     cAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGatat                              >  1:13899/1‑150 (MQ=255)
                                                                                                                        ttCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGAc                           <  2:318392/150‑1 (MQ=255)
                                                                                                                           ctgactgaTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTc                        >  2:380426/1‑150 (MQ=255)
                                                                                                                            tgacagaATTAAAATCAGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCa                       >  2:637964/1‑150 (MQ=255)
                                                                                                                             gactgaTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCAc                      >  1:601918/1‑150 (MQ=255)
                                                                                                                              actgaTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACt                     >  2:1689999/1‑150 (MQ=255)
                                                                                                                              actgaTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGATTCGGCAAACTCATCAAACCGAGATATGACTTCACt                     >  2:1683438/1‑150 (MQ=255)
                                                                                                                                 gaTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAGAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTc                  <  1:1677822/150‑1 (MQ=255)
                                                                                                                                  aTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCa                 >  2:617212/1‑150 (MQ=255)
                                                                                                                                   ttttAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAg                >  2:1385844/1‑150 (MQ=255)
                                                                                                                                    tttAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGt               <  2:548605/150‑1 (MQ=255)
                                                                                                                                     ttAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTa              <  1:217999/150‑1 (MQ=255)
                                                                                                                                      tAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTaa             <  2:1680875/150‑1 (MQ=255)
                                                                                                                                      tAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTaa             >  1:388932/1‑150 (MQ=255)
                                                                                                                                        ctccTGTCACGCTTCCAGCCTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAAt           >  1:433363/4‑150 (MQ=255)
                                                                                                                                        ctccTGTCACGCTTCCAGCCTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACAATGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACAGAGATATGACTTCACTTTCAGTAAAt           >  1:471535/4‑150 (MQ=255)
                                                                                                                                          ccTGTCACGCTTCCAGCCTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATTCCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAAtta         >  1:1286674/2‑149 (MQ=255)
                                                                                                                                           cTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCa        >  2:937903/1‑150 (MQ=255)
                                                                                                                                               cAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACga    >  2:1405385/1‑150 (MQ=255)
                                                                                                                                                 ggCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACgatg  <  1:937903/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCC‑TTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATG  >  NC_000913/1214366‑1214660

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