Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,888,034 A→T 7.1% intergenic (‑43/+27) rnd ← / ← fadD RNase D/fatty acyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,888,0340AT 7.1% 256.4 / 8.3 84intergenic (‑43/+27)rnd/fadDRNase D/fatty acyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base A (37/41);  new base T (2/4);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.52e-01

CGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCGT  >  NC_000913/1887890‑1888183
                                                                                                                                                 |                                                                                                                                                     
cgTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGc                                                                                                                                                  >  2:142296/1‑150 (MQ=255)
 gTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTATTCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCg                                                                                                                                                 <  2:321558/150‑1 (MQ=255)
  tACGAACAAATTCAGTATCCAGGGCTATCGCCGTAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGTTAATCCGGCGt                                                                                                                                                <  1:1086709/150‑1 (MQ=255)
   aCGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTc                                                                                                                                               <  1:503862/150‑1 (MQ=255)
     gAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCt                                                                                                                                              <  2:880724/150‑1 (MQ=255)
     gAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTg                                                                                                                                             >  2:1082901/1‑150 (MQ=255)
        cAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGTTAATCCGGCGttttttt                                                                                                                                          <  2:1665204/150‑7 (MQ=255)
        cAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTgacg                                                                                                                                          <  1:494118/150‑1 (MQ=255)
        cAAATTCAGTATCCAGGGCTATAGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGTTAATCCGGCGttttttt                                                                                                                                          <  1:68865/150‑7 (MQ=255)
         aaaTTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGTTAATCCGGCGtttttttt                                                                                                                                         <  2:1358467/150‑8 (MQ=255)
           aTTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGGCGGATTAACCGGCGTCTGACgact                                                                                                                                       <  2:507793/150‑1 (MQ=255)
              cAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACgactgac                                                                                                                                    <  2:1567791/150‑1 (MQ=255)
               aGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACgactgact                                                                                                                                   <  1:1263088/150‑1 (MQ=255)
               aGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACgactgact                                                                                                                                   <  1:1331142/150‑1 (MQ=255)
               aGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACgactgact                                                                                                                                   <  2:441092/150‑1 (MQ=255)
                  aTCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTaa                                                                                                                                <  2:1662524/150‑1 (MQ=255)
                     cAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGc                                                                                                                             <  2:930512/150‑1 (MQ=255)
                      aGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCt                                                                                                                            >  2:302620/1‑150 (MQ=255)
                        ggCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCa                                                                                                                          >  1:625667/1‑150 (MQ=255)
                         gCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAg                                                                                                                         <  1:610493/150‑1 (MQ=255)
                           tATCGCCGGAAAGGCACGGACGTCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGc                                                                                                                       <  1:895428/150‑1 (MQ=255)
                              cGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCttt                                                                                                                    >  2:1450510/1‑150 (MQ=255)
                               gCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTa                                                                                                                   <  1:1082901/150‑1 (MQ=255)
                               gCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGCCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTa                                                                                                                   <  2:657852/150‑1 (MQ=255)
                                  ggAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTg                                                                                                                <  2:1178148/150‑1 (MQ=255)
                                  ggAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTg                                                                                                                <  2:1557224/150‑1 (MQ=255)
                                  ggAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTg                                                                                                                <  1:796195/150‑1 (MQ=255)
                                      aGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCa                                                                                                            >  2:1550960/1‑150 (MQ=255)
                                       ggCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCAc                                                                                                           <  2:920811/150‑1 (MQ=255)
                                       ggCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCAc                                                                                                           <  1:1023255/150‑1 (MQ=255)
                                       ggCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACAGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCaa                                                                                                           >  2:1532900/1‑149 (MQ=255)
                                        gCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACg                                                                                                                              <  2:342208/130‑1 (MQ=255)
                                        gCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACg                                                                                                                              >  1:342208/1‑130 (MQ=255)
                                        gCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACt                                                                                                          <  1:302620/150‑1 (MQ=255)
                                                 gCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgcgc                                                                                                 >  1:800818/1‑150 (MQ=255)
                                                  cTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgcgcg                                                                                                >  2:937470/1‑150 (MQ=255)
                                                  cTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgcgcg                                                                                                <  1:1092794/150‑1 (MQ=255)
                                                   ttCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCgcgcgc                                                                                               >  2:1355953/1‑150 (MQ=255)
                                                       cacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCg                                                                                           <  2:449269/150‑1 (MQ=255)
                                                       cacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCg                                                                                           <  1:1505960/150‑1 (MQ=255)
                                                         caAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTc                                                                                         >  1:512801/1‑150 (MQ=255)
                                                          aaaGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCa                                                                                        <  2:415718/150‑1 (MQ=255)
                                                            agaagCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACg                                                                                      <  2:565909/150‑1 (MQ=255)
                                                                 ccAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCCGCGTCTTACGACTGACTTAACGCTCAGGCTATATTGTCCACTTTGCCGCGCGCTTCCTCACGTAAtt                                                                                 >  1:658848/1‑150 (MQ=255)
                                                                  cAGGGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATtc                                                                                <  2:711477/150‑1 (MQ=255)
                                                                  cAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTTCCGCGCGCTTCGTCACGTAAtta                                                                                >  2:1310939/1‑149 (MQ=255)
                                                                  cAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATtc                                                                                >  1:1486451/1‑150 (MQ=255)
                                                                     cgcATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCtcg                                                                             <  1:1345813/150‑1 (MQ=255)
                                                                       cATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCtcgtc                                                                           >  1:1326440/1‑150 (MQ=255)
                                                                        aTCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCtcgtcg                                                                          >  1:1530857/1‑150 (MQ=255)
                                                                           gTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACg                                                                                      >  1:316784/1‑135 (MQ=255)
                                                                           gTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACg                                                                                      <  2:316784/135‑1 (MQ=255)
                                                                             ccGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCaaaa                                                                     <  2:1201018/150‑1 (MQ=255)
                                                                              cGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAAt                                                                    >  1:1603078/1‑150 (MQ=255)
                                                                                  gTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAAttttt                                                                >  1:1603594/1‑150 (MQ=255)
                                                                                   tAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTc                                                               <  1:1170819/150‑1 (MQ=255)
                                                                                                aTTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGAttt                                                  >  1:824239/1‑150 (MQ=255)
                                                                                                 ttCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTc                                                 >  1:861806/1‑150 (MQ=255)
                                                                                                  tCAAAATGTTTTCTCTTTAGTGGGCTTAAAAAAAAAAGACGGATTAACCCGAGACAGACGACTGAATTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTATCGTCGCAAAATTTTTCCGACGTTAGATTTCg                                                >  2:1063109/1‑150 (MQ=255)
                                                                                                        tGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACt                                          <  2:948749/150‑1 (MQ=255)
                                                                                                         gTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtc                                         >  1:1364023/1‑150 (MQ=255)
                                                                                                           tttCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtcat                                       >  2:711574/1‑150 (MQ=255)
                                                                                                           tttCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtcat                                       >  2:223397/1‑150 (MQ=255)
                                                                                                              ctctTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCAc                                    <  1:1550960/150‑1 (MQ=255)
                                                                                                               tctTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACg                                   >  1:103872/1‑150 (MQ=255)
                                                                                                                ctTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGa                                  >  1:1226710/1‑150 (MQ=255)
                                                                                                                ctTTAGTGGGCATCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCGTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGa                                  >  2:867612/1‑150 (MQ=255)
                                                                                                                    aGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACt                              <  2:439298/150‑1 (MQ=255)
                                                                                                                         gCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTccacc                         <  2:824239/150‑1 (MQ=255)
                                                                                                                             cAAAAAAAATGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCt                     >  2:1486842/1‑150 (MQ=255)
                                                                                                                             cAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCt                     >  1:119246/1‑150 (MQ=255)
                                                                                                                               aaaaaaaaCGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCtt                    >  1:1484324/1‑150 (MQ=255)
                                                                                                                               aaaaaaaaCGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCtt                    >  2:683009/1‑150 (MQ=255)
                                                                                                                                cgtcagccGCCGGTTAATCCGGCGTCTGACGACTGCCTTACCGCTCAGGCTTTATTGTCCCCTTTGCCGCGCGATTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTc                   >  2:161810/8‑150 (MQ=255)
                                                                                                                                aaaaaaaCGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTc                   >  2:1200566/1‑150 (MQ=255)
                                                                                                                                  aaaaaCGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCgg                 <  1:400564/150‑1 (MQ=255)
                                                                                                                                    agaCGCCGGTTAATCCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTa               >  2:1354521/3‑150 (MQ=255)
                                                                                                                                    aaaCGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCGAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTa               <  1:687644/150‑1 (MQ=255)
                                                                                                                                     aaCGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTAc              >  2:1624656/1‑150 (MQ=255)
                                                                                                                                        gCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACttt           >  1:1288879/1‑150 (MQ=255)
                                                                                                                                         ccGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGGTTCGTAACGTAATTCTCGTCTCAAAATTTTTCCGACGTTAGATTTGGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTg          <  2:81405/150‑1 (MQ=255)
                                                                                                                                           ggATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTa        <  1:510168/150‑1 (MQ=255)
                                                                                                                                            gATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTAt       >  2:1086577/1‑150 (MQ=255)
                                                                                                                                               tAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATccc    <  1:1200566/150‑1 (MQ=255)
                                                                                                                                                aaCCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCg   >  2:26373/1‑150 (MQ=255)
                                                                                                                                                 aCCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCGt  >  2:1492845/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCGT  >  NC_000913/1887890‑1888183

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: