Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,100,813 A→G 5.1% F1566S (TTT→TCT)  U712_15925 ← Dimodular nonribosomal peptide synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,100,8130AG 5.1% 243.3 / 3.3 80F1566S (TTT→TCT) U712_15925Dimodular nonribosomal peptide synthase
Reads supporting (aligned to +/- strand):  ref base A (40/35);  new base G (2/2);  total (42/38)
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.96e-01

GCTGTGCATCTTTTAGCATAAAAGCAATTCGATCAGCCGGATAATCCGGATCAAGAGGGAGATAAGCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAACGGCAATGGCATCAGGACGTAAAGCAGCCTGCTTCT  >  CP006881/3100676‑3100948
                                                                                                                                         |                                                                                                                                       
gCTGTGCATCTTTTAGCATAAAAGCAATTCGATCAGCCGGATAATCCGGATCAAGAGGGAGATAAGCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACaaa                                                                                                                                        <  2:95905/139‑1 (MQ=255)
      cATCTTTTAGCATAAAAGCAATTCGATCAGCCGGATAATCCGGATCAAGAGGGAGATAAGCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTc                                                                                                                                  >  7:16389/1‑139 (MQ=255)
      cATCTTTTAGCATAAAAGCAATTCGATCAGCCGGATAATCCGGATCAAGAGGGAGATAAGCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTc                                                                                                                                  <  2:93435/139‑1 (MQ=255)
             tAGCATTAAAGCAATTCGATCAGCCGGATAATCCGGATCAAGAGGGAGATAAGCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGccc                                                                                                                            <  4:352091/138‑1 (MQ=255)
              aGCATAAAAGCAATTCGATCAGCCGGATAATCCGGATCAAGAGGGAGATAAGCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGa                                                                                                                          >  7:106487/1‑139 (MQ=255)
               gCATAAAAGCAATTCGATCAGCCGGATAATCCGGATCAAGAGGGAGATAAGCTGCCCCGGCCTTTCGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGa                                                                                                                          <  7:277942/138‑1 (MQ=255)
                  tcaaaGCAATTCGATCAGCCGGATAATCCGGATCAAGAGGGAGATAAGCTGAGCCGGCCTTTAGCACCTCCAATTGGCCGAGAGCGATTTCCAACGATCGGGGCATTGCCAATGCCACAAATTGCTGCGGCCCGACAcc                                                                                                                      <  2:296125/137‑1 (MQ=255)
                  tAAAAGCAATTCGATCAGCCGGATAATCCGGATCAAGAGGGAGATAAGCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGAGCCAATGCCACAAATTGCTCCGGCCCGACAcc                                                                                                                      <  6:110687/139‑1 (MQ=255)
                    aaaGCAATTCGATCAGCCGGATAATCCGGATCAAGAGGGAGATAAGCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGACACTGCCACAAATTGCTCCGGCCCGACACCtt                                                                                                                    >  4:404119/1‑139 (MQ=255)
                          aTTCGATCAGCCGGATAATCCGGATCAAGAGGGAGATAAGCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGacac                                                                                                                       >  8:352189/1‑130 (MQ=255)
                          aTTCGATCAGCCGGATAATCCGGATCAAGAGGGAGATAAGCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGacac                                                                                                                       <  7:352189/130‑1 (MQ=255)
                          aTTCGATCAGCCGGATAATCCGGATCAAGAGGGAGATAAGCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTa                                                                                                               >  3:56818/1‑138 (MQ=255)
                           ttCGATCAGCCGGATAATCCGGATCAAGAGGGAGATAAGCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAat                                                                                                             >  7:336892/1‑139 (MQ=255)
                                                      gagGGAGATAAGCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGctct                                                                                  >  6:26431/1‑139 (MQ=255)
                                                      gagGGAGATAAGCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGctct                                                                                  >  4:197648/1‑139 (MQ=255)
                                                       agGGAGATAAGCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAACTTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCtt                                                                                 <  6:10711/139‑1 (MQ=255)
                                                       agGGAGATAAGCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCtt                                                                                 <  8:10430/139‑1 (MQ=255)
                                                            gaTAAGCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGCTCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTttcatt                                                                            <  8:246927/139‑1 (MQ=255)
                                                            gaTAAGCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTttcatt                                                                            <  4:56818/139‑1 (MQ=255)
                                                            gaTAAGCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTttcatt                                                                            <  8:73026/139‑1 (MQ=255)
                                                             aTAAGCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTg                                                                                      <  1:315342/128‑1 (MQ=255)
                                                             aTAAGCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTg                                                                                      >  2:315342/1‑128 (MQ=255)
                                                                 ggcgCCCCGGCCTTAAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGATTGCTCTTTCATTCAGct                                                                       <  7:28092/136‑1 (MQ=255)
                                                                 gCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGct                                                                       <  5:289474/139‑1 (MQ=255)
                                                                 cgctccccGGCCTTTAGCACCGCCAATAGGCCGCCAGGCATTTCCAACGCTCGGGGCAGTGACAATGCCACAAATTGCTCCGGCCAGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGct                                                                       <  5:265820/135‑1 (MQ=255)
                                                                    gCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTg                                                                    >  1:176615/1‑139 (MQ=255)
                                                                    gCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTg                                                                    >  7:203074/1‑139 (MQ=255)
                                                                     ccccGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTg                                                                    >  3:357315/1‑138 (MQ=255)
                                                                     ccccGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTg                                                                    >  8:234034/1‑138 (MQ=255)
                                                                       ccGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCtt                                                                                 >  6:253847/1‑123 (MQ=255)
                                                                       ccGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCtt                                                                                 <  5:253847/123‑1 (MQ=255)
                                                                          gCCTTTAGCCCCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGCCACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGc                                                              >  3:315607/1‑139 (MQ=255)
                                                                          gCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGctc                                                                                   >  6:199920/1‑118 (MQ=255)
                                                                          gCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGctc                                                                                   <  5:199920/118‑1 (MQ=255)
                                                                          gCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGc                                                              >  2:214061/1‑139 (MQ=255)
                                                                           ccTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCAtctgt                                                             <  7:234034/139‑5 (MQ=255)
                                                                               tAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAg                                                                                              <  2:356839/102‑1 (MQ=255)
                                                                               tAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAg                                                                                              >  1:356839/1‑102 (MQ=255)
                                                                               tAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGt                                                         >  4:124154/1‑139 (MQ=255)
                                                                                     cGCCAATAGGCCGACAGCCATTTCCACCGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGa                                                   >  5:88568/1‑139 (MQ=255)
                                                                                     cGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGa                                                   <  8:106487/139‑1 (MQ=255)
                                                                                       ccAATAGGCCGCCAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAGATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGAtt                                                 >  6:364617/1‑139 (MQ=255)
                                                                                       ccAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAGATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGAtt                                                 >  6:370407/1‑139 (MQ=255)
                                                                                               ccGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGctc                                                                                   <  4:129403/97‑1 (MQ=255)
                                                                                               ccGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGctc                                                                                   >  3:129403/1‑97 (MQ=255)
                                                                                               ccGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCatat                                          >  7:36546/1‑138 (MQ=255)
                                                                                               ccGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCAtata                                         >  8:323724/1‑139 (MQ=255)
                                                                                                  aCAGCCATTTCCAACGATCTGGGCAGTGCCAATGCCACAAATAGCTCCGGCCCGACACCTTCACTAATCATAATTCCGCCAAGACGGTTTGCTCTTAAATTAAGATCAAAATACCTCAGATCCTCATTTTCAtataaa                                       >  6:293244/1‑136 (MQ=255)
                                                                                                  aCAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACaa                                      >  1:258681/1‑139 (MQ=255)
                                                                                                  aCAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACaa                                      >  3:239442/1‑139 (MQ=255)
                                                                                                   cAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAAc                                     >  2:117062/1‑139 (MQ=255)
                                                                                                   cAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAAc                                     >  8:181281/1‑139 (MQ=255)
                                                                                                      ccATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGac                                                                                                                         >  1:42977/1‑52 (MQ=255)
                                                                                                      ccATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGac                                                                                                                         <  2:42977/52‑1 (MQ=255)
                                                                                                       cATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAACGGCa                                 <  5:26431/139‑1 (MQ=255)
                                                                                                       cATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAACGGCa                                 <  3:197648/139‑1 (MQ=255)
                                                                                                         tttCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAACGGCaa                                <  3:394822/138‑1 (MQ=255)
                                                                                                         tttCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAACGGCaa                                >  4:394822/1‑138 (MQ=255)
                                                                                                                cGATCGGGGCAGTGCCAATTCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGCGCAATCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAACGGCAATGGCATCa                        <  6:240691/139‑1 (MQ=255)
                                                                                                                cGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAACGGCAATGGCAt                          >  3:147530/1‑137 (MQ=255)
                                                                                                                cGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAACGGCAATGGCAt                          >  8:217132/1‑137 (MQ=255)
                                                                                                                cGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAACGGCAATGGCATCa                        <  4:357315/139‑1 (MQ=255)
                                                                                                                cGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAACGGCAATGGCATCa                        <  8:336892/139‑1 (MQ=255)
                                                                                                                cGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAACGGCAATGGCATCa                        <  2:176615/139‑1 (MQ=255)
                                                                                                                cGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAACGGCAATGGCATCa                        <  7:323724/139‑1 (MQ=255)
                                                                                                                cGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAACGGCAATGGCATCa                        <  8:203074/139‑1 (MQ=255)
                                                                                                                     ggggCAGTGCCAATGCCACAGATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAACGGCAATGGCATCAGGAc                    <  5:370407/138‑1 (MQ=255)
                                                                                                                     ggggCAGTGCCAATGCCACAGATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAACGGCAATGGCATCAGGAc                    <  5:364617/138‑1 (MQ=255)
                                                                                                                            tGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAACGGCAATGGCATCa                        >  3:272991/1‑127 (MQ=255)
                                                                                                                            tGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAACGGCAATGGCATCa                        <  4:272991/127‑1 (MQ=255)
                                                                                                                            tGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAACGGCAATGGCATCAGGACGTAAag              >  5:47171/1‑137 (MQ=255)
                                                                                                                             gCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGct                                                                       >  8:257891/1‑79 (MQ=255)
                                                                                                                             gCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGct                                                                       <  7:257891/79‑1 (MQ=255)
                                                                                                                                 aTGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAACGGCAATGGCATCAGGACGTAAAGCAGCCTg       >  5:286445/1‑139 (MQ=255)
                                                                                                                                 aTGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAACGGCAATGGCATCAGGACGTAAAGCAGCCTg       >  3:402463/1‑139 (MQ=255)
                                                                                                                                     cacaTATTGCTCCGGACCGACACCTTCACTAATGATCATTCGGGCAATCCGGTTTTCTCTTTCATTCAGATCTGCATAGCTCAGTTCCTGATTTTCTTATAGAACGGCATTGGGATCAGGACGTAAAGCAGCCTGCtt  <  8:379231/139‑2 (MQ=255)
                                                                                                                                      acaAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTttcatt                                                                            <  3:247800/65‑1 (MQ=255)
                                                                                                                                      acaAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTttcatt                                                                            >  4:247800/1‑65 (MQ=255)
                                                                                                                                      acaAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCa                                                                  >  4:67075/1‑75 (MQ=255)
                                                                                                                                      acaAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCa                                                                  <  3:67075/75‑1 (MQ=255)
                                                                                                                                      acaAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAACGGCAATGGCATCAGGACGTAAAGCAGCCTGCTtct  >  4:79070/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                       
GCTGTGCATCTTTTAGCATAAAAGCAATTCGATCAGCCGGATAATCCGGATCAAGAGGGAGATAAGCTGCCCCGGCCTTTAGCACCGCCAATAGGCCGACAGCCATTTCCAACGATCGGGGCAGTGCCAATGCCACAAATTGCTCCGGCCCGACACCTTCACTAATCATCATTCGGGCAAGCCGGTTTGCTCTTTCATTCAGCTCTGCATAGCTCAGTTCCTGATTTTCATATACAACGGCAATGGCATCAGGACGTAAAGCAGCCTGCTTCT  >  CP006881/3100676‑3100948

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