Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 2,189,469 C→T 5.6% E235E (GAG→GAA)  U712_10970 ← Tryptophan synthase alpha chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,189,4690CT 5.6% 237.0 / 2.1 71E235E (GAG→GAA) U712_10970Tryptophan synthase alpha chain
Reads supporting (aligned to +/- strand):  ref base C (32/35);  new base T (2/2);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

GCAAATCCAGTCACCTCATTTTAAACTATACAAGCCGCCAAACGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGA  >  CP006881/2189355‑2189605
                                                                                                                  |                                                                                                                                        
gCAAATCCAGTCACCTCATTTTAAACTATACAAGCCGCCAAACGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCa                                                                                                                  <  5:339094/139‑1 (MQ=255)
            aCCTCATTTTAAACTATACAAGCCGCCAAACGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCCGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACg                                                                                                      >  8:92030/1‑139 (MQ=255)
            aCCTCATTTTAAACTATACAAGCCGCCAAACGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACg                                                                                                      >  7:70751/1‑139 (MQ=255)
              cTCATTTTAAACTATACAAGCCGCCAAACGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGac                                                                                                    >  8:220511/1‑139 (MQ=255)
              cTCATTTTAAACTATACAAGCCGCCAAACGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGac                                                                                                    >  7:254981/1‑139 (MQ=255)
                           ataCAAGCCGCCAAACGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAAtt                                                                                       >  6:10562/1‑139 (MQ=255)
                           ataCAAGCCGCCAAACGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAAtt                                                                                       >  3:250189/1‑139 (MQ=255)
                           ataCAAGCCGCCAAACGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAAt                                                                                        >  4:49801/1‑138 (MQ=255)
                            taCAAGCCGCCAAACGCCCTTGCGTAATCTTCAGACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAAttt                                                                                      >  1:118632/1‑139 (MQ=255)
                            taCAAGCCGCCAAACGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAAttt                                                                                      >  1:49304/1‑139 (MQ=255)
                            taCAAGCCGCCAAACGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAAttt                                                                                      >  7:59692/1‑139 (MQ=255)
                            taCAAGCCGCCAAACGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATATCTCGTTTCAGTGCTGATGCGCCGGTCCTGTAGCTCTTCTATTTTTCTCAGGAGCGCACTACCCACTACGAAACCGTCCCAAAtt                                                                                       >  1:140969/1‑138 (MQ=255)
                                 gccgccAAACGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATttcatt                                                                                  >  6:176467/1‑138 (MQ=255)
                                  ccgccACACGCACTTGCGTAATCTTCAAAGTCCTGCAGCGCGTGATCGCTCGGTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATttcattca                                                                                <  6:349810/139‑1 (MQ=255)
                                  ccgccAAACGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATttcattca                                                                                <  5:274881/139‑1 (MQ=255)
                                  ccgccAAACGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATttcattca                                                                                <  2:172316/139‑1 (MQ=255)
                                  ccgccAAACGCCCTTGCGTAAGCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATttcattca                                                                                <  6:9166/139‑1 (MQ=255)
                                    gccAAACGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATc                                                                              >  3:93983/1‑139 (MQ=255)
                                    gccAAACGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATc                                                                              >  7:222339/1‑139 (MQ=255)
                                    gccAAACGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTAATTCATc                                                                              >  6:30402/1‑139 (MQ=255)
                                      ctaaCGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCtt                                                                            >  4:197611/3‑139 (MQ=255)
                                        aaCGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCa                                                                          <  4:250189/139‑1 (MQ=255)
                                           gCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTGTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGa                                                                       <  7:45025/139‑1 (MQ=255)
                                                   gTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTcg                                                               <  8:222339/139‑1 (MQ=255)
                                                    tAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTcgc                                                              >  3:270048/1‑139 (MQ=255)
                                                      aTCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTcgcg                                                             <  7:368311/138‑1 (MQ=255)
                                                      aTCTTCAAACTCCTGCAGCGCTTGATCTCTCGGTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTcgcg                                                             <  1:184085/138‑1 (MQ=255)
                                                       tcctCAAACTCCTGCAGCGCCTGATTTCTCGTTTCAGCGCTGATGAACCGGTCCTTTAATTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACTAATTTCATTCATCTTTATGACCTGTTCACGGt                                                           >  3:372857/4‑139 (MQ=255)
                                                       tCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAgcg                                                                                                                    >  4:115417/1‑82 (MQ=255)
                                                       tCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAgcg                                                                                                                    <  3:115417/82‑1 (MQ=255)
                                                       tCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGt                                                           >  3:194382/1‑139 (MQ=255)
                                                       tCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGt                                                           >  6:120646/1‑139 (MQ=255)
                                                                   tGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCg                                                >  4:203311/1‑138 (MQ=255)
                                                                   tGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCg                                                >  5:316912/1‑138 (MQ=255)
                                                                   tGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCg                                                >  5:298947/1‑138 (MQ=255)
                                                                     cAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTAGCGGTTTGATATACCGAAc                                             >  3:227534/1‑139 (MQ=255)
                                                                           ttGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACa                                       <  4:194382/139‑1 (MQ=255)
                                                                             gATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACa                                       <  4:270048/137‑1 (MQ=255)
                                                                             gATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACa                                       <  3:203311/137‑1 (MQ=255)
                                                                             gATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACa                                       <  8:70751/137‑1 (MQ=255)
                                                                                  ctcGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCg                                >  1:73684/1‑139 (MQ=255)
                                                                                   tcGTTTCCGTGCGGATGAGCCGGTCCTTTCGCTCTTCTTTTTTTCTGCCGAGCGCACTACCCACTACGACACCGACACAAATTTCATTCATCTTCATCCCCTGTCCGCGGTTTGATATACCGAACCCGACAGCAAc                                  <  2:241885/136‑1 (MQ=255)
                                                                                   tcGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAAc                                  <  6:298947/136‑1 (MQ=255)
                                                                                   tcGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAAc                                  <  6:316912/136‑1 (MQ=255)
                                                                                     gTTTCCGGGCGGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGaa                             <  2:49304/139‑1 (MQ=255)
                                                                                     gTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGaa                             <  2:73684/139‑1 (MQ=255)
                                                                                     gTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGaa                             <  5:10562/139‑1 (MQ=255)
                                                                                     gTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGaa                             <  3:197611/139‑1 (MQ=255)
                                                                                     gTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGaa                             <  7:220511/139‑1 (MQ=255)
                                                                                     gTTTCAGCGCTGATGAACCGGTCCTTTAATTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACTAATTTCATTCATCTTTATGACCTGTTCACGGTTTGATATACCGAACCCTACAGCAACCGGaa                             <  4:372857/139‑1 (MQ=255)
                                                                                     gTTCCAGTGCTGATGAGCCGGGCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTCCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGaa                             <  6:308629/139‑1 (MQ=255)
                                                                                        tCAGAGCTGCTGACCCGGCCCTTTAGCTCTTCTATTTTTCTGCCGCGCGCACTTCCCACTACGACACCGTCACAACTTTCATTCCTCTTCATGACCTGTTCGCGGTTTCATATACCGAACCCGACAGCAACCGGAACa                           <  5:26957/138‑1 (MQ=255)
                                                                                         cAGCGCTGATGAGCCGGTCCTTTAATTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACTAATTTCATTCATCTTTATGACCTGTTCACGGTTTGATATACCGAACCCTACAGCAACCGGaa                             >  8:127839/1‑135 (MQ=255)
                                                                                         cAGCGCTGATGAGCCGGTCCTTTAATTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACTAATTTCATTCATCTTTATGACCTGTTCACGGTTTGATATACCGAACCCTACAGCAACCGGaa                             <  7:127839/135‑1 (MQ=255)
                                                                                             gCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGCGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTg                     <  3:9826/139‑1 (MQ=255)
                                                                                             gCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCt                      >  4:188057/1‑138 (MQ=255)
                                                                                             gCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTg                     >  4:9826/1‑139 (MQ=255)
                                                                                              cTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACTACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTg                     <  5:30402/138‑1 (MQ=255)
                                                                                              cTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTg                     <  3:188057/138‑1 (MQ=255)
                                                                                               tgatgaGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATATTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTgag                   <  8:254981/139‑1 (MQ=255)
                                                                                                 atgaGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATttcattca                                                                                <  1:305280/76‑1 (MQ=255)
                                                                                                 atgaGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATttcattca                                                                                >  2:305280/1‑76 (MQ=255)
                                                                                                  tgaGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATc                <  4:93983/139‑1 (MQ=255)
                                                                                                      ccGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATttcattca                                                                                >  4:322453/1‑71 (MQ=255)
                                                                                                      ccGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATttcattca                                                                                <  3:322453/71‑1 (MQ=255)
                                                                                                      ccGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGCACCCGACAGCAACCGGAACAGTGCTGAGATCCTTc            >  1:280060/1‑139 (MQ=255)
                                                                                                      ccGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCtt             >  1:57806/1‑138 (MQ=255)
                                                                                                       cGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCa           <  8:306998/139‑1 (MQ=255)
                                                                                                       cGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCa           <  2:118632/139‑1 (MQ=255)
                                                                                                          tCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATAACGAACCCGACAGCAAACGGAACAGTGCTGAGATCCTTCACTg        >  2:192085/1‑139 (MQ=255)
                                                                                                                 gCTCTTATCTTTTTCTGAAGCGCTCACTTCACACTACGCCACCGTCAAAAATTTCATTCATCTTCCTGACCGGTTCGCGTGTTTATCTACCTAACCCGACAGCACCCGGAACAGTGCTGAGATCCTTCAATGTCCGGa  <  5:158928/138‑1 (MQ=255)
                                                                                                                  |                                                                                                                                        
GCAAATCCAGTCACCTCATTTTAAACTATACAAGCCGCCAAACGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGA  >  CP006881/2189355‑2189605

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑