Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 876,569 A→T 5.6% 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.6% 581.4 / 18.6 160N412I (AAT→ATT) ‡yliFputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (76/75);  new base T (4/5);  total (80/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01

Only 100 of 178 total aligned reads displayed.

GGTAAAGTCTACCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGATTTTGCCCGGCCTGAATTAATTA  >  NC_000913/876421‑876712
                                                                                                                                                    |                                                                                                                                               
 gTAAAGTCTACCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAaat                                                                                                                                                >  1:1133540/1‑149 (MQ=255)
       tCTACCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGt                                                                                                                                          >  1:919659/1‑149 (MQ=255)
         tACCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTcc                                                                                                                                        <  2:1133540/149‑1 (MQ=255)
          aCCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAACAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCt                                                                                                                                       >  4:1552432/1‑149 (MQ=255)
                                  cGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAAt                                                                                                                   >  3:1854257/1‑145 (MQ=255)
                                  cGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTg                                                                                                               >  1:57102/1‑149 (MQ=255)
                                   ggTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATg                                                                                                                  <  2:244161/145‑1 (MQ=255)
                                   ggTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTgg                                                                                                              <  3:1552432/149‑1 (MQ=255)
                                   ggTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTgg                                                                                                              >  2:140142/1‑149 (MQ=255)
                                    gTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTAcc                                                                                                                                    <  2:576040/126‑1 (MQ=255)
                                     tGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCtt                                                                                                                         >  3:528886/1‑136 (MQ=255)
                                       aTTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGa                                                                                                          <  1:140142/149‑1 (MQ=255)
                                        tttATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCATTTGAATGtgc                                                                                                                                        <  4:1713771/118‑3 (MQ=255)
                                        tttATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCATTTGAATGtgc                                                                                                                                        >  3:1713771/1‑116 (MQ=255)
                                        tttATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGAt                                                                                                         >  3:29889/1‑149 (MQ=255)
                                         ttATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGaa                                                                                                                      >  1:1045277/1‑135 (MQ=255)
                                               gggTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTg                                                                                                               <  3:486778/136‑1 (MQ=255)
                                                 gTGGAACGCTGGAAGAGTTGCTATCGATCCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTcc                                                                                                                                        >  2:183/1‑109 (MQ=255)
                                                     aaCGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAaatgaa                                                                                                                                             <  1:839264/100‑1 (MQ=255)
                                                     aaCGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAaatgaa                                                                                                                                             >  2:839264/1‑100 (MQ=255)
                                                     aaCGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGtata                                                                                            >  2:1564147/1‑149 (MQ=255)
                                                     aaCGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAACGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGAt                                                                                                         <  1:585589/136‑1 (MQ=255)
                                                     aaCGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAACGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGAt                                                                                                         >  2:585589/1‑136 (MQ=255)
                                                        gCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAaga                                                                                         <  4:1270428/149‑1 (MQ=255)
                                                        gCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAaga                                                                                         >  4:289371/1‑149 (MQ=255)
                                                            gAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATaa                                                                                           <  3:812746/143‑1 (MQ=255)
                                                            gAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTaa                                                                                     <  1:1564147/149‑1 (MQ=255)
                                                            gAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTaa                                                                                     >  3:913449/1‑149 (MQ=255)
                                                            gAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTaa                                                                                     <  1:577108/149‑1 (MQ=255)
                                                              agagTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGaa                                                                                  >  2:921906/1‑149 (MQ=255)
                                                                  ttGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAaatgaat                                                                                                                                            <  2:77731/88‑1 (MQ=255)
                                                                  ttGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAaatgaat                                                                                                                                            >  1:77731/1‑88 (MQ=255)
                                                                  ttGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGaaaaa                                                                               >  2:627251/1‑149 (MQ=255)
                                                                  ttGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGaaaaa                                                                               <  2:1852108/149‑1 (MQ=255)
                                                                  ttGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAAc                                                                              <  1:921906/149‑1 (MQ=255)
                                                                   tGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAAc                                                                              >  4:1845848/1‑149 (MQ=255)
                                                                     cTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACgg                                                                            >  2:871875/1‑149 (MQ=255)
                                                                       aTCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTaa                                                                                     >  3:584619/1‑138 (MQ=255)
                                                                        tCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTcc                                                                                                                                        <  2:326636/86‑1 (MQ=255)
                                                                        tCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATaa                                                                                           <  1:183524/131‑1 (MQ=255)
                                                                         cGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAgg                                                                        >  4:592490/1‑149 (MQ=255)
                                                                            ttCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTg                                                                                                                        >  4:572927/1‑98 (MQ=255)
                                                                             tCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGcac                                                                    <  4:913449/149‑1 (MQ=255)
                                                                                aaaaGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAaga                                                                                         <  4:1665978/125‑1 (MQ=255)
                                                                                aaaaGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAaga                                                                                         >  3:1665978/1‑125 (MQ=255)
                                                                                 aaaGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAaatgaa                                                                                                                                             >  3:311179/1‑72 (MQ=255)
                                                                                 aaaGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAaatgaa                                                                                                                                             <  4:311179/72‑1 (MQ=255)
                                                                                 aaaGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGAt                                                                <  3:289371/149‑1 (MQ=255)
                                                                                     aaaTCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACgg                                                                            >  1:449390/1‑133 (MQ=255)
                                                                                       aTCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTc                                                                                                                                         <  3:1173876/70‑1 (MQ=255)
                                                                                       aTCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTaa                                                                                     <  4:1479707/122‑1 (MQ=255)
                                                                                       aTCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTaa                                                                                     >  3:1479707/1‑122 (MQ=255)
                                                                                          gTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCa                                                       <  3:1845848/149‑1 (MQ=255)
                                                                                           ttCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAagaga                                                                                       >  4:1247213/1‑115 (MQ=255)
                                                                                           ttCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGt                                                                                       <  3:1247213/116‑1 (MQ=255)
                                                                                           ttCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAg                                                      >  4:1043825/1‑149 (MQ=255)
                                                                                           ttCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAg                                                      >  2:1077244/1‑149 (MQ=255)
                                                                                              aTTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAaga                                                                                         >  1:233746/1‑111 (MQ=255)
                                                                                              aTTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTaa                                                   <  2:458209/149‑1 (MQ=255)
                                                                                                ttCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATc                                                 <  1:17377/149‑1 (MQ=255)
                                                                                                    aGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCatt                                              >  1:1861032/1‑148 (MQ=255)
                                                                                                       tGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTc                                          >  3:1492090/1‑149 (MQ=255)
                                                                                                         ggAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCa                                                       <  3:204442/134‑1 (MQ=255)
                                                                                                         ggAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCa                                                       >  4:204442/1‑134 (MQ=255)
                                                                                                          gAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCa                                       <  2:1633078/149‑1 (MQ=255)
                                                                                                           aaGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAg                                      >  3:651155/1‑149 (MQ=255)
                                                                                                                ttAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCgg                                 >  1:735694/1‑149 (MQ=255)
                                                                                                                ttAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCgg                                 <  3:592490/149‑1 (MQ=255)
                                                                                                                   aTTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCtt                                                                                                                         >  2:1397342/1‑58 (MQ=255)
                                                                                                                   aTTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGc                                                            >  4:240695/1‑119 (MQ=255)
                                                                                                                       aaGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGc                                                                                                  >  4:4301/1‑77 (MQ=255)
                                                                                                                       aaGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGc                                                                                                  <  3:4301/77‑1 (MQ=255)
                                                                                                                         gTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGta                                                                                              >  3:11074/1‑79 (MQ=255)
                                                                                                                         gTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAAGCCGGACAACTGag                        >  4:517216/1‑148 (MQ=255)
                                                                                                                          tACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGGCaaaatcat                       >  3:1452705/1‑144 (MQ=255)
                                                                                                                          tACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGAtt                       <  3:73557/149‑1 (MQ=255)
                                                                                                                           aCCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAg                                                             <  4:82265/110‑1 (MQ=255)
                                                                                                                           aCCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAg                                                             >  3:82265/1‑110 (MQ=255)
                                                                                                                           aCCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCa                                       >  1:1814338/1‑132 (MQ=255)
                                                                                                                           aCCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCa                                       <  2:1814338/132‑1 (MQ=255)
                                                                                                                            ccAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGa                         >  2:193064/1‑145 (MQ=255)
                                                                                                                            ccAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGAtttt                     <  4:1524682/149‑1 (MQ=255)
                                                                                                                            ccAGTATCGGTGTAGCACATTCAAACGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGAtttt                     <  4:651155/149‑1 (MQ=255)
                                                                                                                             cAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCt                                                           >  1:262744/1‑110 (MQ=255)
                                                                                                                                tATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTaa                                                                                     >  1:890154/1‑81 (MQ=255)
                                                                                                                                  tCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGATTTTGCCCgg               <  4:43366/149‑1 (MQ=255)
                                                                                                                                        tAGGACATTCATTTGAATGTGCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGAt                        >  2:592367/1‑134 (MQ=255)
                                                                                                                                          ggacaTTCATTTGAATGTGCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATc                                                 <  4:781232/105‑1 (MQ=255)
                                                                                                                                          gCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGAtt                       >  2:647829/1‑133 (MQ=255)
                                                                                                                                           gacaTTCATTTGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGa                         >  3:1321386/2‑130 (MQ=255)
                                                                                                                                            acaTTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAaga                                                                                         <  4:1352661/65‑1 (MQ=255)
                                                                                                                                             caTTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGATTTTGCCCGGCCTGaattaat    <  1:60099/149‑1 (MQ=255)
                                                                                                                                              aTTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGa                                                                                                       >  4:1111680/1‑49 (MQ=255)
                                                                                                                                              aTTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGATTTTGCCCGGCCTGaattaatt   >  1:660353/1‑149 (MQ=255)
                                                                                                                                               ttCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCa                                       >  4:683269/1‑112 (MQ=255)
                                                                                                                                               ttCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGATTTTGCCCGGCCTGAattaatta  >  2:1303660/1‑149 (MQ=255)
                                                                                                                                                 cAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGATTTTGCCCGGCCTGa          <  1:265799/139‑1 (MQ=255)
                                                                                                                                                 cAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGATTTTGCCCGGCCTGa          >  2:265799/1‑139 (MQ=255)
                                                                                                                                                   aatgaatgTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCatta                                             <  2:907250/102‑1 (MQ=255)
                                                                                                                                                    atgaatgTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGATTTTg                    <  1:642668/126‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GGTAAAGTCTACCGCTTCGGCGGTGACGAATTTGCGGTGATTTATACGGGTGGAACGCTGGAAGAGTTGCTATCGATTCTAAAAGAAATCGTTCATTTCCAGGTGGGAAGCATTAATTTAAGTACCAGTATCGGTGTAGCACATTCAAATGAATGTCCTACCGTCGAACGCTTGAAAATGCTGGCGGATGAGCGGCTGTATAAGAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGATTTTGCCCGGCCTGAATTAATTA  >  NC_000913/876421‑876712

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