Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,172,226 T→C A100A (GCA→GCG U712_10860 ← 3‑methyl‑2‑oxobutanoate hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,172,2260TC98.9% 283.2 / ‑4.7 90A100A (GCA→GCGU712_108603‑methyl‑2‑oxobutanoate hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (33/55);  total (34/56)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

ACTGCGGTGTCAAACCTAAGTGACTGACTACTGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAACTGCTGCATTTTTCAGCGTATCCTCCTTGGACAGGTGATAAGACATAAACGGCATGTCTGTCACAATAAAAGTATTCTGCGCGCCCCTTTTAACAGCTTTTGTATGATGGATCATGTCCGCAACTGT  >  CP006881/2172099‑2172351
                                                                                                                               |                                                                                                                             
acggCGGTGGCAAACCTAAGTGACTGACTACTGGAATGCCTCCAACCGGCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCAttt                                                                                                                      <  2:430515/134‑1 (MQ=255)
      gtgtCAAACCTAAGTGACTGACTACTGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAATGATAGCCGCTGCATTTTTCAGCg                                                                                                               >  4:409621/1‑138 (MQ=255)
      gtgtCAAACCTAAGTGACTGACTACTGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCGGCATTTTTCAGCg                                                                                                               >  3:17057/1‑138 (MQ=255)
        gtCAAACCTAAGTGACTGACTACTGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCACGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTAt                                                                                                            <  5:228531/139‑1 (MQ=255)
        gtCAAACCTAAGTGACTGACTACTGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTAt                                                                                                            <  3:202979/139‑1 (MQ=255)
        gtCAAACCTAAGTGACTGACTACTGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTAt                                                                                                            <  8:381917/139‑1 (MQ=255)
        gtCAAACCTAAGTGACTGACTACTGGAATGCCTCCAAGCGTCAATGCGCGAATGGAGTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTAt                                                                                                            <  8:283354/139‑1 (MQ=255)
        gtCAAACCTAAGTGACTGACTACTGGAAGGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCGTCTCCGCCCTCAAGCTTCAGTGAGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTAt                                                                                                            <  4:285883/139‑1 (MQ=255)
        ggcAAACCTAAGTGAATGACTACTGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTAt                                                                                                            <  2:221575/137‑1 (MQ=255)
         tCAAACCTAAGTGACTGACTACTGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTAt                                                                                                            <  3:68939/138‑1 (MQ=255)
         tCAAACCTAAGTGACTGACTACTGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTAt                                                                                                            <  5:152603/138‑1 (MQ=255)
         tCAAACCTAAGTGACTGACTACTGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTAt                                                                                                            <  7:309924/138‑1 (MQ=255)
         tCAAACCTAAGTGACTGACTACTGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTAt                                                                                                            <  2:262093/138‑1 (MQ=255)
         tCAAACCTAAGTGACTGACTACTGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTAt                                                                                                            <  5:45378/138‑1 (MQ=255)
         tCAAACCTAAGTGACTGACTACTGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTAt                                                                                                            <  8:256610/138‑1 (MQ=255)
         tCAAACCTAAGTGACTGACTACTGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTAt                                                                                                            <  8:305210/138‑1 (MQ=255)
           aaaCCTAAGTGACTGACTTCTGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCAGCTTTCCTGAACGATAGCCTCTGCATTTTTCAGCGTAtctt                                                                                                         <  5:290310/139‑3 (MQ=255)
                   gTGACTGACTACTGGAATGCCTACAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTAt                                                                                                            <  4:395785/128‑1 (MQ=255)
                   gTGACTGACTACTGGAATGCCTACAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTAt                                                                                                            >  3:395785/1‑128 (MQ=255)
                      actgactACTGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAgg                                                                                              >  4:104008/1‑139 (MQ=255)
                      actgactACTGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAg                                                                                               <  1:435829/138‑1 (MQ=255)
                        tgactACTGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGt                                                                                             >  7:311782/1‑138 (MQ=255)
                          actactGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGAt                                                                                          <  8:297431/139‑1 (MQ=255)
                          actactGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGAt                                                                                          <  8:157870/139‑1 (MQ=255)
                          actactGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGAt                                                                                          <  7:203989/139‑1 (MQ=255)
                          actactGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGAt                                                                                          <  6:423511/139‑1 (MQ=255)
                          actactGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGAt                                                                                          <  6:92858/139‑1 (MQ=255)
                          actactGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCCTCCTTAGACAGGTGAt                                                                                          <  5:9788/139‑1 (MQ=255)
                          actactGGAATGCCTCCAAGCGGCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGAt                                                                                          <  4:138783/139‑1 (MQ=255)
                                  aaggCCACCAAGCGTCAATGCGCGAATGTATTCAAACCCCCCTTCTCCGCCCTCACGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACTATCGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAACACATa                                                                                  <  4:17057/136‑1 (MQ=255)
                                  aaTGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATa                                                                                  <  3:404160/139‑1 (MQ=255)
                                    tGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCg                                                                                                               <  1:228847/108‑1 (MQ=255)
                                    tGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCg                                                                                                               >  2:228847/1‑108 (MQ=255)
                                    tGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATaca                                                                                <  8:195321/139‑3 (MQ=255)
                                    tGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATaaa                                                                                >  7:195321/1‑139 (MQ=255)
                                    tGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATaaa                                                                                <  8:139194/139‑1 (MQ=255)
                                      ccTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTAt                                                                                                            <  6:432535/109‑1 (MQ=255)
                                      ccTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTAt                                                                                                            >  5:432535/1‑109 (MQ=255)
                                      ccTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGa                                                                                                  <  3:57146/119‑1 (MQ=255)
                                      ccTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGa                                                                                                  >  4:57146/1‑119 (MQ=255)
                                      ccTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACg                                                                              >  2:110044/1‑139 (MQ=255)
                                      ccTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAACTGCTGCATTTTTCAGCGTATCCTCCTTGGCCAGGTGATAAGACATAAACg                                                                              >  1:270325/1‑139 (MQ=255)
                                      ccTCAAAGCTTCCATGCGCTAATGTCTTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCGTTCCTGCACGCTAGCCGCTGCATTTTTCGGCGTCGCTTCCTTAGACAGGGGAGAAGACCTAAACg                                                                              >  8:198812/1‑139 (MQ=255)
                                          cacgCGTCAATGCGCGAATGGATTCAACCACGCCTCCTCCGCCCTCAAGCTTCAGGGCGTCAGCTCCGCTTACCTGAACGATAGCCGCGGCATTTTTCAGCGTATCTTCCTTAGACAGGTCATAAGACATAAACGGc                                                                            <  1:178758/134‑1 (MQ=255)
                                           aaGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCAta                                                                         <  8:311782/139‑2 (MQ=255)
                                            aGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTAt                                                                                                            <  7:148264/103‑1 (MQ=255)
                                            aGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTAt                                                                                                            >  8:148264/1‑103 (MQ=255)
                                            aGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGa                                                                                                  >  7:303940/1‑113 (MQ=255)
                                            aGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGa                                                                                                  <  8:303940/113‑1 (MQ=255)
                                                tCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACa                                                                                                >  4:313396/1‑111 (MQ=255)
                                                tCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACa                                                                                                <  3:313396/111‑1 (MQ=255)
                                                tCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCAtatctgt                                                                    <  1:110044/139‑1 (MQ=255)
                                                tCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCAtatctgt                                                                    <  4:225639/139‑1 (MQ=255)
                                                 cAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCGTTCCTGAACGATAGCCGCTGCATTTTTCAGAGTATCTTCCTTTGACAGGTGATAAGACATAAACGGCAtatatgt                                                                    >  4:53643/1‑138 (MQ=255)
                                                      cgcgAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACaa                                                               <  3:448413/138‑1 (MQ=255)
                                                        cgAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGCACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGCCAGGTGATAAGACATAACCGGCATATCTGTCACAATa                                                             >  5:181797/1‑138 (MQ=255)
                                                            tGGATTCAAACACGCCTTCTCCGGCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTATCAGCGAATCTTCCTTAGACAGGTGATAAGACATTATCGGCATATCTGTCAGAATAAAGGt                                                        >  7:149158/1‑139 (MQ=255)
                                                                ttCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTc                                                    >  8:368801/1‑139 (MQ=255)
                                                                     acacGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGgcgc                                               >  3:401696/1‑139 (MQ=255)
                                                                      cacGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTCGACAGGTGATAAGACATAAACGGCATATCTGTCACATTAAAGGTATTCGgcgca                                              <  3:53643/139‑2 (MQ=255)
                                                                      cacGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGgcgca                                              <  7:127998/139‑2 (MQ=255)
                                                                      cacGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGgcgca                                              <  5:372872/139‑2 (MQ=255)
                                                                         gCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGgcgcac                                             <  8:205909/137‑3 (MQ=255)
                                                                         gCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGgcgcac                                             >  7:205909/1‑135 (MQ=255)
                                                                               tcCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATaa                                                            >  4:81741/1‑116 (MQ=255)
                                                                               tcCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATaa                                                            <  3:81741/116‑1 (MQ=255)
                                                                                   cccTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGAt                                                                                          <  1:127374/82‑1 (MQ=255)
                                                                                   cccTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGAt                                                                                          >  2:127374/1‑82 (MQ=255)
                                                                                   cccTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGGCGCACCCCTTTTAAcgg                                 >  5:328960/1‑137 (MQ=255)
                                                                                   cccTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGGCGCACCCCTTTTAAcgg                                 >  7:319107/1‑137 (MQ=255)
                                                                                     cTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGGCGCACCCCTTTTAACGGCt                               >  3:76920/1‑139 (MQ=255)
                                                                                        aaGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGGCGCACCCCTTTTAACGGCtttt                            <  6:171692/139‑1 (MQ=255)
                                                                                              cAGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGGCGCACCCCTTTTAACGGCTTTTGTatga                      >  3:204112/1‑139 (MQ=255)
                                                                                               aGTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGGCGCACCCCTTTTAACGGCTTTTGTatga                      <  3:104008/138‑1 (MQ=255)
                                                                                                gTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGGCGCACCCCTTTTAACGGCtttt                            <  2:38903/131‑1 (MQ=255)
                                                                                                gTGCGTCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGGCGCACCCCTTTTAACGGCtttt                            >  1:38903/1‑131 (MQ=255)
                                                                                                    ggcAGCTCCGCTTTCCTGAACGATAGCGGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGGCGCACCCCTTTTAACGGCTTTTGTATGATGGATc                <  3:66028/137‑1 (MQ=255)
                                                                                                    gTCAGCTCCGCTTTCCTGAATGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGGCGCACCCCTTTTAACGGCTTTTGTATGATGGATc                <  3:409621/139‑1 (MQ=255)
                                                                                                     tCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATaaa                                                                                <  1:100934/74‑1 (MQ=255)
                                                                                                     tCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATaaa                                                                                >  2:100934/1‑74 (MQ=255)
                                                                                                     tCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATa                                                             <  8:421639/93‑1 (MQ=255)
                                                                                                     tCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATa                                                             >  7:421639/1‑93 (MQ=255)
                                                                                                     tCAGCTCCGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGGCGCACCCCTTTTAACGGCTTTTGTATGATGGATCa               >  2:312644/1‑139 (MQ=255)
                                                                                                           ccGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACATGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGGCGCACCCCTTTTAACGGCTTTTGTATGATGGATCATGTCCg         >  2:142799/1‑139 (MQ=255)
                                                                                                            cGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGGCGCACCCCTTTTAACGGCTTTTGTATGATGGATCATGTCCGc        <  7:368801/139‑1 (MQ=255)
                                                                                                            cGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGGCGCACCCCTTTTAACGGCTTTTGTATGATGGATCATGTCCGc        <  8:319107/139‑1 (MQ=255)
                                                                                                            cGCTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGGCGCACCCCTTTTAACGGCTTTTGTATGATGGATCATGTCCGc        <  4:295495/139‑1 (MQ=255)
                                                                                                              cTTTCCTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGGCGCACCCCTTTTAACGGCTTTTGTATGATGGATCATGTCCGCaa      >  8:376015/1‑139 (MQ=255)
                                                                                                                   cTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGGCGCACCCCTTTTAACGGCTTTTGTATGATGGATCATGTCCGCAACTGt  >  6:23278/1‑138 (MQ=255)
                                                                                                                   cTGAACGATAGCCGCTGCATTTTTCAGCGTATCTTCCTTAGACAGGTGATAAGACATAAACGGCATATCTGTCACAATAAAGGTATTCGGCGCACCCCTTTTAACGGCTTTTGTATGATGGATCATGTCCGCAACTGt  <  4:76920/138‑1 (MQ=255)
                                                                                                                               |                                                                                                                             
ACTGCGGTGTCAAACCTAAGTGACTGACTACTGGAATGCCTCCAAGCGTCAATGCGCGAATGGATTCAAACACGCCTTCTCCGCCCTCAAGCTTCAGTGCGTCAGCTCCGCTTTCCTGAACGATAACTGCTGCATTTTTCAGCGTATCCTCCTTGGACAGGTGATAAGACATAAACGGCATGTCTGTCACAATAAAAGTATTCTGCGCGCCCCTTTTAACAGCTTTTGTATGATGGATCATGTCCGCAACTGT  >  CP006881/2172099‑2172351

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: