Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,517,438 G→A 5.5% T132M (ACG→ATG)  yfeA ← putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,517,4380GA 5.5% 371.6 / 8.3 109T132M (ACG→ATG) yfeAputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (53/50);  new base A (3/3);  total (56/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 112 total aligned reads displayed.

GTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAA  >  NC_000913/2517293‑2517586
                                                                                                                                                 |                                                                                                                                                    
gttgtAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCa                                                                                                                                                   <  3:745943/149‑1 (MQ=255)
 ttgtAAATCAGCACAGCGGTGGAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATAGCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCac                                                                                                                                                  >  2:280133/1‑149 (MQ=255)
 ttgtAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAAGCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCac                                                                                                                                                  <  3:346141/149‑1 (MQ=255)
      aaTCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCaccaaac                                                                                                                                             >  1:500167/1‑149 (MQ=255)
       aTCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCaccaaacc                                                                                                                                            >  3:474056/1‑149 (MQ=255)
       aTCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCaccaaacc                                                                                                                                            >  2:435604/1‑149 (MQ=255)
            cacaGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGACGCCGATCGGCATCaaat                                                                                                                                                >  1:955456/1‑137 (MQ=255)
            cacaGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCGCCCCACAAGATACATGCTGCATTTGACGCCGATCGGCATCaaat                                                                                                                                                <  2:955456/140‑4 (MQ=255)
                     gAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGAc                                                                                                                              <  4:474056/149‑1 (MQ=255)
                          ggCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGAc                                                                                                                              >  2:874092/1‑144 (MQ=255)
                          ggCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGAc                                                                                                                              <  1:874092/144‑1 (MQ=255)
                          ggCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGc                                                                                                                         >  1:602145/1‑149 (MQ=255)
                           gCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGtt                                                                                                                           >  3:1115698/1‑146 (MQ=255)
                           gCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGtt                                                                                                                           <  4:1115698/146‑1 (MQ=255)
                           gCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGACGCCGATCGGCATCaaat                                                                                                                                                <  1:1342492/125‑4 (MQ=255)
                           gCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGACGCCGATCGGCATCaaat                                                                                                                                                >  2:1342492/1‑122 (MQ=255)
                             ttAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCaccaaa                                                                                                                                              >  4:268758/1‑125 (MQ=255)
                             ttAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCaccaaa                                                                                                                                              <  3:268758/125‑1 (MQ=255)
                             ttAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAg                                                                                                                      >  4:802328/1‑149 (MQ=255)
                             ttAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAg                                                                                                                      >  4:84937/1‑149 (MQ=255)
                             ttAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGACGCCGATCGGCATCaaat                                                                                                                                                <  1:204326/123‑4 (MQ=255)
                             ttAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGACGCCGATCGGCATCaaat                                                                                                                                                >  2:204326/1‑120 (MQ=255)
                               aGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGAc                                                                                                                    <  3:845701/149‑1 (MQ=255)
                                gCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCa                                                                                                                                      >  2:795568/1‑130 (MQ=255)
                                 cAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGtt                                                                                                                           >  2:1266859/1‑140 (MQ=255)
                                 cAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGtt                                                                                                                           <  1:1266859/140‑1 (MQ=255)
                                  aaaTCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACat                                                                                                                  <  4:381996/148‑1 (MQ=255)
                                  aaaTCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACat                                                                                                                  >  3:381996/1‑148 (MQ=255)
                                   aaTCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACatat                                                                                                                >  4:780195/1‑149 (MQ=255)
                                    aTCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGGACAGACGTTGCCAGACATATc                                                                                                               >  4:560109/1‑149 (MQ=255)
                                       aaCGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGc                                                                                                            >  4:1285334/1‑149 (MQ=255)
                                        aCGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGGACAGACGTTGCCAGACATATCGGCt                                                                                                           <  3:560109/149‑1 (MQ=255)
                                        aCGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCt                                                                                                           <  3:84937/149‑1 (MQ=255)
                                        aCGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCt                                                                                                           >  3:1458539/1‑149 (MQ=255)
                                           aCCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCa                                                                                                                                      >  1:181959/1‑119 (MQ=255)
                                           aCCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCa                                                                                                                                      <  2:181959/119‑1 (MQ=255)
                                            ccGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTc                                                                                                       <  4:1458539/149‑1 (MQ=255)
                                              gTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCaa                                                                                                     >  4:304388/1‑149 (MQ=255)
                                                gAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGAATTTGATGCCGATCGGCGTCACCAAACCTAGCCAGAACAGACGTTGCCAGACAAATCGGCTGGTCAATc                                                                                                   <  2:500167/149‑1 (MQ=255)
                                                   aaTGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAAt                                                                                                    <  1:775035/145‑1 (MQ=255)
                                                   aaTGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAAt                                                                                                    >  2:775035/1‑145 (MQ=255)
                                                     tGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATaa                                                                                              <  3:780195/149‑1 (MQ=255)
                                                         aTCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGc                                                                                                                         <  3:488180/118‑1 (MQ=255)
                                                            cGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAAc                                                                                       <  2:602145/149‑1 (MQ=255)
                                                             gCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGc                                                                                                                                        >  2:955056/1‑99 (MQ=255)
                                                             gCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGc                                                                                                                                        <  1:955056/99‑1 (MQ=255)
                                                                tCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATcc                                                                                                  >  2:957779/1‑134 (MQ=255)
                                                                            gTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGAc                                                                                                                    <  2:397505/104‑1 (MQ=255)
                                                                            gTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGAc                                                                                                                    >  1:397505/1‑104 (MQ=255)
                                                                            gTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCa                                                                       >  4:474376/1‑149 (MQ=255)
                                                                                           gggAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGt                                                                                                        >  2:875366/1‑101 (MQ=255)
                                                                                           gggAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGt                                                                                                        <  1:875366/101‑1 (MQ=255)
                                                                                           gggAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCa                                                        <  3:304388/149‑1 (MQ=255)
                                                                                             gAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAgg                                                      <  1:429941/149‑1 (MQ=255)
                                                                                               aaTCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCa                                                                              >  2:540882/1‑123 (MQ=255)
                                                                                               aaTCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCa                                                                              <  1:540882/123‑1 (MQ=255)
                                                                                                 tCAAGGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGcaca                                                  >  2:1325593/1‑149 (MQ=255)
                                                                                                 tCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGc                                                                                                                         <  2:45079/78‑1 (MQ=255)
                                                                                                 tCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGc                                                                                                                         >  1:45079/1‑78 (MQ=255)
                                                                                                 tCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGcaca                                                  >  1:949028/1‑149 (MQ=255)
                                                                                                 tCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGcaca                                                  >  4:669457/1‑149 (MQ=255)
                                                                                                   aaaGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGcacaca                                                >  4:929011/1‑149 (MQ=255)
                                                                                                     agaagCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACacaac                                              <  3:1285334/149‑1 (MQ=255)
                                                                                                       aagCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAacac                                            <  2:949028/149‑1 (MQ=255)
                                                                                                        agCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGt                                                                                   <  1:516745/109‑1 (MQ=255)
                                                                                                        agCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGt                                                                                   >  2:516745/1‑109 (MQ=255)
                                                                                                           tCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACa                                        >  2:954900/1‑149 (MQ=255)
                                                                                                             cccACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTg                                                            >  3:1377903/1‑127 (MQ=255)
                                                                                                             cccACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTg                                                            <  4:1377903/127‑1 (MQ=255)
                                                                                                               cacaAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACg                                    >  4:588454/1‑149 (MQ=255)
                                                                                                                acaAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGt                                   <  1:435604/149‑1 (MQ=255)
                                                                                                                acaAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGt                                   >  2:166520/1‑149 (MQ=255)
                                                                                                                 caAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTc                                  <  1:1325593/149‑1 (MQ=255)
                                                                                                                 caAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTc                                  <  3:669457/149‑1 (MQ=255)
                                                                                                                  aaGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCt                                 >  3:929052/1‑149 (MQ=255)
                                                                                                                          tgctgcATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCg                         >  2:1291877/1‑149 (MQ=255)
                                                                                                                           gctgcATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGa                        >  3:477761/1‑149 (MQ=255)
                                                                                                                             tgcATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAacac                                            <  3:333529/127‑1 (MQ=255)
                                                                                                                             tgcATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAacac                                            >  4:333529/1‑127 (MQ=255)
                                                                                                                             tgcATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGc                      <  2:487242/149‑1 (MQ=255)
                                                                                                                              gcATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAg                                                       >  2:875217/1‑115 (MQ=255)
                                                                                                                              gcATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAg                                                       <  1:875217/115‑1 (MQ=255)
                                                                                                                               cATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCg                                                                    >  3:287439/1‑101 (MQ=255)
                                                                                                                               cATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCg                                                                    <  4:184015/101‑1 (MQ=255)
                                                                                                                               cATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCg                                                                    <  4:287439/101‑1 (MQ=255)
                                                                                                                               cATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCg                                                                    >  3:184015/1‑101 (MQ=255)
                                                                                                                               cATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTgcg                                                          <  3:549527/111‑1 (MQ=255)
                                                                                                                               cATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTgcg                                                          >  4:549527/1‑111 (MQ=255)
                                                                                                                                 tttGATGCCGATCGGCGTCACCAAGCCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGa                        <  1:950393/143‑1 (MQ=255)
                                                                                                                                 tttGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGa                  >  1:655475/1‑149 (MQ=255)
                                                                                                                                      tGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCa             <  3:929011/149‑1 (MQ=255)
                                                                                                                                       gCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAg            >  2:923304/1‑149 (MQ=255)
                                                                                                                                         cGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTc          <  3:588454/149‑1 (MQ=255)
                                                                                                                                          gATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCt         <  2:655475/149‑1 (MQ=255)
                                                                                                                                            tCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACacaaca                                             >  2:936743/1‑111 (MQ=255)
                                                                                                                                            tCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACacaaca                                             <  1:936743/111‑1 (MQ=255)
                                                                                                                                               gCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATAACGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGa                        >  1:712235/1‑129 (MQ=255)
                                                                                                                                                 gTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTgg                               >  3:1218925/1‑120 (MQ=255)
                                                                                                                                                 gTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTgg                               <  4:1218925/120‑1 (MQ=255)
                                                                                                                                                 gTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTaaaaa  >  1:64083/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAA  >  NC_000913/2517293‑2517586

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