Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,231,259 G→A 7.7% intergenic (+20/‑406) ygjK → / → fadH glycoside hydrolase/2,4‑dienoyl‑CoA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,231,2590GA 7.7% 361.4 / 13.7 116intergenic (+20/‑406)ygjK/fadHglycoside hydrolase/2,4‑dienoyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base G (62/45);  new base A (4/5);  total (66/50)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 123 total aligned reads displayed.

TTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGAT  >  NC_000913/3231111‑3231408
                                                                                                                                                    |                                                                                                                                                     
      gccgATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTa                                                                                                                                                >  1:1013252/1‑150 (MQ=255)
      gccgATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTa                                                                                                                                                <  1:1488794/150‑1 (MQ=255)
         gATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAAc                                                                                                                                             <  1:1004017/150‑1 (MQ=255)
         gATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAAc                                                                                                                                             >  1:736529/1‑150 (MQ=255)
             gACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCt                                                                                                                                         >  2:15620/1‑150 (MQ=255)
              aCCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCtt                                                                                                                                        >  2:931503/1‑150 (MQ=255)
               ccAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTa                                                                                                                                       <  2:746995/150‑1 (MQ=255)
                 aaTTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATc                                                                                                                                     <  1:406585/150‑1 (MQ=255)
                  aTTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                    >  1:675931/1‑150 (MQ=255)
                  aTTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                    >  1:1537815/1‑150 (MQ=255)
                    tCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGc                                                                                                                                  >  2:1110753/1‑150 (MQ=255)
                     cAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGcc                                                                                                                                 <  2:114993/150‑1 (MQ=255)
                      aGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGccc                                                                                                                                <  2:813537/150‑1 (MQ=255)
                      aGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGccc                                                                                                                                >  1:1261527/1‑150 (MQ=255)
                       ggAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATAAGGCGTTTACGCCGCATCCGccag                                                                                                                               <  2:149327/150‑3 (MQ=255)
                        gAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTAATCCGCCCTa                                                                                                                              >  2:94700/1‑150 (MQ=255)
                         aaaaTTACACCCCGCTGACCGGCGCACAGCAAGGCGCACCAAATTTCTCCTCCCGCGCCGCGCATTTGTATATGTTCTATAACGATTTTTTCCGTAAGCAAAAATTAGCAACTGCTGGCGGAGGCGGCGTAAACGCCTTATCCGCCCTAc                                                                                                                             <  2:1234031/150‑1 (MQ=255)
                            aTTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACAtg                                                                                                                          <  1:867360/150‑1 (MQ=255)
                            aTTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACAtg                                                                                                                          <  1:1110753/150‑1 (MQ=255)
                             ttACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACAtgt                                                                                                                         <  2:736529/150‑1 (MQ=255)
                             ttACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACAtgt                                                                                                                         <  1:512348/150‑1 (MQ=255)
                                  aCCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAAt                                                                                                                    >  1:116578/1‑150 (MQ=255)
                                  aCCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAAt                                                                                                                    >  1:1114231/1‑150 (MQ=255)
                                  aCCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAAt                                                                                                                    <  1:184149/150‑1 (MQ=255)
                                    ccGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAAtcc                                                                                                                  <  1:1158935/150‑1 (MQ=255)
                                           cAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGc                                                                                                           >  2:16948/1‑150 (MQ=255)
                                           cAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATAAGGCGTTTACGCCGCATCCGccagca                                                                                                                             >  1:824885/1‑128 (MQ=38)
                                           cAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATAAGGCGTTTACGCCGCATCCGccagca                                                                                                                             <  2:824885/132‑5 (MQ=38)
                                             ggCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCg                                                                                                         <  1:48052/150‑1 (MQ=255)
                                             ggCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCg                                                                                                         >  2:762181/1‑150 (MQ=255)
                                             ggCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCg                                                                                                         >  1:1351060/1‑150 (MQ=255)
                                             ggCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCg                                                                                                         >  2:847025/1‑150 (MQ=255)
                                             ggCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTCGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCg                                                                                                         >  1:548043/1‑150 (MQ=255)
                                              gcgcACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCgg                                                                                                        <  2:199834/150‑1 (MQ=255)
                                                gcACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGAt                                                                                                      >  1:1356121/1‑150 (MQ=255)
                                                   caGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATaaa                                                                                                   >  1:1511118/1‑150 (MQ=255)
                                                     gCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAAc                                                                                                 >  1:181033/1‑150 (MQ=255)
                                                         ggCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATGCAAcgctt                                                                                             <  1:246992/150‑3 (MQ=255)
                                                          gcgcACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGc                                                                                            <  1:16948/150‑1 (MQ=255)
                                                          gcgcACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGc                                                                                            <  2:341756/150‑1 (MQ=255)
                                                          gcgcACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGc                                                                                            <  2:81707/150‑1 (MQ=255)
                                                              aCCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGc                                                                                        >  2:1084470/1‑150 (MQ=255)
                                                               ccAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCg                                                                                       >  2:1317121/1‑150 (MQ=255)
                                                                  aaTTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGTATAAAACGCGGCAAGCGTTg                                                                                    <  1:656316/150‑1 (MQ=255)
                                                                    tttCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCa                                                                                  >  2:327065/1‑150 (MQ=255)
                                                                       ctcCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATcc                                                                               >  2:1291393/1‑150 (MQ=255)
                                                                        tcCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCg                                                                              >  1:412971/1‑150 (MQ=255)
                                                                          cTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGc                                                                            >  1:935063/1‑150 (MQ=255)
                                                                           tGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCa                                                                           >  2:1493480/1‑150 (MQ=255)
                                                                           tGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCa                                                                           >  1:1492470/1‑150 (MQ=255)
                                                                               gTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTg                                                                       >  2:1180437/1‑150 (MQ=255)
                                                                                tGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGc                                                                      <  1:1152785/150‑1 (MQ=255)
                                                                                   cgcgCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCAcc                                                                   >  2:481624/1‑150 (MQ=255)
                                                                                       cATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCgcgc                                                               >  1:228742/1‑150 (MQ=255)
                                                                                        aTTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGcc                                                              >  1:125531/1‑150 (MQ=255)
                                                                                         tttGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCa                                                             <  2:1537815/150‑1 (MQ=255)
                                                                                            gTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTg                                                          >  1:341763/1‑150 (MQ=255)
                                                                                             tataTGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTgg                                                         >  1:1310856/1‑150 (MQ=255)
                                                                                             tataTGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTgg                                                         >  1:1163452/1‑150 (MQ=255)
                                                                                                atGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCgg                                                      >  2:311098/1‑150 (MQ=255)
                                                                                                atGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCgg                                                      >  1:217978/1‑150 (MQ=255)
                                                                                                atGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCgg                                                      >  2:512453/1‑150 (MQ=255)
                                                                                                 tgttgtATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGa                                                     >  1:825797/1‑150 (MQ=255)
                                                                                                  gttgtATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGAt                                                    >  2:936005/1‑150 (MQ=255)
                                                                                                  gttgtATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGAt                                                    <  1:931503/150‑1 (MQ=255)
                                                                                                   ttgtATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATg                                                   >  1:88862/1‑150 (MQ=255)
                                                                                                       ataACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATgcggc                                               >  1:246483/1‑150 (MQ=255)
                                                                                                        taACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATgcggcg                                              >  1:1057452/1‑150 (MQ=255)
                                                                                                         aaCGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATACGCCCTACATGTGCAATCCTCATGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAAATGCACCGCGCCACTGGCGGATGCGGCGt                                             >  2:665789/1‑150 (MQ=255)
                                                                                                           cGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGa                                           <  2:672574/150‑1 (MQ=255)
                                                                                                             aTTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAAc                                         >  1:1403947/1‑150 (MQ=255)
                                                                                                             aTTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAAc                                         >  2:380585/1‑150 (MQ=255)
                                                                                                                   tCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTa                                   >  1:149157/1‑150 (MQ=255)
                                                                                                                       tAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCt                                     >  2:146953/1‑144 (MQ=255)
                                                                                                                       tAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCt                                     <  1:146953/144‑1 (MQ=255)
                                                                                                                       tAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCg                               <  2:1514910/150‑1 (MQ=255)
                                                                                                                         aGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATAAGGCGTTCACGCCGCATCCGcc                             <  1:1428664/150‑1 (MQ=18)
                                                                                                                          gCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGccc                            >  1:1191488/1‑150 (MQ=255)
                                                                                                                           cAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCt                           >  1:576574/1‑150 (MQ=255)
                                                                                                                           cAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCt                           >  1:703427/1‑150 (MQ=255)
                                                                                                                            aataatTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTa                          >  2:1263492/1‑150 (MQ=255)
                                                                                                                               aatTTGCAACTGCTGGCGGATGTGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACAt                       >  1:219163/1‑150 (MQ=255)
                                                                                                                                  ttGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACAtgtg                    <  2:1351060/150‑1 (MQ=255)
                                                                                                                                    gCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACAtgtgtg                  <  1:62447/150‑1 (MQ=255)
                                                                                                                                    gCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACAtgtgtg                  <  2:125531/150‑1 (MQ=255)
                                                                                                                                     cAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATgtgtgt                 >  1:931357/1‑150 (MQ=255)
                                                                                                                                      aaCTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGtt                >  1:1282967/1‑150 (MQ=255)
                                                                                                                                        ctgctgGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTcc              <  1:154721/150‑1 (MQ=255)
                                                                                                                                          gctgGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCg            <  1:15620/150‑1 (MQ=255)
                                                                                                                                           ctgGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGt           <  2:217978/150‑1 (MQ=255)
                                                                                                                                           agggCGGATAAGGCGTTTACGCCGCATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATgc                                                  >  1:1530345/3‑111 (MQ=17)
                                                                                                                                           agggCGGATAAGGCGTTTACGCCGCATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAAc                                         >  2:476736/3‑120 (MQ=17)
                                                                                                                                           agggCGGATAAGGCGTTTACGCCGCATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAAc                                         <  1:476736/118‑1 (MQ=17)
                                                                                                                                           agggCGGATAAGGCGTTTACGCCGCATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGt           >  2:56740/3‑150 (MQ=17)
                                                                                                                                              gCGGATGCGGCGTAAACGCCGTATCCGCCCTACATGGGCAAGCCTCCTGGCCGGATAAAACGCGGCAAGCGGTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAgg        <  2:1254445/150‑1 (MQ=255)
                                                                                                                                                ggATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTc      <  1:986807/150‑1 (MQ=255)
                                                                                                                                                  aTGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCgg    <  1:93347/150‑1 (MQ=255)
                                                                                                                                                    gcggcgTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAGAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGAt  <  2:1026994/150‑1 (MQ=255)
                                                                                                                                                    gcggcgTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAACACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGAt  >  1:523872/1‑150 (MQ=255)
                                                                                                                                                    gcggcgTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGAt  <  2:95950/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATATGTTGTATAACGATTTTTTCCGTAAGCAATAATTTGCAACTGCTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGAT  >  NC_000913/3231111‑3231408

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