Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 9,907 A→C 8.8% intergenic (+14/+21) mog → / ← satP molybdochelatase incorporating molybdenum into molybdopterin/succinate‑acetate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009139,9070AC 8.8% 545.1 / 39.9 171intergenic (+14/+21)mog/satPmolybdochelatase incorporating molybdenum into molybdopterin/succinate‑acetate transporter
Reads supporting (aligned to +/- strand):  ref base A (73/83);  new base C (6/9);  total (79/92)
Fisher's exact test for biased strand distribution p-value = 7.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.11e-01

Only 100 of 172 total aligned reads displayed.

AACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATT  >  NC_000913/9759‑10051
                                                                                                                                                     |                                                                                                                                                
 aCGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAg                                                                                                                                                 >  4:630272/1‑149 (MQ=255)
   gTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCg                                                                                                                                               >  3:565400/1‑149 (MQ=255)
   gTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTCGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTCGCGAATAAAAACAT‑CCCCCCGAGCg                                                                                                                                               >  3:674665/1‑149 (MQ=255)
    ttGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCgg                                                                                                                                              >  1:314390/1‑149 (MQ=255)
     tgtgGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCggg                                                                                                                                             <  1:644699/149‑1 (MQ=255)
       tgGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCggggg                                                                                                                                           <  3:588842/149‑1 (MQ=255)
       tgGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCggggg                                                                                                                                           <  3:392889/149‑1 (MQ=255)
            cacGGTATTTTTGTCAGCGTACCGTACTGCATTCAGTTGCTGCAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGCTCGGGGGGAttt                                                                                                                                      <  4:983693/149‑3 (MQ=255)
            cacGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCGCCCGAGCGGGGGGAtct                                                                                                                                      >  2:1244108/1‑149 (MQ=255)
             acGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGCTCGGGGGGAtttt                                                                                                                                     <  4:37410/149‑4 (MQ=255)
             acGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGCTCGGGGGGAtttt                                                                                                                                     <  4:1147871/149‑4 (MQ=255)
              cGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCa                                                                                                                                    >  1:890130/1‑149 (MQ=255)
                 tATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCaaaa                                                                                                                                 <  1:309472/149‑1 (MQ=255)
                 tATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCaaaa                                                                                                                                 <  2:706856/149‑1 (MQ=255)
                    ttttGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACaa                                                                                                                              <  4:10507/149‑1 (MQ=255)
                    ttttGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACaa                                                                                                                              <  2:596504/149‑1 (MQ=255)
                        gCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAg                                                                                                                          >  1:1360793/1‑149 (MQ=255)
                         ccAGCGTACCGTACTGCATTCAGTTGCTGGAGGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGt                                                                                                                         <  3:603722/149‑1 (MQ=255)
                         ccAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGt                                                                                                                         <  2:127639/149‑1 (MQ=255)
                           aGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTgg                                                                                                                       <  3:448435/149‑1 (MQ=255)
                           aGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTgg                                                                                                                       >  1:372436/1‑149 (MQ=255)
                           aGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTgg                                                                                                                       <  2:1118289/149‑1 (MQ=255)
                           aGCGGACCGTAGTGTATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGATGAGTGCAAGACGCCACGTTAGCGAATAAAAAAAT‑CCCCACGAGCGGGGGGAGCTCAAAACAATTAGTgg                                                                                                                       <  4:902373/149‑1 (MQ=255)
                                 ccgtacTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATGTCAAAACAATTAGTGGGATTCa                                                                                                                 <  2:747185/149‑1 (MQ=255)
                                 ccgtacTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCa                                                                                                                 <  1:650769/149‑1 (MQ=255)
                                 ccgtacTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCa                                                                                                                 <  2:162619/149‑1 (MQ=255)
                                    tacTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCa                                                                                                              >  2:855428/1‑149 (MQ=255)
                                      cTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAAt                                                                                                            >  4:288947/1‑149 (MQ=255)
                                       tGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATc                                                                                                           <  2:807472/149‑1 (MQ=255)
                                       tGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATc                                                                                                           <  3:1037445/149‑1 (MQ=255)
                                          aTTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGc                                                                                                        >  1:1276244/1‑149 (MQ=255)
                                          aTTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGc                                                                                                        <  2:138279/149‑1 (MQ=255)
                                          aTTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGc                                                                                                        <  3:48473/149‑1 (MQ=255)
                                           ttCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCa                                                                                                       <  3:139966/149‑1 (MQ=255)
                                              aGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGaa                                                                                                    <  3:155815/149‑1 (MQ=255)
                                               gTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAAc                                                                                                   >  4:642027/1‑149 (MQ=255)
                                                ttGCTGGAAGGGCCATACGTTGAAAGGGCACCGGAAGTGGTTGAAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACg                                                                                                  >  4:1041474/1‑149 (MQ=255)
                                                ttGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACg                                                                                                  >  3:396890/1‑149 (MQ=255)
                                                  gCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGt                                                                                                >  1:792119/1‑149 (MQ=255)
                                                  gCTGGAAGGGACATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCCATCGGCAGAACGGt                                                                                                <  4:702726/149‑1 (MQ=255)
                                                   cTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTg                                                                                               <  3:257509/149‑1 (MQ=255)
                                                     ggAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCg                                                                                             >  4:1388569/1‑149 (MQ=255)
                                                         gggCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCa                                                                                         <  1:1226696/149‑1 (MQ=255)
                                                          ggCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCaa                                                                                        >  1:554849/1‑149 (MQ=255)
                                                            ccATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAAc                                                                                      >  1:18221/1‑149 (MQ=255)
                                                            ccATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAAc                                                                                      >  2:271762/1‑149 (MQ=255)
                                                            ccATACGGTGTAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAAc                                                                                      >  2:1150641/1‑149 (MQ=255)
                                                              aTACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGCTCgggggg                                                                                                                                          >  4:1402266/1‑95 (MQ=255)
                                                              aTACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTg                                                                                    <  3:399428/149‑1 (MQ=255)
                                                               tACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGc                                                                                   >  4:1302745/1‑149 (MQ=255)
                                                                 cGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTc                                                                                 <  4:1092698/149‑1 (MQ=255)
                                                                    tGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGtt                                                                              <  3:182969/149‑1 (MQ=255)
                                                                      aaaCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCa                                                                            <  2:1360793/149‑1 (MQ=255)
                                                                        aCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGt                                                                          <  3:928022/149‑1 (MQ=255)
                                                                        aCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGt                                                                          <  2:1276244/149‑1 (MQ=255)
                                                                          ggCACCGGAAGTGGGTGCAGCATACAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CGCCCCGAGCGGGGGGCTCTCAAAACAATTAGTGGGATTCGCCAATCGGCAGAACGGTGCGACCAAACTACTCGTTCAGTAc                                                                        >  2:908425/1‑149 (MQ=255)
                                                                           gCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCCGCTCGGGGGGAttt                                                                                                                                      <  4:500108/87‑3 (MQ=255)
                                                                           gCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACt                                                                       <  3:1302745/149‑1 (MQ=255)
                                                                             aCCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGAt                                                                                                                                        <  3:1018825/82‑1 (MQ=255)
                                                                             aCCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTc                                                                     >  3:1375541/1‑149 (MQ=255)
                                                                              ccGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACaa                                                                                                                              >  3:694597/1‑91 (MQ=255)
                                                                              ccGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACaa                                                                                                                              <  4:694597/91‑1 (MQ=255)
                                                                                 gAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAccc                                                                 >  1:470769/1‑149 (MQ=255)
                                                                                        ttGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCTAGCGGCAGAACGGTTCGACCAAACTGCTCGTTCAGTACTTCATCCATCGCCa                                                          >  4:445987/1‑149 (MQ=255)
                                                                                        ttGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCa                                                          <  2:1298258/149‑1 (MQ=255)
                                                                                        ttGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCa                                                          <  1:842933/149‑1 (MQ=255)
                                                                                           cagcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagat                                                       <  4:396890/149‑1 (MQ=255)
                                                                                            agcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCaa                                                                                                             <  3:1136992/94‑1 (MQ=255)
                                                                                                ttCAGGCCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGAtt                                                  >  4:1179919/1‑149 (MQ=255)
                                                                                                ttCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGAtt                                                  <  3:1388569/149‑1 (MQ=255)
                                                                                                    gACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTgcgc                                              >  4:78379/1‑149 (MQ=255)
                                                                                                     aCCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCt                                             <  1:1063019/149‑1 (MQ=255)
                                                                                                     aCCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCt                                             >  2:890471/1‑149 (MQ=255)
                                                                                                           gagTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTACTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCAcc                                       >  4:839768/1‑149 (MQ=255)
                                                                                                             gTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGc                                     >  4:1268744/1‑149 (MQ=255)
                                                                                                                  aGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTc                                                                             >  1:272714/1‑104 (MQ=255)
                                                                                                                  aGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTc                                                                             <  2:272714/104‑1 (MQ=255)
                                                                                                                  aGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCGGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATc                                >  1:1116180/1‑149 (MQ=255)
                                                                                                                     cgcgACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGc                             <  4:1037892/149‑1 (MQ=255)
                                                                                                                     cgcgACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGc                             <  3:1179919/149‑1 (MQ=255)
                                                                                                                     cgcgACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGc                             <  2:500616/149‑1 (MQ=255)
                                                                                                                     cgcgACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGc                             <  2:1102641/149‑1 (MQ=255)
                                                                                                                             tAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCa                     <  3:78379/149‑1 (MQ=255)
                                                                                                                                  aaTAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCgg                >  2:780525/1‑149 (MQ=255)
                                                                                                                                      aaaaaaT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCaaa            >  4:571494/1‑149 (MQ=255)
                                                                                                                                       aaaaaT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAg           >  4:253410/1‑149 (MQ=255)
                                                                                                                                       aaaaaT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAg           >  1:887223/1‑149 (MQ=255)
                                                                                                                                        aaaaT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAga          >  2:1280594/1‑148 (MQ=255)
                                                                                                                                        aaaaT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGt          <  1:1475140/149‑1 (MQ=255)
                                                                                                                                        aaaaG‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGt          <  4:167652/149‑1 (MQ=255)
                                                                                                                                          aaT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTgg        >  1:1306461/1‑149 (MQ=255)
                                                                                                                                           aT‑CCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGa                                  <  1:1276801/122‑1 (MQ=255)
                                                                                                                                           aT‑CCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGga       >  1:711600/1‑149 (MQ=255)
                                                                                                                                           aT‑CCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGga       >  1:988751/1‑149 (MQ=255)
                                                                                                                                           aT‑CCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGga       >  2:442227/1‑149 (MQ=255)
                                                                                                                                            t‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGgat      >  3:246869/1‑149 (MQ=255)
                                                                                                                                               ccccnGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGgatga    >  1:520246/1‑149 (MQ=255)
                                                                                                                                                ccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGgatgat   >  3:457634/1‑149 (MQ=255)
                                                                                                                                                ccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGgatgat   >  3:766164/1‑149 (MQ=255)
                                                                                                                                                 cccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGAtt  >  2:584440/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
AACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAT‑CCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATT  >  NC_000913/9759‑10051

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