Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,517,427 T→C 5.3% I136V (ATC→GTC)  yfeA ← putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,517,4270TC 5.3% 395.5 / 6.7 113I136V (ATC→GTC) yfeAputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (55/52);  new base C (3/3);  total (58/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 = 9.96e-01

Only 100 of 115 total aligned reads displayed.

TAGTAGAAAAGCATGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCA  >  NC_000913/2517279‑2517575
                                                                                                                                                    |                                                                                                                                                    
tagtagAAAAGCATGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGAt                                                                                                                                                      <  4:196282/149‑1 (MQ=255)
 agtagAAAAGCATGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATg                                                                                                                                                     >  4:346141/1‑149 (MQ=255)
  gtagAAAAGCATGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGc                                                                                                                                                    <  2:1004705/149‑1 (MQ=255)
          gCATGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGc                                                                                                                                            <  1:1121076/149‑1 (MQ=255)
          gCATGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGc                                                                                                                                            <  1:402591/149‑1 (MQ=255)
              gttgtAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCa                                                                                                                                        <  3:745943/149‑1 (MQ=255)
               ttgtAAATCAGCACAGCGGTGGAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATAGCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCac                                                                                                                                       >  2:280133/1‑149 (MQ=255)
               ttgtAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAAGCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCac                                                                                                                                       <  3:346141/149‑1 (MQ=255)
                          cacaGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGACGCCGATCGGCATCaaat                                                                                                                                     >  1:955456/1‑137 (MQ=255)
                          cacaGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCGCCCCACAAGATACATGCTGCATTTGACGCCGATCGGCATCaaat                                                                                                                                     <  2:955456/140‑4 (MQ=255)
                                   gAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGAc                                                                                                                   <  4:474056/149‑1 (MQ=255)
                                     aaaGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGAt                                                                                                                                                <  1:77089/118‑1 (MQ=255)
                                     aaaGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGAt                                                                                                                                                >  2:77089/1‑118 (MQ=255)
                                       aGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATc                                                                                                                                               >  1:1321810/1‑117 (MQ=255)
                                       aGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATc                                                                                                                                               <  2:1321810/117‑1 (MQ=255)
                                        ggCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGAc                                                                                                                   <  1:874092/144‑1 (MQ=255)
                                        ggCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGAc                                                                                                                   >  2:874092/1‑144 (MQ=255)
                                        ggCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGc                                                                                                              >  1:602145/1‑149 (MQ=255)
                                         gCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGtt                                                                                                                <  4:1115698/146‑1 (MQ=255)
                                         gCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGtt                                                                                                                >  3:1115698/1‑146 (MQ=255)
                                         gCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGACGCCGATCGGCATCaaat                                                                                                                                     >  2:1342492/1‑122 (MQ=255)
                                         gCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGACGCCGATCGGCATCaaat                                                                                                                                     <  1:1342492/125‑4 (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                                                                                                                                     >  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                                                                                                       >  3:381996/1‑148 (MQ=255)
                                                aaaTCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACat                                                                                                       <  4:381996/148‑1 (MQ=255)
                                                 aaTCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACatat                                                                                                     >  4:780195/1‑149 (MQ=255)
                                                 aaTCAACGACCATGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGAt                                                                                                                                                >  1:1446574/1‑106 (MQ=255)
                                                 aaTCAACGACCATGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGAt                                                                                                                                                <  2:1446574/106‑1 (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)
                                                              gAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGAATTTGATGCCGATCGGCGTCACCAAACCTAGCCAGAACAGACGTTGCCAGACAAATCGGCTGGTCAATc                                                                                        <  2:500167/149‑1 (MQ=255)
                                                                 aaTGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAAt                                                                                         <  1:775035/145‑1 (MQ=255)
                                                                       aTCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGc                                                                                                              >  4:488180/1‑118 (MQ=255)
                                                                       aTCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGc                                                                                                              <  3:488180/118‑1 (MQ=255)
                                                                          cGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAAc                                                                            <  2:602145/149‑1 (MQ=255)
                                                                           gCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGc                                                                                                                             <  1:955056/99‑1 (MQ=255)
                                                                             aTCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAAcgcg                                                                         >  4:145921/1‑149 (MQ=255)
                                                                              tCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATcc                                                                                       >  2:957779/1‑134 (MQ=255)
                                                                                          gTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGAc                                                                                                         >  1:397505/1‑104 (MQ=255)
                                                                                          gTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGAc                                                                                                         <  2:397505/104‑1 (MQ=255)
                                                                                          gTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCa                                                            >  4:474376/1‑149 (MQ=255)
                                                                                                 tCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATc                                                                                                                                               <  2:331066/59‑1 (MQ=255)
                                                                                                 tCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATc                                                                                                                                               >  1:331066/1‑59 (MQ=255)
                                                                                                  cTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAata                                                    >  2:429941/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)
                                                                                                           gAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAgg                                           <  3:802328/149‑1 (MQ=255)
                                                                                                             aaTCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCa                                                                   <  1:540882/123‑1 (MQ=255)
                                                                                                             aaTCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCa                                                                   >  2:540882/1‑123 (MQ=255)
                                                                                                               tCAAGGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGcaca                                       >  2:1325593/1‑149 (MQ=255)
                                                                                                               tCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGc                                                                                                              <  2:45079/78‑1 (MQ=255)
                                                                                                               tCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGc                                                                                                              >  1:45079/1‑78 (MQ=255)
                                                                                                               tCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGcaca                                       >  4:669457/1‑149 (MQ=255)
                                                                                                               tCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGcaca                                       >  1:949028/1‑149 (MQ=255)
                                                                                                                 aaaGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGcacaca                                     >  4:929011/1‑149 (MQ=255)
                                                                                                                 aaaGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGcacaca                                     >  1:487242/1‑149 (MQ=255)
                                                                                                                   agaagCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACacaac                                   <  3:1285334/149‑1 (MQ=255)
                                                                                                                     aagCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAacac                                 <  2:949028/149‑1 (MQ=255)
                                                                                                                      agCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGt                                                                        >  2:516745/1‑109 (MQ=255)
                                                                                                                      agCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGt                                                                        <  1:516745/109‑1 (MQ=255)
                                                                                                                         tCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACa                             >  2:954900/1‑149 (MQ=255)
                                                                                                                           cccACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTg                                                 <  4:1377903/127‑1 (MQ=255)
                                                                                                                           cccACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTg                                                 >  3:1377903/1‑127 (MQ=255)
                                                                                                                             cacaAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACg                         >  4:588454/1‑149 (MQ=255)
                                                                                                                              acaAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGt                        <  1:435604/149‑1 (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                                 >  4:333529/1‑127 (MQ=255)
                                                                                                                                           tgcATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAacac                                 <  3:333529/127‑1 (MQ=255)
                                                                                                                                           tgcATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGc           <  2:487242/149‑1 (MQ=255)
                                                                                                                                            gcATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAg                                            <  1:875217/115‑1 (MQ=255)
                                                                                                                                            gcATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAg                                            >  2:875217/1‑115 (MQ=255)
                                                                                                                                             cATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCg                                                         >  3:184015/1‑101 (MQ=255)
                                                                                                                                             cATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCg                                                         >  3:287439/1‑101 (MQ=255)
                                                                                                                                             cATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCg                                                         <  4:287439/101‑1 (MQ=255)
                                                                                                                                             cATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCg                                                         <  4:184015/101‑1 (MQ=255)
                                                                                                                                             cATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTgcg                                               >  4:549527/1‑111 (MQ=255)
                                                                                                                                             cATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTgcg                                               <  3:549527/111‑1 (MQ=255)
                                                                                                                                             cATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCAt         <  3:474376/149‑1 (MQ=255)
                                                                                                                                               tttGATGCCGATCGGCGTCACCAAGCCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGa             <  1:950393/143‑1 (MQ=255)
                                                                                                                                               tttGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGa             >  2:950393/1‑143 (MQ=255)
                                                                                                                                               tttGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGa       >  1:655475/1‑149 (MQ=255)
                                                                                                                                                    tGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCa  <  3:929011/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TAGTAGAAAAGCATGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCA  >  NC_000913/2517279‑2517575

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: