Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 2,188,141 C→G 5.3% P29P (CCG→CCC U712_10960 ← Prephenate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,188,1410CG 5.3% 255.0 / 4.7 76P29P (CCG→CCCU712_10960Prephenate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (34/38);  new base G (2/2);  total (36/40)
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.38e-01

ATTACTGTGGCTGCTTTTTTTGCACCGCTAATAAACGAATCAGCACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACC  >  CP006881/2188019‑2188273
                                                                                                                          |                                                                                                                                    
aTTACTGTGGCTGCTTTTTTTGCACCGCTATTTAACGAATCAGCACGCTCGTGTATCACTCCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGAGTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGa                                                                                                                      >  7:142562/1‑139 (MQ=255)
aTTACTGTGGCTGCTTTTTTTGAACCGCTTATAAACGAATCAGCACCCTCGTATATCACACCTAATTTTAATGCGGCTACTACCTGTTCATAAGAGATGTCTATTCAGATAATCCGTTTGCCCGGATGCTCTTTTTTGa                                                                                                                      >  8:149781/1‑139 (MQ=255)
     tgtgGCTGCTTTTTTTGAAAATCTAATAAACGAATCCGCACGCGCGTCTATCACACCGAATTTTAATGCGGCGACTAAATGTTCATAAGATATGTAGATTCCGATAATCCGTTTTACCGTATGCTTTTTTTTTATCGCt                                                                                                                 >  4:153187/1‑139 (MQ=255)
         gctgctTTTTTTGCACCGCTAATAAACGAATCAGCACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAgcg                                                                                                             >  5:245491/1‑139 (MQ=255)
         gctgctTTTTTTGCACCGCTAATAAACGAATCAGCACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAgcg                                                                                                             >  8:316626/1‑139 (MQ=255)
         gctgctTTTTTTGCACCGCTAATAAACGAATCAGCACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAgcg                                                                                                             >  8:177132/1‑139 (MQ=255)
         gctgctTTTTTTGCACCGCTAATAAACGAATCAGCACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAgcg                                                                                                             >  2:37079/1‑139 (MQ=255)
         gctgcgTTTTTTGCACCGCTAAAAAACGAATCAGCACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCAGTTCATCAGAGATGTCGATTCCGATAATACCTTTGCCCGGATGCTCTTTTTTGATCGCTTgcg                                                                                                             >  2:225485/1‑139 (MQ=255)
          ctgctTTTTTTGCACCGCTAATAAACGAATCAGCACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCgg                                                                                                                                    >  1:180669/1‑115 (MQ=255)
          ctgctTTTTTTGCACCGCTAATAAACGAATCAGCACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCgg                                                                                                                                    <  2:180669/115‑1 (MQ=255)
                 ttttGCACCGCTAATAAACGAATCAGCACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGa                                                                                                      <  5:283194/138‑1 (MQ=255)
                         cGCTAATAAACGAATCAGCACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTAt                                                                                             <  4:218419/139‑1 (MQ=255)
                         cGCTAATAAACGAATCAGCACGATCGTCTATCACGCCTAATTTTAATGCCGCTACCGCCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATGATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTa                                                                                              >  4:39701/1‑138 (MQ=255)
                                     aaTCAGCACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGAcc                                                                                 >  7:52838/1‑139 (MQ=255)
                                     aaTCAGCACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGAcc                                                                                 >  4:104402/1‑139 (MQ=255)
                                      aTCAGCACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCg                                                                                >  2:309110/1‑139 (MQ=255)
                                      aTCAGCACGCTCGTCTATCACAACTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCTATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCg                                                                                <  1:225485/139‑1 (MQ=255)
                                         aGCACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGc                                                                              <  5:167697/138‑1 (MQ=255)
                                         aGCACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGc                                                                              <  4:269724/138‑1 (MQ=255)
                                         aGCACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGc                                                                              <  7:149781/138‑1 (MQ=255)
                                         aGCACGCTCGTCTATCACACCTAATTATAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGc                                                                              <  2:278555/138‑1 (MQ=255)
                                           cACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAg                                                                           <  7:313766/139‑1 (MQ=255)
                                           cACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAg                                                                           <  7:174335/139‑1 (MQ=255)
                                           cACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAg                                                                           <  7:177132/139‑1 (MQ=255)
                                           cACGCTCGTCTATCACAACTAATTTTAATGCGTCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCACTTATTCTCTTTTTTGATCGCTATAGCAATCTAACCGCCTATCAATCCGAGACCGGagag                                                                           <  1:23238/139‑1 (MQ=255)
                                            aCGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGc                                                                          <  3:104402/139‑1 (MQ=255)
                                                tCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGAGTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTTATCGCTAGCGCAATCGCACCGCCTAGCAATCCGAGACCGGCGAGCAat                                                                       >  6:321978/1‑138 (MQ=255)
                                                  gTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATAtt                                                                    <  3:199711/139‑1 (MQ=255)
                                                       tCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGt                                                                  <  6:245491/136‑1 (MQ=255)
                                                            ccTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCa                                                          >  7:286406/1‑139 (MQ=255)
                                                               aaTTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTc                                                       >  2:296114/1‑139 (MQ=255)
                                                                aTTTTAATGCGGCTACTTCGTGTTCATCTGAGATGTCGATTCCTATCATCCGTTGGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCAtt                                                        <  5:257589/137‑1 (MQ=255)
                                                                aTTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCAtt                                                        <  5:162115/137‑1 (MQ=255)
                                                                aTTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCAtt                                                        <  3:180084/137‑1 (MQ=255)
                                                                     aaTGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGATAGCGCAATCGAACCGCCTATCAAACCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTc                                                  >  4:12063/1‑138 (MQ=255)
                                                                      aTGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCa                                                 <  1:37079/138‑1 (MQ=255)
                                                                      aTGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCAt                                                <  2:291207/139‑1 (MQ=255)
                                                                           gCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCTAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGtc                                           <  1:309112/139‑1 (MQ=255)
                                                                           gCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGtc                                           <  8:286406/139‑1 (MQ=255)
                                                                           gCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGtc                                           <  1:309110/139‑1 (MQ=255)
                                                                           gCTACTACCTGTTCATCAGAGATGTCGATTCCAATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGt                                            >  8:92608/1‑138 (MQ=255)
                                                                           gCTAATACCTGTTCATAATAGATGTCGATTCCTATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTATCTCAATATAACCGCCTATCAATCCGAGACCGGCGATCAATATTGTATCTTTCATTCGATTCATACgga                                           <  3:153187/139‑3 (MQ=255)
                                                                            ctaccgccTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATGATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGtca                                          >  5:55903/7‑139 (MQ=255)
                                                                                cgccTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCGGGATGATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCAc                                      <  3:39701/137‑1 (MQ=255)
                                                                                  ccTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGTCATCAc                                      <  3:92894/137‑1 (MQ=255)
                                                                                      ttCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGc                                                                                                                                >  7:37866/1‑43 (MQ=255)
                                                                                      ttCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGc                                                                                                                                <  8:37866/43‑1 (MQ=255)
                                                                                      ttCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAg                                                                           >  5:162871/1‑96 (MQ=255)
                                                                                      ttCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAg                                                                           <  6:162871/96‑1 (MQ=255)
                                                                                      ttCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCAtt                                                        >  6:111499/1‑115 (MQ=255)
                                                                                      ttCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCAtt                                                        <  5:111499/115‑1 (MQ=255)
                                                                                      ttCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGAtt                                                   >  6:153666/1‑120 (MQ=255)
                                                                                      ttCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGAtt                                                   <  5:153666/120‑1 (MQ=255)
                                                                                      ttCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCTTACGTCATCACCTGGtt                                >  8:165007/1‑139 (MQ=255)
                                                                                      ttCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGTCATCACCTGGtt                                >  1:23685/1‑139 (MQ=255)
                                                                                      ttCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGTCATCACCTGGtt                                >  5:173446/1‑139 (MQ=255)
                                                                                      ttCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGTCATCACCAGGtt                                >  1:115423/1‑139 (MQ=255)
                                                                                      ttCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGTCATCAACTGGtt                                >  4:333561/1‑139 (MQ=255)
                                                                                      ttCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCCAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGTAATCACCTGGtt                                >  5:37460/1‑139 (MQ=255)
                                                                                          tcaGAGATGTCGATTCCGATAATCCGTTTGCCGGGATGATTTTTTTTGATGGCTAGGGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTTGATTCATACGTCATCACCTGGTTtata                            <  6:55903/139‑1 (MQ=255)
                                                                                          tcaGAGATGTCGATTCCAATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGTCATCACCTGGTTtata                            <  7:92608/139‑1 (MQ=255)
                                                                                          tcaGAGATGTCGATTACGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGTCATCACCTGGTTtata                            <  8:75096/139‑1 (MQ=255)
                                                                                           caGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCg                                                                                <  4:276975/86‑1 (MQ=255)
                                                                                           caGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCg                                                                                >  3:276975/1‑86 (MQ=255)
                                                                                                aTGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGTCATCACCTGGTTTATAGACTa                       >  7:72/1‑138 (MQ=255)
                                                                                                 tggcGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGTCATCACCTGGTTTATAGACTaaa                     <  2:99149/136‑1 (MQ=255)
                                                                                                 tGTCGATTCCGAGAATCCCTTTGCCCGGATGCTCGTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATACTGTATCTTTCATTCGAATCATACGTCATCACCTGGTTTATAGACTaaa                     <  8:142562/139‑1 (MQ=255)
                                                                                                   tCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAg                                                                           <  3:203419/83‑1 (MQ=255)
                                                                                                   tCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAg                                                                           >  4:203419/1‑83 (MQ=255)
                                                                                                      aTTACGATAATTCGCTTGACCGGATGATCTTTTTTGATCGCTAGCGCAATCGTACAGCCTATCAATCCGAGACCGGCGAGCTATATTGTAGCTTTCATTCGATTCATACGTCATCACCTGGGTTATAGACTAAAGAGGt                <  7:225278/139‑1 (MQ=255)
                                                                                                        tCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGCCGAGCAATATTGTATCTTTCATTCGATTCATACGTCATCACCTGGTTAATAGACTAAAGAggata              >  8:225278/1‑136 (MQ=255)
                                                                                                          cGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGTCATCACCTGGttt                               >  3:230120/1‑120 (MQ=255)
                                                                                                          cGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGTCATCACCTGGttt                               <  4:230120/120‑1 (MQ=255)
                                                                                                                   gTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGACTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCAc   <  5:204531/139‑1 (MQ=255)
                                                                                                                    tttGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCAcc  >  4:88357/1‑139 (MQ=255)
                                                                                                                    tttGCCCGGATGCACTTTTTTGATCTCTAGCGCAAACGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGAATCTTTCATTCGATTCATACGTCATCAACTGGTTTATATACTAAAGAGGTGATATTTGATCAcc  >  7:146139/1‑139 (MQ=255)
                                                                                                                          |                                                                                                                                    
ATTACTGTGGCTGCTTTTTTTGCACCGCTAATAAACGAATCAGCACGCTCGTCTATCACACCTAATTTTAATGCGGCTACTACCTGTTCATCAGAGATGTCGATTCCGATAATCCGTTTGCCCGGATGCTCTTTTTTGATCGCTAGCGCAATCGAACCGCCTATCAATCCGAGACCGGCGAGCAATATTGTATCTTTCATTCGATTCATACGTCATCACCTGGTTTATAGACTAAAGAGGTGATATTTGATCACC  >  CP006881/2188019‑2188273

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: