Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 876,569 A→T 5.5% N412I (AAT→ATT) ‡ yliF → putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913876,5690AT 5.5% 397.8 / 10.0 109N412I (AAT→ATT) ‡yliFputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (51/52);  new base T (3/3);  total (54/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01

Only 100 of 125 total aligned reads displayed.

GGTAAAGTCTACCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGATTTTGCCCGGCCTGAATTAA  >  NC_000913/876421‑876709
                                                                                                                                                    |                                                                                                                                            
ggTAAAGTCTACCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCaaa                                                                                                                                              <  1:507842/149‑1 (MQ=255)
ggTAAAGTCTACCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCaaa                                                                                                                                              <  2:855538/149‑1 (MQ=255)
 gTAAAGTCTACCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCaaa                                                                                                                                              >  3:407259/1‑148 (MQ=255)
     aGTCTACCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAaatgaat                                                                                                                                         <  4:406479/149‑1 (MQ=255)
      gTCTACCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAaatgaatg                                                                                                                                        >  4:250228/1‑149 (MQ=255)
         tACCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTcc                                                                                                                                     >  3:1156755/1‑149 (MQ=255)
              cTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAaat                                                                                                                                             >  4:521027/1‑136 (MQ=255)
              cTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAaat                                                                                                                                             <  3:521027/136‑1 (MQ=255)
                     ggTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCaaa                                                                                                                                              >  1:565708/1‑128 (MQ=255)
                                 gCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCt                                                                                                             <  4:386687/149‑1 (MQ=255)
                                   ggTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTgg                                                                                                           <  4:1156755/149‑1 (MQ=255)
                                                            gAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTgg                                                                                                           <  3:415413/124‑1 (MQ=255)
                                                            gAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTgg                                                                                                           >  4:415413/1‑124 (MQ=255)
                                                            gAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTaa                                                                                  >  1:894578/1‑149 (MQ=255)
                                                            gAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTaa                                                                                  <  2:894578/149‑1 (MQ=255)
                                                                 gTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAaga                                                                                      >  3:231139/1‑140 (MQ=255)
                                                                      tATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGa                                                                                                                    >  3:90350/1‑105 (MQ=255)
                                                                      tATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGa                                                                                                                    <  4:90350/105‑1 (MQ=255)
                                                                       aTCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGaaaaa                                                                            <  2:1451021/144‑1 (MQ=255)
                                                                       aTCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGaaaaa                                                                            >  1:1451021/1‑144 (MQ=255)
                                                                        tCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCa                                                                       <  2:493023/148‑1 (MQ=255)
                                                                        tCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCa                                                                       >  1:493023/1‑148 (MQ=255)
                                                                           aTTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGc                                                                   >  3:1278219/1‑149 (MQ=255)
                                                                               tAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAaga                                                                                      >  4:309469/1‑126 (MQ=255)
                                                                               tAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAaga                                                                                      <  3:309469/126‑1 (MQ=255)
                                                                               tAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAg                                                               >  2:488581/1‑149 (MQ=255)
                                                                                 aaaGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGc                                                                                                          >  2:730294/1‑104 (MQ=255)
                                                                                 aaaGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGc                                                                                                          <  1:730294/104‑1 (MQ=255)
                                                                                 aaaGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAggg                                                                    <  1:909589/142‑1 (MQ=255)
                                                                                 aaaGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAggg                                                                    <  3:1057653/142‑1 (MQ=255)
                                                                                 aaaGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAggg                                                                    >  2:909589/1‑142 (MQ=255)
                                                                                 aaaGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAggg                                                                    >  4:1057653/1‑142 (MQ=255)
                                                                                 aaaGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAgg                                                                     <  2:270252/141‑1 (MQ=255)
                                                                                 aaaGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAgg                                                                     >  1:270252/1‑141 (MQ=255)
                                                                                  aaGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGAt                                                                                                      >  4:301835/1‑107 (MQ=255)
                                                                                  aaGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGAt                                                                                                      <  3:301835/107‑1 (MQ=255)
                                                                                  aaGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGAtt                                                            >  4:1072855/1‑149 (MQ=255)
                                                                                   aGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAaatgaat                                                                                                                                         <  3:283303/71‑1 (MQ=255)
                                                                                   aGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAaatgaat                                                                                                                                         >  4:283303/1‑71 (MQ=255)
                                                                                   aGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAaga                                                                                      <  3:482247/122‑1 (MQ=255)
                                                                                   aGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAaga                                                                                      >  4:482247/1‑122 (MQ=255)
                                                                                    gAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAg                                                          >  1:1169434/1‑149 (MQ=255)
                                                                                        tCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTgg                                                      <  1:488581/149‑1 (MQ=255)
                                                                                           ttCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAaatgaat                                                                                                                                         >  2:367750/1‑63 (MQ=255)
                                                                                           ttCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAaatgaat                                                                                                                                         <  1:367750/63‑1 (MQ=255)
                                                                                           ttCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCtt                                                                                                                      >  3:1221847/1‑82 (MQ=255)
                                                                                           ttCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCtt                                                                                                                      <  4:1221847/82‑1 (MQ=255)
                                                                                           ttCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAg                                                                      >  4:614058/1‑130 (MQ=255)
                                                                                           ttCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAg                                                                      <  3:614058/130‑1 (MQ=255)
                                                                                           ttCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGc                                                                   <  4:527923/133‑1 (MQ=255)
                                                                                           ttCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTgg                                                      >  3:531658/1‑146 (MQ=255)
                                                                                              aTTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGaaaaa                                                                            >  1:203734/1‑121 (MQ=255)
                                                                                              aTTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGaaaaa                                                                            <  2:203734/121‑1 (MQ=255)
                                                                                               tttCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTg                                                                                                            <  2:576054/88‑1 (MQ=255)
                                                                                               tttCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTg                                                                                                            >  1:576054/1‑88 (MQ=255)
                                                                                               tttCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGaaaaa                                                                            >  4:395926/1‑120 (MQ=255)
                                                                                                 tCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCa                                             <  3:1072855/149‑1 (MQ=255)
                                                                                                    aGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTg                                                                                                            >  2:1332838/1‑83 (MQ=255)
                                                                                                    aGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGAACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTg                                                                                                            <  1:1332838/83‑1 (MQ=255)
                                                                                                     ggTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCa                                                                       <  3:6055/119‑1 (MQ=255)
                                                                                                     ggTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCa                                                                       >  4:6055/1‑119 (MQ=255)
                                                                                                     ggTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCattat                                         <  2:1169434/149‑1 (MQ=255)
                                                                                                      gTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCattatt                                        >  1:979449/1‑149 (MQ=255)
                                                                                                          gAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCa                                    <  2:979449/149‑1 (MQ=255)
                                                                                                               aTTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAg                                                                      <  2:494294/110‑1 (MQ=255)
                                                                                                               aTTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAg                                                                      >  1:494294/1‑110 (MQ=255)
                                                                                                               aTTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTaa                                                <  2:916691/132‑1 (MQ=255)
                                                                                                               aTTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTaa                                                >  1:916691/1‑132 (MQ=255)
                                                                                                               aTTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCatt                                           >  4:1013977/1‑137 (MQ=255)
                                                                                                                ttAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCgg                              >  2:522828/1‑149 (MQ=255)
                                                                                                                   aTTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGAt                                                                                                      >  4:1139732/1‑74 (MQ=255)
                                                                                                                   aTTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATaa                                                                                        >  2:306333/1‑88 (MQ=255)
                                                                                                                   aTTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATaa                                                                                        <  1:306333/88‑1 (MQ=255)
                                                                                                                               gTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGt                                                                                    <  3:382720/80‑1 (MQ=255)
                                                                                                                               gTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGt                                                                                    >  4:382720/1‑80 (MQ=255)
                                                                                                                               gTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGAtt                    <  1:1165653/144‑1 (MQ=255)
                                                                                                                               gTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGAtt                    >  2:1165653/1‑144 (MQ=255)
                                                                                                                                    ggTGTAGCACATACAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCa                                             >  4:984249/1‑114 (MQ=255)
                                                                                                                                      tgtAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGAc                            >  3:1464442/1‑129 (MQ=255)
                                                                                                                                      tgtAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGAc                            <  4:1464442/129‑1 (MQ=255)
                                                                                                                                        tAGGACATTCATTTGAATGTGCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGc                                                         >  2:1286428/1‑98 (MQ=38)
                                                                                                                                        tAGGACATTCATTTGAATGTGCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGc                                                         <  1:1286428/98‑1 (MQ=39)
                                                                                                                                          ggacaTTCATTTGAATGTGCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAg                                                                                 >  4:176657/3‑72 (MQ=38)
                                                                                                                                          ggacaTTCATTTGAATGTGCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAg                                                                                 <  3:176657/70‑1 (MQ=37)
                                                                                                                                          ggacaTTCATTTGAATGTGCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCatt                                           >  3:194103/3‑110 (MQ=255)
                                                                                                                                          ggacaTTCATTTGAATGTGCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCatt                                           <  4:194103/108‑1 (MQ=255)
                                                                                                                                          gCACATTCAAATGAATGTCCTACCGTCGAACGCtt                                                                                                                      >  1:584005/1‑35 (MQ=255)
                                                                                                                                          gCACATTCAAATGAATGTCCTACCGTCGAACGCtt                                                                                                                      <  2:584005/35‑1 (MQ=255)
                                                                                                                                          gCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATc                                              <  4:1461952/107‑1 (MQ=255)
                                                                                                                                          gCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCa                                    >  3:1077102/1‑117 (MQ=255)
                                                                                                                                              aTTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTa                                                           >  4:1035372/1‑90 (MQ=255)
                                                                                                                                              aTTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTa                                                           <  3:1035372/90‑1 (MQ=255)
                                                                                                                                               ttCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAaga                                                                                      <  4:1176781/62‑1 (MQ=255)
                                                                                                                                               ttCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAaga                                                                                      >  3:1176781/1‑62 (MQ=255)
                                                                                                                                               ttCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTgg                                                      <  1:199011/94‑1 (MQ=255)
                                                                                                                                               ttCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCa                                             >  2:231642/1‑103 (MQ=255)
                                                                                                                                               ttCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCa                                             <  1:231642/103‑1 (MQ=255)
                                                                                                                                               ttCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGc                                     >  2:296778/1‑111 (MQ=255)
                                                                                                                                               ttCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGc                                     <  1:296778/111‑1 (MQ=255)
                                                                                                                                               ttCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGAATAGCTgg                                                      >  2:199011/1‑94 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
GGTAAAGTCTACCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGATTTTGCCCGGCCTGAATTAA  >  NC_000913/876421‑876709

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