Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,504,875 T→G 6.3% intergenic (+92/+60) frlR → / ← yhfS putative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,504,8750TG 6.3% 404.0 / 13.1 122intergenic (+92/+60)frlR/yhfSputative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base T (66/47);  new base G (4/4);  total (70/52)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.99e-01

Only 100 of 133 total aligned reads displayed.

GAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTC  >  NC_000913/3504727‑3505021
                                                                                                                                                    |                                                                                                                                                  
gAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCctct                                                                                                                                                    >  1:1467024/1‑149 (MQ=255)
 aaCTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCtctc                                                                                                                                                   <  2:334932/149‑1 (MQ=255)
 aaCTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCtctc                                                                                                                                                   <  4:834370/149‑1 (MQ=255)
  aCTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAAGAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTcc                                                                                                                                                  <  3:96689/149‑1 (MQ=255)
    tcttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTcccc                                                                                                                                                >  2:993139/1‑149 (MQ=255)
    tcttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTcccc                                                                                                                                                >  3:333238/1‑149 (MQ=255)
    tcttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCGcccc                                                                                                                                                >  2:695139/1‑149 (MQ=255)
       tctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCaaa                                                                                                                                             <  3:148630/149‑1 (MQ=255)
         tGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAgg                                                                                                                                           <  4:908391/149‑1 (MQ=255)
          gCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAggg                                                                                                                                          <  3:381833/149‑1 (MQ=255)
          gCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCGcccctttggg                                                                                                                                          <  1:427797/149‑7 (MQ=255)
           ccGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAgggg                                                                                                                                         <  1:1055604/149‑1 (MQ=255)
               gcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAg                                                                                                                                     >  1:700507/1‑149 (MQ=255)
                ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAgg                                                                                                                                    <  2:700507/149‑1 (MQ=255)
                ccAATCGCATCACCTTAACTAATGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGAGAgg                                                                                                                                    <  1:1466503/149‑1 (MQ=38)
                 cAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAggg                                                                                                                                   <  2:621209/149‑1 (MQ=255)
                 cAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAggg                                                                                                                                   <  2:1216983/149‑1 (MQ=255)
                  aaTCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAgggg                                                                                                                                  >  4:582460/1‑149 (MQ=255)
                             cTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTa                                                                                                                       <  1:594389/149‑1 (MQ=255)
                              ttAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTAc                                                                                                                      <  4:760885/149‑1 (MQ=255)
                                      tGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTa                                                                                                              <  3:865628/149‑1 (MQ=255)
                                       gATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAg                                                                                                             >  1:1026738/1‑149 (MQ=255)
                                        ataataAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAgg                                                                                                                                    >  4:462691/1‑125 (MQ=255)
                                        ataataAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAgg                                                                                                                                    <  3:462691/125‑1 (MQ=255)
                                        ataataAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGc                                                                                                            <  2:1029302/149‑1 (MQ=255)
                                        ataataAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGc                                                                                                            >  2:1457657/1‑149 (MQ=255)
                                         taataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCt                                                                                                           <  2:453279/149‑1 (MQ=255)
                                          aataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCtt                                                                                                          >  2:182608/1‑149 (MQ=255)
                                             aaaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGAc                                                                                                       <  3:1296432/149‑1 (MQ=255)
                                              aaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACg                                                                                                      >  2:592962/1‑149 (MQ=255)
                                               aaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTa                                                                                                                       <  1:935669/131‑1 (MQ=255)
                                               aaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACgg                                                                                                     >  2:1230840/1‑149 (MQ=255)
                                                 gACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTg                                                                                                   <  4:333238/149‑1 (MQ=255)
                                                 gACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTg                                                                                                   >  3:1434413/1‑149 (MQ=255)
                                                  aCATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCGCCCTATCCCTCCCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGc                                                                                                  <  4:1213420/149‑1 (MQ=255)
                                                    aTTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGAGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTa                                                                                                >  1:1308235/1‑149 (MQ=255)
                                                       aaCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATcc                                                                                             >  4:1144109/1‑149 (MQ=255)
                                                           gTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCa                                                                                         <  3:912240/149‑1 (MQ=255)
                                                              ggCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTa                                                                                      >  4:1414716/1‑149 (MQ=255)
                                                               gCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTaa                                                                                     >  2:1093174/1‑149 (MQ=255)
                                                                ccGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTAACCAGCATTaaa                                                                                    <  1:1099373/149‑1 (MQ=255)
                                                                   gATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAg                                                                                                                                            <  1:706198/90‑1 (MQ=255)
                                                                   gATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGGCCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATAc                                                                                 >  3:1277024/1‑149 (MQ=255)
                                                                     tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                                               <  4:150455/149‑1 (MQ=255)
                                                                      aaGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                              <  4:1421417/149‑1 (MQ=255)
                                                                        gATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCa                                                                            >  1:604664/1‑149 (MQ=255)
                                                                        gATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCa                                                                            >  1:843331/1‑149 (MQ=255)
                                                                         aTGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCaa                                                                           <  2:38289/149‑1 (MQ=255)
                                                                             gcCAGCATCGCCGGATGCGATGCTGGCGCATTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATAct                                                                       >  1:272759/1‑149 (MQ=38)
                                                                             gcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAgg                                                                                                                                    <  3:1398195/88‑1 (MQ=255)
                                                                             gcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATAct                                                                       >  2:976073/1‑149 (MQ=255)
                                                                              ccAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATActc                                                                      >  2:914459/1‑149 (MQ=255)
                                                                               cAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATActct                                                                     >  2:765991/1‑149 (MQ=255)
                                                                                aGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACtctc                                                                    >  3:624187/1‑149 (MQ=255)
                                                                                 gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                              <  1:175483/138‑1 (MQ=255)
                                                                                 gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCa                                                                   <  1:993139/149‑1 (MQ=255)
                                                                                  catcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                  <  2:758492/149‑1 (MQ=255)
                                                                                  catcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                  >  2:146320/1‑149 (MQ=255)
                                                                                      gcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCaa                                                              >  3:755227/1‑149 (MQ=255)
                                                                                        atcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCaaaa                                                            >  1:1033287/1‑149 (MQ=255)
                                                                                        atcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCaaaa                                                            >  3:719897/1‑149 (MQ=255)
                                                                                            ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCg                                                        >  3:641108/1‑149 (MQ=255)
                                                                                                  gCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCAc                                                  >  1:80887/1‑149 (MQ=255)
                                                                                                   cTGGCGCGTTTAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCGCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTTACGCAAAACCCGCAGCAcc                                                 <  2:582062/149‑1 (MQ=255)
                                                                                                      gcgcGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGtc                                              >  1:786700/1‑149 (MQ=255)
                                                                                                         cgTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTctct                                           <  2:843331/149‑1 (MQ=255)
                                                                                                          gTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCtt                                          <  4:719897/149‑1 (MQ=255)
                                                                                                                ttttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTCCGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGc                                    >  3:885163/1‑149 (MQ=255)
                                                                                                                ttttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGc                                    >  4:1217273/1‑149 (MQ=255)
                                                                                                                ttttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTATCACGCAAAATCCGCAGCACCGTCTCTTCACCGc                                    >  2:500595/1‑149 (MQ=255)
                                                                                                                      aTCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGAt                              >  1:908310/1‑149 (MQ=255)
                                                                                                                         ccGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                                               <  4:583742/97‑1 (MQ=255)
                                                                                                                         ccGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCg                           >  4:1375291/1‑149 (MQ=255)
                                                                                                                         ccGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCgg                          >  2:360681/1‑149 (MQ=255)
                                                                                                                         ccGTACGTTCCCCTCACCCGAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTGAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTAGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCg                           >  4:919865/1‑149 (MQ=255)
                                                                                                                          cGTACGTTCCCCTCACCCTGACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCAATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCgg                          <  2:700757/149‑1 (MQ=255)
                                                                                                                          cGTACGTTCCCCTCACCCTAACCCTCTCCCCATAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCgg                          <  4:1106995/149‑1 (MQ=255)
                                                                                                                          cGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCgg                          >  2:1119151/1‑149 (MQ=255)
                                                                                                                           gTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCggg                         >  2:1067608/1‑149 (MQ=255)
                                                                                                                               gTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTaa                     >  3:1163110/1‑149 (MQ=255)
                                                                                                                               gTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTaa                     >  2:263079/1‑149 (MQ=255)
                                                                                                                               gTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTaa                     >  2:1136779/1‑149 (MQ=255)
                                                                                                                                    ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATa                                                                         <  2:1355894/92‑1 (MQ=255)
                                                                                                                                    ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATa                                                                         >  1:1355894/1‑92 (MQ=255)
                                                                                                                                    ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGa                >  2:839939/1‑149 (MQ=255)
                                                                                                                                     cTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa               >  4:190084/1‑149 (MQ=255)
                                                                                                                                       cACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATc             <  2:272759/149‑1 (MQ=255)
                                                                                                                                        aCCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGCTTAATGCGAATCg            >  2:28674/1‑149 (MQ=255)
                                                                                                                                          ccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGca          <  1:765991/149‑1 (MQ=255)
                                                                                                                                          ccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGca          >  1:713063/1‑149 (MQ=255)
                                                                                                                                          ccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGca          >  3:626789/1‑149 (MQ=255)
                                                                                                                                          ccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGCATTAACCCGAATCGca          <  4:755227/149‑1 (MQ=255)
                                                                                                                                           cTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcac         >  1:580678/1‑149 (MQ=255)
                                                                                                                                             tccccTCGCCCCTTTGGGGAGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACaa       >  2:76663/3‑149 (MQ=255)
                                                                                                                                               cccTCGCCCCTTTGGGGAGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATg     >  2:1452465/1‑149 (MQ=255)
                                                                                                                                               cccTCGCCCCTTTGGGGAGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATg     >  1:199688/1‑149 (MQ=255)
                                                                                                                                                 ctctCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGtt   >  4:21650/1‑149 (MQ=255)
                                                                                                                                                 ctctCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGtt   <  3:464303/149‑1 (MQ=255)
                                                                                                                                                 ctctCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAACACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGtt   >  1:121450/1‑149 (MQ=255)
                                                                                                                                                  tctcCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc  >  3:785710/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTC  >  NC_000913/3504727‑3505021

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: