Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 986,184 C→A 25.1% intergenic (+50/‑200) U712_04810 → / → U712_04815 putative oxidoreductase yhdF/putative amino acid permease YhdG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881986,1840CA25.1% 146.6 / 6.2 74intergenic (+50/‑200)U712_04810/U712_04815putative oxidoreductase yhdF/putative amino acid permease YhdG
Reads supporting (aligned to +/- strand):  ref base C (12/41);  new base A (14/4);  total (29/45)
Fisher's exact test for biased strand distribution p-value = 5.00e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGACCATTG  >  CP006881/986058‑986318
                                                                                                                                 |                                                                                                                                      
tgctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAAt                                                                                                                            >  1:35225/1‑139 (MQ=255)
  ctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTTATTTAAACTTATAAAACAAAAAAAACCAAAATAGCTTGGTTTTTGTTTGGTGAAAAAATTAAAAAAAAAAAAatt                                                                                                                          >  7:96786/1‑139 (MQ=255)
  ctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTTTTTAAACGTAATCAACGAAAAACCAGCTCAATAGCTGGTTTTTTTTGTGGTGTAAAATTTAAAAAAAAAAAAatt                                                                                                                          >  3:311092/1‑139 (MQ=255)
     gCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTTTTTAAAAGTAATAAAACAAAAAACCACCCAAAAGCTGGGTTTTTTTTGTGGTGAAAAATTTAAAAATAAAAAAttttt                                                                                                                        >  2:239075/1‑137 (MQ=255)
     gCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTGTTATTTAAACGTAATAAACCAAAAAACACCGCAAGAGCGGGGTTTTTGTGTGGTGAAAAATTTAAAAAAAAAAAAtttttg                                                                                                                       >  3:120355/1‑136 (MQ=255)
        tCTGACGAATCTTCCTATATGACATGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGATGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGAt                                                                                                                    <  1:80838/139‑1 (MQ=255)
          tGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTGTTTTTTAAACGGTATAAAAGAAAAAACCCCACAAAAGGTGGGTTTTTTGTTGGGGTAAAAATTAAAAATAAAAAATTTTTGGAttt                                                                                                                  >  2:90479/1‑139 (MQ=255)
            aCGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAACCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGAtttt                                                                                                                 <  3:125685/138‑1 (MQ=255)
                         tatGACAGGGCAGACCATTCATGTGAATGGCGGCCGTGTTATTTGAACCTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTTTGGTGAAAAATTTAAAAATAAAAAATTTTTCGATTTTTATAACAAtttttt                                                                                                   >  3:381765/1‑137 (MQ=255)
                           tGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTTTTTAAACGTAATCAACGAAAAAACAGCTCAAGAGCTGGGTTTTTTTGTGGTGAAAAATTTAAAAATAAAAAATTTTTGGATTTTTATTAAAAtttttttt                                                                                                 >  8:52520/1‑139 (MQ=255)
                           tGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTAAACGTAATAAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTTTCGATTTTTCTAAAAAATtttttt                                                                                                 >  6:156530/1‑139 (MQ=255)
                               aGGGCAGACCATTCATGTGAATGGCGGCCGTTTTTTTTAAACGTAATCAACCAAAAAACACCTCAAGAGCGGGTTTTTTTTGTTGTGTAAAAATTAAAAATAAAAAATTTTTTGATTTTTATAAAAAATTTTTTTTAAt                                                                                             >  4:165623/1‑139 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCGTTTTTTTTAAACGTAATAAACGAAAAACCAGATTAATAGCTGGTTTTTTGTTGGTTGAAAAATTTCAAAATAAAAAATATGTGGATTTTTATTACAATTTTGttttgtttt                                                                                           >  6:88376/1‑134 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATAAACGAAAAACCACCTCAAGAGCTGGTTTTTTTTTTGGTGAAAAATTTAAAAATAAAAAATTTTTCGATTTTTTTAAAAATTTTTTTTTAAttt                                                                                           >  5:351301/1‑139 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCGTTTTATTTAAACGTAATAAACCAAAAAACACCTCAAGAGCGGGGTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTTTAGATTTTTATAACAATTTTTTTTTAAttt                                                                                           >  2:120217/1‑139 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCGTGTTATTTTAACGTATTAAAAGAAAAAACACCACAACCGCGGGTTTTTTGTGGGGGGTTAAAATTAAAAATAAAAAATTTTTGGATTTTTTTAACAATTTTTTTTTGAttt                                                                                           >  5:283420/1‑139 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCGTGTTATTTGAACCGAAAAAACGAAAAACCAGCACCAGAGCGGGTTTTTTTTGTGGTGTAAAATTTCAAAATAAAAAATTTTTCGATTTTTCTTACAAATTTTTTTTCAttt                                                                                           >  6:161364/1‑139 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCGTGTTATTTCAACGTAATCAACGAAAAAACAGCTCAAGAGCTGGTTTTTTGTGTGGTGTAAAATTTCAAAATAAAAAAATTGTCGATTTTTCTAACAATTTTGTTTTGAttt                                                                                           >  7:316779/1‑139 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCGTGTTATTTAAACGTAAAAAACCAAAAACCAGCAACAGAGCGGGTTTTTTTTGTGGGTTAAAATTTTAAAATAAAAAATTTTTTGATTTTTTTAAAAATTTTTTTTTAAttt                                                                                           >  6:6693/1‑139 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCCTTTTATTTCAACCTAATAAACCAAAAAACACCTCAATAGCTGGTTTTTTTTTTGGTGAAAAAATTAAAAATAAAAAAATTGAGGATTATTATAAAAATTTTGTTTTCAatt                                                                                           >  8:250348/1‑139 (MQ=255)
                                   cAGACCATTCATGTGAATGGCGGCCGTTTTTTTTTAACCTAATAAAACAAAAACAAACGAAAGAGCTGGTGTTTTTTGTGGGGGTAAAATTTAAAAAAAAAAATTTTTGTGTTTTTTTAAAAAATTTTTTTTTAttttt                                                                                         >  5:312357/1‑139 (MQ=255)
                                    aGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTg                                                                                        <  6:38509/139‑1 (MQ=255)
                                    aGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTg                                                                                        <  2:35225/139‑1 (MQ=255)
                                    aGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTg                                                                                        <  7:43389/139‑1 (MQ=255)
                                    aGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTg                                                                                        <  5:156530/139‑1 (MQ=255)
                                    aGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTg                                                                                        <  8:359021/139‑1 (MQ=255)
                                    aGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTg                                                                                        <  3:68904/139‑1 (MQ=255)
                                      aCCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTgta                                                                                      <  1:307350/139‑3 (MQ=255)
                                      aCCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCa                                                                                      <  8:96786/139‑1 (MQ=255)
                                      aCCATTCATGTGAATGGCGGCCGTGTTTTTTTAAAGGAATAAACCAAAAAACAACACAAACGCGGGGTTTTTTTGTTGGTTTAAATTTAAAAAAAAAAAAATTTTAGATTTTTTTTAAAAATTTTTTTTTAttttttg                                                                                       >  4:209285/1‑136 (MQ=255)
                                      aCCATTCATGTGAATGGCGGCCGTGTTATTTAAACGTAATCAACGAAAAAACCCCGCAAGAGCGGGTTTTTTTTGTGGTGTAAAATTTTTAAAAAAAAAATTTGTAGAATTTTTTAACAATTTTTTTTTGAttttttg                                                                                       >  2:153524/1‑136 (MQ=255)
                                         aTTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtt                                                                                   <  7:131281/139‑1 (MQ=255)
                                         aTTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtt                                                                                   <  3:97803/139‑1 (MQ=255)
                                           tCATGTGAATGGCGGCCGTGTTATTTAAACATAAAAAAACAAAAAAAACCAGAAAAACAGGGTTTTTTTTTGGTGAAAAAATTAAAAAAAAAAAAATTTTGTATTTTTATAAAAAATTTTTTTTAATTTTTTAAAtttt                                                                                 >  4:283278/1‑139 (MQ=255)
                                            cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                 <  8:211386/138‑1 (MQ=255)
                                             aTGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtat                                                                               <  5:161364/139‑1 (MQ=255)
                                                tgAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATATCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                            <  1:243712/139‑1 (MQ=255)
                                                 gtattgcgggCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                            <  6:351301/133‑1 (MQ=255)
                                                  aaTGGCGGCCGTTTTATTTAAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTAAAAAAAAAAAATTTGTCGATTTTTATAAAAATTTTTTTTTAATTTTTTGAATTttatttat                                                                           >  2:6252/1‑138 (MQ=255)
                                                      gcggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                       <  8:292678/138‑1 (MQ=255)
                                                      gcggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                       <  3:362748/138‑1 (MQ=255)
                                                      gcggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                       <  6:126282/138‑1 (MQ=255)
                                                      gcggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAAATTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                       <  4:10055/138‑1 (MQ=255)
                                                       cggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATAtt                                                                     <  3:362444/139‑1 (MQ=255)
                                                         gcCGTTTTATTTGAACGTAATAAACGAAAAACCAGCTAAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTGGATTTTTCTAACAATTTTTTTTTCATTTTTGCAATTTTATATTTAATTTTGa                                                                   >  1:97151/1‑139 (MQ=255)
                                                                 tATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGga                                                             <  2:84908/138‑1 (MQ=255)
                                                                   tttCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGgaga                                                           <  1:239075/138‑1 (MQ=255)
                                                                   tttCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGgaga                                                           <  5:6693/138‑1 (MQ=255)
                                                                        aCGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGaa                                                     >  7:188183/1‑139 (MQ=255)
                                                                         cGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAAt                                                    >  2:256110/1‑139 (MQ=255)
                                                                               acAACGAAAAACCAGCTCAACAGCTGCATTTTTCAGTGGTGAAAAATATCAAAATAAAAAATATCTCGATTTTTCTAACAATTATGTTATCAATATTGCAATTTTAAATAAAATATTGAGAATCAGAAAAGAATGATTc                                               >  8:47536/2‑139 (MQ=255)
                                                                                     aagaaCCAGCTCAATAGCTGGTTTTTTGTTTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTATAAAAATTTTTTTTTATTTTTTGCATTTTTATATATTATATTGAGTAGGAGAAATGAATGATTCGCGAcc                                         <  7:133713/136‑1 (MQ=255)
                                                                                     aaacaCCAGCTCCAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                         <  8:316779/139‑1 (MQ=255)
                                                                                     aaaaaCCCGCTCAAGAGGTGGTTTTTTGTGTGGTTAAAAATTTAAAAATAAAAAATTTGTCGATTTTTCTAAAAATTTTTTTTTAATTTTTGCAATTTTATATATAATATTGAGTAGGAGAACGGAATGATTCGCGAcc                                         <  3:209285/139‑1 (MQ=255)
                                                                                     aaaaaCCCGCTAAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTTTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                         <  6:162626/139‑1 (MQ=255)
                                                                                     aaaaaCCAGCTCAATAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                         <  4:120355/139‑1 (MQ=255)
                                                                                     aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                         <  4:311092/139‑1 (MQ=255)
                                                                                     aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                         <  3:165623/139‑1 (MQ=255)
                                                                                     aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGAAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                         <  7:52520/139‑1 (MQ=255)
                                                                                     aaaaaCCAGCTCAAGAGATGGTTTTTTGTGTGGTGAAAAATATCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTAATTTTTGCAATTTTATATATAATATTCAGTAGGAGAAAGGAATGATTCGCGAcc                                         <  4:307102/139‑1 (MQ=255)
                                                                                     aaaaaCCAGCTAAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTAAAAATAAAAATTTTGTCTATTTTTCTAACATTTTTTTTTTCTTTTTTGAAATTTTATATATAATTTTGTGTAGGAGAAAGGAATGATTCGCGAcc                                         <  3:17888/139‑1 (MQ=255)
                                                                                     aaaaaCCAGCTAAAGAGATGGTTTTTTTTGTGGTTAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAAAATTTTTGTTTTCATTTTTTCAATTTTATATATAATATTGAGTAGTATAAAGGATTGATTCGCGAcc                                         <  3:163299/139‑1 (MQ=255)
                                                                                     aaaaaCAATATCAATAGCTGGTTTTTTGTGTGGTGAAAAATTTAAAAATAAAAAATTTTTCGATTTTTCTAAAAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGATTAGGAGAAAGTATTGATTCGCGAcc                                         <  7:250348/139‑1 (MQ=255)
                                                                                      aaaaCCAGCTCAAGAGCTGTTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                         <  1:120217/138‑1 (MQ=255)
                                                                                       aaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTTAAAAATTTAAAAAAAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGAAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAg                                       <  1:381347/139‑1 (MQ=255)
                                                                                          ccAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTAc                                    <  3:78156/139‑1 (MQ=255)
                                                                                                   gagCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTTAGTAGGAGAAAGGAATGATTAGCGACCAGTACAAAtatatc                           >  8:133713/1‑137 (MQ=255)
                                                                                                          tttttttGTGTGGTGAAAAATTTCAAAATAAAAAATTTTTCTATTTTTCTAACAATTTTGTTTTATTTTTTGCAATTTTATATATAATATTGAGTAGGAGAACGGAATGATTCGCGACCAGTACAAATATAGTAGaaa                      <  5:384930/137‑1 (MQ=255)
                                                                                                              tttGTGTTGTTAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCAAACAATTTTGTTTAAATTTTTGAAATTTAATATAAAATATTGAGTCGGACAAAGGAATGATGCGCGACCAGTACAAATATAGTAGaaaaaaa                  <  6:349473/138‑1 (MQ=255)
                                                                                                                  tgtgGTGAAAATTTTCAAAATAAAAAATTTGTCTTTTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATTTATAATATTGAGTAGTAGAAAGGAATGATTCGCGACCAGTACAAATATATTAGAAAAAAAAGaat             <  7:402048/139‑1 (MQ=255)
                                                                                                                    tggttAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGaata            <  1:6252/138‑1 (MQ=255)
                                                                                                                    tggtgAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTAATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGaata            <  1:256110/138‑1 (MQ=255)
                                                                                                                       tgAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGa        >  1:254670/1‑139 (MQ=255)
                                                                                                                             aTTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGACCATTg  <  1:90479/139‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                      
TGCTGGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGACCATTG  >  CP006881/986058‑986318

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: