Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,189,655 T→C L173L (CTA→CTG U712_10970 ← Tryptophan synthase alpha chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,189,6550TC98.7% 240.3 / ‑4.8 77L173L (CTA→CTGU712_10970Tryptophan synthase alpha chain
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (43/32);  total (43/34)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

ATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCTAGAGAAGATACACAGTAGACGAAGCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTTACCTCATGGTTTTTGCAATCCGCTTGAAGGCTGGCGCTTTCT  >  CP006881/2189527‑2189790
                                                                                                                                |                                                                                                                                       
atCTTTATGACCTGTTCATGGTTTGATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATa                                                                                                                               <  6:56339/139‑1 (MQ=255)
atCTTTATGACCTGTTCACGGTTTGATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATa                                                                                                                               <  7:443294/139‑1 (MQ=255)
atCTTTATGACCTGTTCACGGTTTGATATAACGAACCCTACAGCAACCGGAACAGTGCTGAGATTGTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACAGGTCACACCCAGAGAAGATa                                                                                                                               <  1:295546/139‑1 (MQ=255)
  cTTTATGACCTGTTCACGGTTTGATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATa                                                                                                                               <  3:174271/137‑1 (MQ=255)
  cTTTATGACCTGTTCACGGTTTGATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATa                                                                                                                               >  4:174271/1‑137 (MQ=255)
  cTTTATGACCTGTTCACGGTTTGATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATAGCGGACACCGCTCACACCCAGAGAAGATaca                                                                                                                             >  3:400628/1‑139 (MQ=255)
           cTGTTCACGGTTTGATATACCCAACCCTACACCCACCGGAACAGTGCCGAGATGCTTCACAGAAGGAATGAACGGGTACACGGATGAATGGCACTCTTTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACg                                                                                                                    <  6:14576/139‑1 (MQ=255)
            tGTTCACGGTTTGATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATa                                                                                                                               >  1:310461/1‑127 (MQ=255)
            tGTTCACGGTTTGATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATa                                                                                                                               <  2:310461/127‑1 (MQ=255)
                 acgGTTTGATATACCGAACCCTACAGCAACCGGAAGAGTGCTGAGATTGTTCTCAGTACGATTGAACGGGTACACGTATGAATTGAACTCATTGCGGAGACCGGTCACAGGAAGAGAAGATACACAGTAGACGAAcccc                                                                                                              <  7:289158/138‑1 (MQ=255)
                   ggTTTGATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGCCACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTc                                                                                                            >  7:251090/1‑139 (MQ=255)
                      ttGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCTAGAGAAGATACACAGTAGACGAAGCCCTCGGc                                                                                                         <  6:80128/139‑1 (MQ=255)
                       tGATATCCCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGc                                                                                                         <  7:32619/138‑1 (MQ=255)
                       tGATATACCGAACCCTACAGCAACCGGACCAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGc                                                                                                         <  6:282578/138‑1 (MQ=255)
                       tGATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGc                                                                                                         <  8:71465/138‑1 (MQ=255)
                       tGATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGc                                                                                                         <  1:231628/138‑1 (MQ=255)
                       tGATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGc                                                                                                         <  1:198319/138‑1 (MQ=255)
                       tGATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCt                                                                                                        >  5:259535/1‑139 (MQ=255)
                       tGATAAACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACCGTACGAATGACCGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGc                                                                                                         <  4:381088/138‑1 (MQ=255)
                        gATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTCCACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCtt                                                                                                       >  6:227369/1‑139 (MQ=255)
                        gATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCtt                                                                                                       >  7:126009/1‑139 (MQ=255)
                        gATATACCGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCtt                                                                                                       >  8:285306/1‑139 (MQ=255)
                               cGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACCCCGGTCACACCCAGAGAAGATACACAGTAGACGACCCCCTCGGCTTGTTCaa                                                                                                 >  1:160186/1‑138 (MQ=255)
                               cGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCaa                                                                                                 >  1:131987/1‑138 (MQ=255)
                               cGAACCCTACAGCAACCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCaa                                                                                                 >  4:391786/1‑138 (MQ=255)
                                             aCCGGAACAGTGCTGAGATTCTTCACAGTACGAATGCACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCCCACCCAGAGAAGATACACAGTAGACGAACCCCACGGCTAGTTCAATAATGGTTTTCAAAc                                                                                  >  5:152654/1‑139 (MQ=255)
                                             aCCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAgaga                                                                                                                                    <  3:30180/89‑1 (MQ=255)
                                             aCCGGAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAgaga                                                                                                                                    >  4:30180/1‑89 (MQ=255)
                                                ggAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATa                                                                                                                               >  4:294760/1‑91 (MQ=255)
                                                ggAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATa                                                                                                                               <  3:294760/91‑1 (MQ=255)
                                                ggAACAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAACCGGc                                                                               >  3:150293/1‑139 (MQ=255)
                                                   aCAGTGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCttt                                                                            >  6:364196/1‑139 (MQ=255)
                                                       tGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCt                                                                        >  4:90471/1‑139 (MQ=255)
                                                       tGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCt                                                                        >  4:303512/1‑139 (MQ=255)
                                                       tGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCt                                                                        >  7:7724/1‑139 (MQ=255)
                                                       tGCTGAGATTCTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCt                                                                        >  2:33223/1‑139 (MQ=255)
                                                               atcTTCACACTACGCATGAACGGGGACCCGGATGAATTGAACTCATTGCGGACACCGGGCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGGTCAATAATGGTTTTCAAACGGCTTTCGCGTGTCGgcg                                                                <  4:235366/137‑1 (MQ=255)
                                                                 cTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAACCGGCTTTCGCTTGTCGgcgc                                                               >  5:76160/1‑138 (MQ=255)
                                                                 cTTCACAGTACGAATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGgcgc                                                               >  5:33938/1‑138 (MQ=255)
                                                                        gtacgaaTGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAg                                                       <  8:126009/133‑1 (MQ=255)
                                                                         tacgaaTGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACCCAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCACACGGCTTTCGCTTGTCGGCGCAACTAAAg                                                       >  1:414863/6‑138 (MQ=255)
                                                                         tacgaaTGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGCCGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAg                                                       >  7:299851/6‑138 (MQ=255)
                                                                                 aaCGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAg                                                                                                                          >  4:396437/1‑63 (MQ=255)
                                                                                 aaCGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAg                                                                                                                          <  3:396437/63‑1 (MQ=255)
                                                                                 aaCGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACg                                              >  2:401316/1‑139 (MQ=255)
                                                                                 aaCGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACg                                              >  5:193147/1‑139 (MQ=255)
                                                                                  aggggTTCACGGATGAATTGAACTCATTGCGGAAACCGGTCACACCCAGATAAGATACACAGTAGTCGAACCCCTCGGCTTGTTAAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGt                                             <  4:400628/137‑1 (MQ=255)
                                                                                  aCGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGt                                             <  3:90471/139‑1 (MQ=255)
                                                                                   cGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCag                                                                                                                                       >  5:175654/1‑48 (MQ=255)
                                                                                   cGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCag                                                                                                                                       <  6:175654/48‑1 (MQ=255)
                                                                                   cGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGtc                                            >  6:233864/1‑138 (MQ=255)
                                                                                         cacGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTc                                       <  8:92306/138‑1 (MQ=255)
                                                                                          acGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCAt                                     <  6:259535/139‑1 (MQ=255)
                                                                                          acGGATGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCAt                                     <  7:285306/139‑1 (MQ=255)
                                                                                              aTGAATTGAACTCATTGGGGACACCGGTCACACCCAGAGAAGATAcaca                                                                                                                           <  2:287142/49‑1 (MQ=255)
                                                                                              aTGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATAcaca                                                                                                                           >  1:287142/1‑49 (MQ=255)
                                                                                              aTGAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATggct                                 >  2:114206/1‑137 (MQ=255)
                                                                                                gAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCATTATTGGGTTTCAAACGGCTTTGGCTTGTCGGCGCTCCTAAAGCATTATACGTCTCCTCATGGcttg                               >  2:416212/1‑138 (MQ=255)
                                                                                                gAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGt                                             >  2:423505/1‑125 (MQ=255)
                                                                                                gAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGcttt                               >  2:49695/1‑139 (MQ=255)
                                                                                                gAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGcttt                               >  3:280436/1‑139 (MQ=255)
                                                                                                gAATTGAACTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGAAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGt                                             <  1:423505/125‑1 (MQ=255)
                                                                                                      aaCTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCAt                                     <  5:321970/127‑1 (MQ=255)
                                                                                                      aaCTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCAt                                     >  6:321970/1‑127 (MQ=255)
                                                                                                      aaCTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACatt                         >  6:126814/1‑137 (MQ=255)
                                                                                                      aaCTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTACACGAACCCCTCTGCTTGTTCAATAATGGTTATCAAACGGCTTTCGCTTGTCGGCTCAACTAAAGAAATATACGTCACCTCATGTCTTTTACat                          >  7:228640/1‑137 (MQ=255)
                                                                                                       aCTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTc                        <  5:364196/139‑1 (MQ=255)
                                                                                                        cTCATTGCGGACACCGGTCACCCCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTACAAACGGCTATCGCATGTCGGCGCAACTAAAGAAATATACGTCACCTCATGCCTTTTACATTcc                       >  3:394411/1‑139 (MQ=255)
                                                                                                        cTCATTGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATTGTTTTCAAACGGCTTTCGCTTGTCGGTGCAACTAAAGAAATATACGTCACCTCATGGCATTTACATTcc                       >  5:227986/1‑139 (MQ=255)
                                                                                                             tGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTcatctt                   >  4:342461/1‑133 (MQ=255)
                                                                                                             tGCGGACACCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCtt                   >  4:374980/1‑138 (MQ=255)
                                                                                                                    acCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCTTGAAGGCTg           <  5:120515/139‑1 (MQ=255)
                                                                                                                    acCGGTCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCTTGAAGGCTg           <  3:374980/139‑1 (MQ=255)
                                                                                                                         tCACACCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCTTGAAGGCTgttgct      <  4:280436/139‑6 (MQ=255)
                                                                                                                          cacaCCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCTTGAAGGCTGTTGCtt     <  2:17270/139‑1 (MQ=255)
                                                                                                                           acacCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCTTGAAGGCTGTTGCttt    <  1:421212/139‑1 (MQ=255)
                                                                                                                             acCCAGAGAAGATACACAGTAGACGAACCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTCACCTCATGGCTTTTACATTCCTCTTGAAGGCTGTTGCTttct  >  4:419154/1‑139 (MQ=255)
                                                                                                                                |                                                                                                                                       
ATCTTCATGACCTGTTCGCGGTTTGATATACCGAACCCGACAGCAACCGGAACAGTGCTGAGATCCTTCACTGTCCGGATGAACGGGTACACGGATGAATTGAACTCATTGCGGACACCGGTCACACCTAGAGAAGATACACAGTAGACGAAGCCCTCGGCTTGTTCAATAATGGTTTTCAAACGGCTTTCGCTTGTCGGCGCAACTAAAGAAATATACGTTACCTCATGGTTTTTGCAATCCGCTTGAAGGCTGGCGCTTTCT  >  CP006881/2189527‑2189790

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: