Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 2,189,532 C→T 5.3% M214I (ATG→ATA)  U712_10970 ← Tryptophan synthase alpha chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,189,5320CT 5.3% 244.7 / 2.1 75M214I (ATG→ATA) U712_10970Tryptophan synthase alpha chain
Reads supporting (aligned to +/- strand):  ref base C (26/45);  new base T (2/2);  total (28/47)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AACGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCTAGAGAAGATACAC  >  CP006881/2189395‑2189668
                                                                                                                                         |                                                                                                                                        
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)
               tCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGt                                                                                                                          >  3:194382/1‑139 (MQ=255)
               tCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGt                                                                                                                          >  6:120646/1‑139 (MQ=255)
                           tGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCg                                                                                                               >  5:298947/1‑138 (MQ=255)
                           tGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCg                                                                                                               >  4:203311/1‑138 (MQ=255)
                           tGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCg                                                                                                               >  5:316912/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:316912/136‑1 (MQ=255)
                                           tcGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAAc                                                                                                 <  6:298947/136‑1 (MQ=255)
                                             gTTTCCGGGCGGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGaa                                                                                            <  2:49304/139‑1 (MQ=255)
                                             gTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGaa                                                                                            <  2:73684/139‑1 (MQ=255)
                                             gTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGaa                                                                                            <  7:220511/139‑1 (MQ=255)
                                             gTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGaa                                                                                            <  3:197611/139‑1 (MQ=255)
                                             gTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGaa                                                                                            <  5:10562/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                                                                                            <  7:127839/135‑1 (MQ=255)
                                                 cAGCGCTGATGAGCCGGTCCTTTAATTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACTAATTTCATTCATCTTTATGACCTGTTCACGGTTTGATATACCGAACCCTACAGCAACCGGaa                                                                                            >  8:127839/1‑135 (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)
                                                          tgaGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATc                                                                               <  4:93983/139‑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)
                                                                            cttctATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGaa                                                             >  6:272986/1‑139 (MQ=255)
                                                                             ttctATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGa                                                              <  5:176467/137‑1 (MQ=255)
                                                                              tctATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACg                                                           >  1:280137/1‑139 (MQ=255)
                                                                              tctATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACg                                                           >  3:74817/1‑139 (MQ=255)
                                                                              tctATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCACGGTTTGATATACCGAACCCGACAGCAACCGGTACAGTGCTGAGATCCTTCACTGTCCGCATGAACg                                                           >  1:329035/1‑139 (MQ=255)
                                                                              gctctttttCTGCCGAGCGCACTTCCACCGACGCCACCGGCACAAATTTCATTCATCTTCCTGACCTTGTCGCGGTGTGCTATACAGAACCCGACAGCACCCGGAACAGTGCTGAGATCCTTCACTGTCCGGATTAACg                                                           >  6:158928/5‑139 (MQ=255)
                                                                                          aCGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACgg                                                  <  4:74817/136‑1 (MQ=255)
                                                                                          aCGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACgg                                                  <  2:280137/136‑1 (MQ=255)
                                                                                          aCGAGCGCAATTCCAACTACGACACCGTAAAAAATTTCATTCATCTTCATGACCTGTTCGGGGTTTTATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATTAACGGGTACACGGAt                                                <  8:307201/138‑1 (MQ=255)
                                                                                               cgcACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACGGATGAAtt                                           >  2:11873/1‑138 (MQ=255)
                                                                                               cgcACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGAGGTTTGATATACCGAAACCGACAGCAACCGGAACAGTTCTGAGATCCTTCACTGTCCGTATGAACGGGTACACGGATGAAtt                                           >  7:307201/1‑138 (MQ=255)
                                                                                                         aCTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGaa                                                                                            >  5:248655/1‑79 (MQ=255)
                                                                                                         aCTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGaa                                                                                            <  6:248655/79‑1 (MQ=255)
                                                                                                               acacCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACGGATGAATTGAACTCATTGCGGacac                          >  5:6391/1‑139 (MQ=255)
                                                                                                                   cGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACa                                                                                          <  7:128886/71‑1 (MQ=255)
                                                                                                                   cGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACa                                                                                          >  8:128886/1‑71 (MQ=255)
                                                                                                                   cGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACGGAt                                                <  5:302545/113‑1 (MQ=255)
                                                                                                                   cGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACGGAt                                                >  6:302545/1‑113 (MQ=255)
                                                                                                                   cGTCAAAAATTTCATTCATCTTCATGACGTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGt                      <  5:272986/139‑1 (MQ=255)
                                                                                                                     tCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACa                                                                                                      <  1:197264/57‑1 (MQ=255)
                                                                                                                     tCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACa                                                                                                      >  2:197264/1‑57 (MQ=255)
                                                                                                                     tCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTca                    <  2:57806/139‑1 (MQ=255)
                                                                                                                     tCACAAATTTCATTCATCTTCATGACCTGGTCGCGGTTTGATATCCCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTca                    <  2:280060/139‑1 (MQ=255)
                                                                                                                     tCAAAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTTATATACCGAACCCGACAGCAACCGGAAGAGTGCTGAGATACTTCACTGTCCGGATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTca                    <  4:227534/139‑1 (MQ=255)
                                                                                                                            tttCATTCATCTTCATGACATGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCTag             <  3:49801/139‑1 (MQ=255)
                                                                                                                            gttctttcaGCTTCTTGCAGTGTTCACGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCTag             <  2:329035/138‑1 (MQ=255)
                                                                                                                            cttcattcaTCTTCACGACCTGTGCGCGGTTTCATATCCCGAACCCGCCAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACGGATGACTTGAACTCATTGCCGACACCGGTCACACCTag             <  4:272081/138‑1 (MQ=255)
                                                                                                                               cattcatCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCTagag           >  3:338336/1‑138 (MQ=255)
                                                                                                                                attcatCTTCATGACCTGTTCGCGGTTTGATATCCCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCTAgaga          <  5:120646/138‑1 (MQ=255)
                                                                                                                                    atCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCTAGAGAAGATa     <  1:11873/139‑1 (MQ=255)
                                                                                                                                        tCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAACAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCTAGAGAAGATacac  <  8:59692/138‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
AACGCCCTTGCGTAATCTTCAAACTCCTGCAGCGCTTGATCTCTCGTTTCAGTGCTGATGAGCCGGTCCTTTAGCTCTTCTATTTTTCTGACGAGCGCACTTCCCACTACGACACCGTCACAAATTTCATTCATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCTAGAGAAGATACAC  >  CP006881/2189395‑2189668

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: ‑