Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 969,533 G→T 17.7% D509Y (GAT→TAT)  U712_04735 → Aerobic glycerol‑3‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881969,5330GT17.7% 185.1 / 4.4 90D509Y (GAT→TAT) U712_04735Aerobic glycerol‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (33/39);  new base T (15/3);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 7.09e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCTGAAAGATGAGGCCGCGAAACGCAACATTCCGGTTCATATTCTTGCTGAGGCAGAATACAGTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGAAAAACAAACATACAGAAAACCTCAACAAGCTTTTACACGATGCGGTCGTACCGCTTG  >  CP006881/969400‑969668
                                                                                                                                     |                                                                                                                                       
aGCGCTGAAAGATGAGGCCGCGAAACGCAACATTCCGGTTCATATTCTTGCTGAGGCAGAATACAGTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGatat                                                                                                                                     >  5:73399/1‑138 (MQ=255)
   gcTGAAAGATGAGGCCTCGAAACGCAACATTCCGGTTCATATTCTTGCTGAGGCATAATACAGTATAGAAGATGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATc                                                                                                                                    <  4:394409/136‑1 (MQ=255)
   gcTGAAAGATGAGGCCGCGAAACGCAACATTCGGTTTCATATTCTTGCTGAGGCATAATACAGTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAAt                                                                                                                                 <  1:83922/139‑1 (MQ=255)
    cTGAAAGATGAGGCCGCGAAACGCAACATTCCGGTTCATATTCTTGCTGAGGCAGAATACAGTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGATGTTTATTTTTTGATATCAAtt                                                                                                                                >  5:235313/1‑139 (MQ=255)
    cTGAAAGATGAGGCCGCGAAACGCAACATTCCGGTTCATATTCTTGCTGAGGCAGAATACAGTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAAtt                                                                                                                                >  7:193603/1‑139 (MQ=255)
       aaaGATGAGGCCGCGAAACGCAACATTCCGGTTCATATTCTTGCTGAGGCAGAATACAGTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATTTCAAtt                                                                                                                                >  4:366147/1‑136 (MQ=255)
       aaaGATGAGGCCGCGAAACGCAACATTCCGGTTCATATTCTTGCTGAGGCAGAATACAGTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAAtt                                                                                                                                >  4:373864/1‑136 (MQ=255)
       aaaGATGAGGCCGCGAAACGCAACATTCCGGTTCATATTCTTGCTGAGGCAGAATACAGTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTGTTTTTTTTGATATCAAtt                                                                                                                                >  4:365600/1‑136 (MQ=255)
         aGATGAGGCCGCGAAACGCAACATTCCGGTTCATATTCTTGCTGAGGCAGAATACAGTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTa                                                                                                                           <  2:130337/139‑1 (MQ=255)
           ttGAGTCCGCGAAACTCAAATTTACGGTTCTTATTTTTGCTTAGTAATAATACTGTATCCACGAAGCGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAg                                                                                                                         <  1:263785/138‑1 (MQ=255)
           aggAGGCCGCGAAACTCAACATTCCGGTTAATTTTCTGGCTGAGTCATAATACAGTATCGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAg                                                                                                                         <  2:156977/137‑1 (MQ=255)
           aTGAGGCCGCGAAAGGCAAAATTCCGGTTCTTATTCTTGCTGAGGAATAATACTTTATAGAAGAAGACATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGTACGGGACGTTTATTTTTTGATATCAATTGGGt                                                                                                                            <  3:121140/136‑1 (MQ=255)
           aTGAGGCCGCGAAACGCAACATTCCGGTTCATATTCTTGCTGAGGCAGAATACATTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAg                                                                                                                         <  8:193603/139‑1 (MQ=255)
           aTGAGGCCGCGAAACGCAACATTCCGGTTCATATTCTTGCTGAGGCAGAATACAGTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAg                                                                                                                         <  6:91123/139‑1 (MQ=255)
           aTGAGGCCGCGAAACGCAACATTCCGGTTCATATTCTTGCTGAGGCAGAATACAGTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAg                                                                                                                         <  1:32629/139‑1 (MQ=255)
                         cgaaaCATTCCGTTTAATATTCTTGCTGAGGCATAATAAAGTATACAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATg                                                                                                           <  6:73399/136‑1 (MQ=255)
                                   cGGTCCATATTCTTGCTCAGGCAGAATACAGAAAACAAGACGCGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGATGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTAtt                                                                                                   <  6:235313/137‑1 (MQ=255)
                                          tatTCTTGCTGAGGCAGAATACAGTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTGTATTTTTTTTTATCATATGGGTAATAACATATAAAGAATCCGGTATTGATTTTAt                                                                                           >  1:371956/1‑138 (MQ=255)
                                            tacTTGCTGCGGCAGAATACAGTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATgag                                                                                        <  6:193682/137‑1 (MQ=255)
                                                tGCTGAGGCAGAATACAGTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTGTATTTTTTGTTATTAATTGGGTTAGAACATATAAAAATGCCGTTTTTTTTTTTATggtcgag                                                                                    >  1:213886/1‑139 (MQ=255)
                                                 gCTGAGTCACAATACAGTATACAAGAAGACATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATgagcgag                                                                                    <  7:367257/138‑1 (MQ=255)
                                                 gCTGAGGCAGAATACAGTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATgagcgag                                                                                    <  6:87059/138‑1 (MQ=255)
                                                 ccTGAGCCAGAATACAGTAACGACGAAGAGATGCCTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGGTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATgagc                                                                                       <  4:202422/134‑1 (MQ=255)
                                                   tGAGGCAGAATACAGTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTTTTTTTTGATTTCAATAGGGTAAGAACATATATATATGCCATTATTGATTTTATtagcgag                                                                                    >  3:394371/1‑136 (MQ=255)
                                                   tGAGGCAGAATACAGTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTTATATCAATTTGGTTATAAAATATAAATATTCCGTTATTTATTTTATtagcgag                                                                                    >  3:56692/1‑136 (MQ=255)
                                                   tGAGGCAGAATACAGTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTGTATTTTTTGATATCAATTGGGTTAGAACATATTAATATGCCGTTATTGATTTTATgatctat                                                                                    >  3:51049/1‑131 (MQ=255)
                                                        cAGAATACAGTATATATCAATACATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGAAGTTTATTTTTTGATATCAATTGTGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCa                                                                            <  1:136832/139‑1 (MQ=255)
                                                           aaTACAGTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATg                                                                         >  8:269717/1‑139 (MQ=255)
                                                           aaTACAGTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCAGCAGAACGGGACGGTTATTTTTTGAGATCAATTGGGTTAAAACATATAAAGATGCCGGTATTGATTTTTTGATCGAGCGATTCAAATg                                                                         >  8:332667/1‑139 (MQ=255)
                                                             tACAGTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTTATATAAATTGGGTTAGAACATATAAATATTCCGTTATTGATTTTTTGAGCGAGCGATTCCAATggg                                                                       >  8:371066/1‑139 (MQ=255)
                                                              aCAGTAAAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTGTATTTTTTGTTATAAATTTGGGTAGAAAAAATAAATATGCCCTTATTTAATTTTAGATCGAGCGATAACAAtc                                                                         >  4:184520/1‑135 (MQ=255)
                                                                 gTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATGAATTGGGTAAGAACATATAAAGATGCCGTGATTGATTTTATGAGCGCGCGATTCCCATGGGAt                                                                     >  7:354905/1‑137 (MQ=255)
                                                                  tataGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTTATATCAATTTGGTTAGAACATATAAATATTCCGTTTTTGATTTTATTAGGGATCGATTCCAATGGGATGAg                                                                  >  6:273327/1‑139 (MQ=255)
                                                                  tataGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTTGGTAAGAACATTTAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAg                                                                  >  5:69781/1‑139 (MQ=255)
                                                                  tataGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATAAAAAATGGGTAAGAACATATAAAGAAGCCGATATTGATTTTATGAGCGAGCGATTCCAATCGGATGAg                                                                  >  8:11713/1‑139 (MQ=255)
                                                                     agaagaagaGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTTTTTTTTGATATCAATTTGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGATCGAGCGATTCCCATTGGATGAGCAg                                                               >  8:32261/1‑139 (MQ=255)
                                                                     agaagaagaGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTGTATTTTTTTATATCAATTGGGTTAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAg                                                               >  2:34960/1‑139 (MQ=255)
                                                                       aagaagaGATGACTGCAACCCCTGGTGACTTCTTTGTCCGCAGAACGGGACGGTTTTTTGTTGTTAAGAAGTGGGTAAGAACAAAGAATCCTGGCGGTAGTGGTTTTATGCGCGCGGGCTTCCACTGGGATGAGcctgc                                                             >  2:18998/1‑135 (MQ=255)
                                                                       aagaagaGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAAt                                                                                                                                 <  8:403904/71‑1 (MQ=255)
                                                                       aagaagaGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAAt                                                                                                                                 >  7:403904/1‑71 (MQ=255)
                                                                        agaagaGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTTATATCAATTGGGTAAGAACATATAAATATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAg                                                               >  4:212794/1‑136 (MQ=255)
                                                                        agaagaGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATTTCAATTGGGTAATAAAATATTAAGAAGCCGTTATTTATTTTATTAGCCAGCGATTACAATGGGATGAGctg                                                               >  4:248209/1‑134 (MQ=255)
                                                                        agaagaGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTGTTTTTTTTGATTTCAATTGGGTTAGAACATATAAAGATGCCGTTATTGTTTTTATGAGCGAGCGATTACAATGGGGTGAGCAg                                                               >  4:416926/1‑136 (MQ=255)
                                                                         gaagaGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTGTATTTTTTTATATCAATTGGGGAAGAACATATTAAGATGGCGTTATTGTATTTATGGGCGAGGGGATGCAATGGGGTGTGCAGGCg                                                            >  5:288163/1‑138 (MQ=255)
                                                                         acccccaTGACTGCAACCCCTGCTGTCTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGa                                                           <  2:213886/133‑1 (MQ=255)
                                                                          aagaGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATTGGATGAGCAGGCGa                                                           >  8:334386/1‑138 (MQ=255)
                                                                          aagaGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAAATGGGGTAGAACATATAAAGATTCCGTTATTTATTTTATGATCGAGCGATTGCAATGGGATGAGCAGGc                                                             >  4:165576/1‑136 (MQ=255)
                                                                          aagaGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTGTATTTTTTGATATAAATTGGGTAAGAACATATAAAGATGCCGGTTTTGATTTTATGAGCGAGCGGTTTCAATTGGATGAGCAGGCGa                                                           >  5:252756/1‑138 (MQ=255)
                                                                            gagaTGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTTATATAAATTTGGTAATAACATATAAATAAGCCTTTATTGATTTTATGAGCGATCGAGTCAAATAGGATTTGCAGGCGaaa                                                         >  8:331651/1‑138 (MQ=255)
                                                                            gagaTGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTTAGAACATATTAAGATGGCGTTATTTATTTTATGAGCGAGCGTTTCCAATGGGATGAGCAGGCGaaa                                                         >  1:54360/1‑138 (MQ=255)
                                                                            gagaTGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACata                                                                                                                   >  5:335863/1‑80 (MQ=255)
                                                                            gagaTGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACata                                                                                                                   <  6:335863/80‑1 (MQ=255)
                                                                            gagaTGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCa                                                                            <  8:202596/119‑1 (MQ=255)
                                                                            gagaTGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCa                                                                            >  7:202596/1‑119 (MQ=255)
                                                                            gagaTGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGaaa                                                         >  5:365881/1‑138 (MQ=255)
                                                                            gagaTGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGaaa                                                         >  6:180415/1‑138 (MQ=255)
                                                                            gagaTGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTGTTTTTTTTTATATCAATTTGGGTAGAACATATAAAGATGCCGGTATTGATTTTATTAGCGAGCGATTACAATGGGATGAGCAGGCGaaa                                                         >  6:241602/1‑138 (MQ=255)
                                                                            gagaTGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTGTATTTTTTTATATAAAATGGGTAATAACATATAAACATGCCGGTTTTGATTTTATGAGCGAGCGAATACAATGGGATGAGCAGGCGaaa                                                         >  5:309427/1‑138 (MQ=255)
                                                                            gagaTGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTGTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGaaa                                                         >  2:294152/1‑138 (MQ=255)
                                                                            cagaTGACTGCAACCCATTCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTTATATGAATATGGTAAAAAAATATAAAGATGCCGTTTATGATTTTTTGACCGAGCGATTACAATGGGATGAGCAGGCGaaa                                                         >  7:115484/2‑138 (MQ=255)
                                                                               aTGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTTATATCAATTGGGTAAGAACATATAAAGATGCCGTTTTTGATTTTATGAGCGAGCGATTCCAAAGGGAAGAGCAGGCGaaa                                                         >  4:188000/1‑135 (MQ=255)
                                                                               aTGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTa                                                                                                     >  5:160357/1‑91 (MQ=255)
                                                                               aTGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTa                                                                                                     <  6:160357/91‑1 (MQ=255)
                                                                               aTGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGaaa                                                         >  4:416149/1‑135 (MQ=255)
                                                                               aTGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTGTATTTTTTTATATCAAATGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGaaa                                                         >  3:412343/1‑135 (MQ=255)
                                                                               aTGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTGTATTTTTTGATATAAATTGGGTAAGCACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGTTGAGCAGGCGaaa                                                         >  3:292006/1‑135 (MQ=255)
                                                                                tGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGaaaaa                                                       <  4:117261/136‑1 (MQ=255)
                                                                                  aCTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGTGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGaaaaa                                                       <  3:373864/134‑1 (MQ=255)
                                                                                   cTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGTGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGAAaaacaa                                                    <  3:365600/136‑1 (MQ=255)
                                                                                   cTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGAAAaacaaaca                                                 <  7:371066/139‑1 (MQ=255)
                                                                                   cTGCAACACCTGCTTACTTCTTTGTCCGCAGAACGGGTCGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGAAAaacaaaca                                                 <  7:331651/139‑1 (MQ=255)
                                                                                    tGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTTATATCaa                                                                                                                                  >  5:232524/1‑57 (MQ=255)
                                                                                    tGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCaa                                                                                                                                  <  6:232524/57‑1 (MQ=255)
                                                                                     acatccccTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTTATATCTATTGGGTAAGAACATATAAAGATGCGGTTTTTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGAAaaacaaac                                                  <  3:36767/132‑1 (MQ=255)
                                                                                      cAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGAAAaacaaaca                                                 <  3:416149/136‑1 (MQ=255)
                                                                                      cAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGAAAAACAAACATa                                               >  6:367858/1‑138 (MQ=255)
                                                                                                      cttTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGAAAAACAAACATACAGAAAACCTca                                   <  4:412343/134‑1 (MQ=255)
                                                                                                              gCAGAACGGGACGTTTATTTTTTTAAATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGAAAAACAAACATACAGAAAACCTCAACAAGCtttt                         <  4:19040/136‑1 (MQ=255)
                                                                                                                 gAACGTGAAGTTTATTTTTTTATATCATTTGGGTAAGAAAATATAAAGATGCCGTTATTGATTTTATGAGAGAGCGATTCCAATGGGATGAGCAGGCGAAAAACAAACATACTGAAAACCTCAACAAGCTTTTACACGa                   <  8:115484/139‑1 (MQ=255)
                                                                                                                   aCGTGACGTTTATTTTATGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGAAAAACAAACATACAGAAAACCTCAACAAGCTTTTACACGATg                 <  2:54360/139‑1 (MQ=255)
                                                                                                                   aCGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGAAAAACAAACATACAGAAAACCTCAACTAGCTTTTACACGATg                 <  7:11713/139‑1 (MQ=255)
                                                                                                                   aCGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGAAAAACAAACATACAGAAAACCTCAACAAGCTTTTACACGATg                 <  7:334386/139‑1 (MQ=255)
                                                                                                                     gggACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGAAAAACAAACATACAGAAAACCTCAACAAGCTTTTACACGATGCg               <  6:309427/139‑1 (MQ=255)
                                                                                                                         cGTTTATTTTTTGATATCAATTGGGTAAGAACATATAACGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGAAAAACAAACATACAGAAAACCTCAACAAGCTTTTACACGATGCGGTCg           <  6:288163/139‑1 (MQ=255)
                                                                                                                          gTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGAAAAACAAACATACAGAAAACCTCAACAAGCTTTTACACGATGCGGTCGt          <  6:168352/139‑1 (MQ=255)
                                                                                                                          gTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGAAAAACAAACATACAGAAAACCTCAACAAGCTTTTACACGATGCGGTCGt          <  5:180415/139‑1 (MQ=255)
                                                                                                                           tttATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTTATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGAAAAACAAACATACAGAAAACCTCAACAAGCTTTTACACGATGCGGt             <  3:237230/135‑1 (MQ=255)
                                                                                                                             tATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGAAAAACAAACATACAGAAAACCTCAACAAGCTTTTACACGATGCGGTCg           <  3:165576/135‑1 (MQ=255)
                                                                                                                                 ttttTATATAAATTTGGTAAGAACATATAAAGATGCCGTTATTTATTTTTTGATCGAGCGATTCCAATGGGATGAGCAGGCGAAAAACAAACATACAGAAAACCTCAACAAGCTTTTACACGAAGCGGTCGTACCGcgt   >  8:421711/1‑137 (MQ=255)
                                                                                                                                   ttGATATCAAATGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGAAAAACAAACATACAGAAAACCTCAACAAGCTTTTACACGATGCGGTCGTACCGCTTg  <  2:371956/138‑1 (MQ=255)
                                                                                                                                     gATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGAAAAACAAACATACAGAAAACCTCAACAAGCTTTTACACGATGCGGTCGTACCGCtt   <  3:366147/135‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                       
AGCGCTGAAAGATGAGGCCGCGAAACGCAACATTCCGGTTCATATTCTTGCTGAGGCAGAATACAGTATAGAAGAAGAGATGACTGCAACCCCTGCTGACTTCTTTGTCCGCAGAACGGGACGTTTATTTTTTGATATCAATTGGGTAAGAACATATAAAGATGCCGTTATTGATTTTATGAGCGAGCGATTCCAATGGGATGAGCAGGCGAAAAACAAACATACAGAAAACCTCAACAAGCTTTTACACGATGCGGTCGTACCGCTTG  >  CP006881/969400‑969668

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