Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,256,0690GC28.9% 77.5 / 10.7 47P133R (CCG→CGG) plsBglycerol‑3‑phosphate 1‑O‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (12/20);  new base C (13/0);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 1.15e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCC‑‑AGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTT  >  NC_000913/4255838‑4256318
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                
aaaCGGGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCC‑‑AGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCAc                                                                                                                                                                                                                                             <  1:101423/251‑1 (MQ=255)
                        cGTTGACGGGCAAAGGTTATTCGCCCCGCGCTCGCC‑‑AGTTCCTGAGCGATAGTTTCACCCGGGCCGTGTCCATCCGCCATACGGCGCAGCGAACCGAACGGCGAGAAACCCACAAAACTGTCGCGACCGGGCCACTGTACAGCTAAAAATTTCTGTACGCCGTTATGCATACGCAGCTGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAACCATCACCGACACTGGCACCATCTGCACATc                                                                                                                                                                                                                    <  1:98933/251‑1 (MQ=255)
                           tgacggtcaaaggacctccgcccctcgcgcgCC‑‑ATTTTCTAAACAATAGTTTTATCCGTGCCGTGTTCATCCGCCCTCCGGCCCAGCTAAACTGACGGCGAGAAACGCACAAAACTGTCGCGCCCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCTGCGCGCGACCAAACATCACCTACACTGGCACCATCTGCACATCCa                                                                                                                                                                                                                  <  1:296093/225‑1 (MQ=255)
                             gagtggtcaaggggcagcgtcgaCGCGCGCCCGAGTTTCAGAGCTATAGTTTTACCCGTGACGGTGTCACCCGCCATACGGCGCCGGGAAACTGACGGCGAGAAACGCAAAAAACTTTCGAGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCATTGCGCCTATTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAg                                                                                                                                                                                                                <  2:690454/229‑1 (MQ=255)
                                ggaaaaGTGCATCCGCCCCACGCGCGCC‑‑ATTTTCGGAGCTAGTGTTTTAGCCGTGCCGTGTCCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAACCGCACAAAACTGGCGCGCCCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCTTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTg                                                                                                                                                                                                           <  1:407459/248‑1 (MQ=255)
                                        ccaacccccccaggcgcccC‑‑ATTTTCAGAGCGAAAGTTTTACCGGTGCCGTTTCAACCCGACATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAACCTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTa                                                                                                                                                                                                   <  1:368279/238‑1 (MQ=255)
                                            atacgagccctccgcctatttttttGAGCGTTAGTTTTACCTGTGCCGTGTTCTTCCGCCAAACGGCGCAGGGAAACTGACGGGGAAAACCGCACAAACCTGTGGCGACCGAGCCACAGTCCAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACg                                                                                                                                                                                                 <  2:784337/227‑1 (MQ=255)
                                             cgcgCCCCGCGCGCA‑‑TGTTGCTGAGCGTTATTTTTATCCTTGCCGTGTTCATCCGCCATACGGCGCAGCTAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCCCTTCGCCTTTTTCACGCCCCGGCGCGCGACTAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTg                                                                                                                                                                                              <  1:515808/251‑1 (MQ=255)
                                                       gcgcC‑‑AGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCTTTTTTACGCCCCCGCGCGCGACCAAACAACACCCACACTGGCACCATCTGCACCTCCACATTTGGGGTGCTACGGTGCaaaatggg                                                                                                                                                                                     >  1:722575/1‑245 (MQ=255)
                                                               gTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTTACTTCGCCTTTTTTACCCCCCCGCGCGCGGCCCAAAAACACCGACACTGGCACCCTCTTCACATCCCGATTTGGGGTTCTTCGGTGCAAATCGTGATAGTCg                                                                                                                                                                              >  1:374587/1‑251 (MQ=255)
                                                                     gagcattttttttacccggcccgtgTGCACCCGCCATACGGCGCAGCAAAACTGAGGGCGAAAACCGCAAAAACCTGTCGCGCCCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGaa                                                                                                                                                                         <  2:437789/230‑1 (MQ=255)
                                                                        ggATAGTTTTATCCGTGCCGGGTACACCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACGGTCGCGCCCGAGCCACAGTACAGCGAAAAATTGCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTGCCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGc                                                                                                                                                                     <  1:101577/250‑1 (MQ=255)
                                                                        cGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCTTTTTCCCGCCCCCGCGCGCGGCCAAACATCACCCACCCTGGCACCATCTGCACCTCCAGATTTTGGTTTCTACCGTGCAAAACGAGGTAGTCGTGGAACAGc                                                                                                                                                                     >  1:8387/1‑251 (MQ=255)
                                                                          aTAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCTAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTTACTTCCCCTTTTTTACCCCCCGGGGCGCGGACAAACCTCACCGCACATGGCACCATCTGCACACCCACATTTGAGTTGCTACAGAGCACAACAAGATAGTCGGGGaaaagatt                                                                                                                                                                   >  1:529503/1‑245 (MQ=255)
                                                                           tAGTTTTATCCGTGACGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTCCAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCCCATCTCCTTTTTTAACGCCCAGCGCGCGACCAACAATAAACGCACCTTGCACCATCTGCACAACCAAATATAGGCGGTCACGGTGCAAATAGAGCTAGTCGTGGAACAGcat                                                                                                                                                                   >  1:132750/1‑248 (MQ=255)
                                                                            aGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTTACTTCGCCTTTTTCACGCCCCCGCGCGCGGCCCAACATCACCGCCACTGGCACCCACTGCACATCCAAATTTGGGTTGCTACGGGGCCAAACGCGATAGTCTTGGGACAGCTTa                                                                                                                                                                 >  1:156980/1‑250 (MQ=255)
                                                                              ttttATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCTAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCCCTTTCCCTTTTTTACGCCCCGGCCCGCGGCCAAACATCACCCCCACTGGCACCCTCTCCACCCCCCGAATTTGGTTGCTACGGGGCCAAACGCGGTAGTCGTGGGACAGCtttat                                                                                                                                                                >  1:4354/1‑247 (MQ=255)
                                                                              ttttATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCCCCTTTTTTCCGCCCCCGCGCGCGACCAAACATCACCGCACCTCGCCCCCTCTGCCCATCCAAGTTTTGGTTGCTACGGTTCAAATCGAGATAGTCGTGGAACAACTTAAt                                                                                                                                                               >  1:187127/1‑250 (MQ=255)
                                                                                   cccgggccgtgTCCATCCGCCATACGGCGCAGCAAAACTGACGGCGAGAACCGCAAAAAACTGTCGCGCCCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGTTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGActc                                                                                                                                                         <  1:709955/246‑1 (MQ=255)
                                                                                    cccgtgccggtgtaaTCCGCCATACGGCCCGGCGAAACTAACGGCGAGAAACGCCCAAAACTGACGCGACCTAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGCCCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGActc                                                                                                                                                         <  2:292275/237‑1 (MQ=255)
                                                                                     cgtgccgtgTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTTACCCCCCGCCGCCCGGCCAAACATCCCCCACACTGGCCCCCTCTTCACAACCAGAATTTGGTTGCTCCGGTGCAAAACGAGATAGTCGTGGAACAGCTTAATAAACtctt                                                                                                                                                       >  1:236380/1‑251 (MQ=255)
                                                                                         ccgtgTCCACCCGCCATACGGCGCGGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTGTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCttt                                                                                                                                                   <  1:257611/251‑1 (MQ=255)
                                                                                            tgtTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTTACGCCCCCGCGCGCGACCAAACAACACCGACACCGGGACCATCTGCACATCCAGATTTTGGTTGCTACGGGGCAAAACGAGAAAGTCGGGGAAAAGCTTAATAGACTCTTCttttgg                                                                                                                                                >  1:276080/1‑248 (MQ=255)
                                                                                            gggttacaccccccATGCGGCGCAGGAAATCGACGGGCGAGAACCGCACAATATTGTCGGGCCCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTGAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTc                                                                                                                                                  <  2:761695/239‑1 (MQ=255)
                                                                                                 cATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGGGCGGGACCAAACATCACCGACACTGGCACCATTTGCACATCCAGATTTTGGGTGCTACCGAGCAAAACCAGAAAGTCGTGGAACAGCCCAATAAACTCCTCTTTCGGCggg                                                                                                                                            >  1:315756/1‑249 (MQ=255)
                                                                                                            aCGGCGCAGCGAAACGGACGGCGAGAAACGCACAAAACGGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCATGTAATAGGTg                                                                                                                                  <  1:105821/250‑1 (MQ=255)
                                                                                                                cgcAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCATTTTTACGCCCCCGGGCGCGACCAAACATCAACCACCCCGGCACCATATGACCAGCCAAGATTTGGTTGGTACGGTGCAAATCTGGCTAGTAGAGGAACAAGTTAATAAACTCTTATTTCGGCATGGAATAGGGGaaaac                                                                                                                            >  1:105513/1‑248 (MQ=255)
                                                                                                                        aaCTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGcc                                                                                                                    <  1:188203/251‑1 (MQ=255)
                                                                                                                            gACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCTCCCTTTTCACGCCCCGGCGCGCGGCCAAACATCACCCACACCGGCACCCTCTCCCCATCCAGATTTGGGGTGCTACGGGGCAAATCGCGAAAGGCGTGGGACAGCTTAATCGACTCTTCTTTTGGCGGGTAATAGGTGAACACACGCcgcccgcc                                                                                                                >  1:697427/1‑251 (MQ=255)
                                                                                                                                    agaAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGATTCACTTCCCCTTTTTTACGCCCCCGCGCGCGACCCAAACAACCCCACACTGGCACCCTCTGCACACCCAGCTTTGGGGTGCTACCGTGCAAATCCAGGCAGTCGTGGAAGCGCATAATAAAATCTTCTTTCGGCGGGTAATAGGGGAGACACCGCGGCCCGCCGagaatg                                                                                                         >  1:67425/1‑251 (MQ=255)
                                                                                                                                                        cgACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTTACGCCCCCGCGCGCGACCCAACAACACCCACACTGGCACCATCTGCACATCCCGATTTTGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCCGCGGGGAATAGGTGAACACACGCGGCCCGCCGCGGATGAACACAAAGCGCCGCAgg                                                                                       >  1:557552/1‑248 (MQ=255)
                                                                                                                                                                ccACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTTACGCCCCCGCGCGCGCCCAAAAAACACCGACACTGGCACCATCTGCACAACCCGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGCGGGACAGCTTAATAGACTCTTCTTTCCGCGTGGAATAGGTGAACACACGCGGCCCGCCGTGAAAGAAAACATAAATCGGCAGTAGCGTGCCgc                                                                             >  1:380089/1‑250 (MQ=255)
                                                                                                                                                                            gAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTTCCTTCTCCTTTTTCTCGCCCCCGCCCGCGCCCCAACATCACCCCCACTGGCACCATCTTCACATCCAGATTTTGGTTGCTCCGGGGCACATCGAGGTAACCGTGGAAAAGCTTTATAGAATTTTTTTTCGACGTGACATATGGGCACAAAAGCGGCCCGCCGCGAATGGAAACATAGCGCGCCAGAAGCGTGCCGGCGATTTTCAGCg                                                                 >  1:373130/1‑251 (MQ=255)
                                                                                                                                                                                      tGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTTCCCCCCCCGCGCGCGACCAAAAATCACCCACACTGGGAACCTCTGCACCACCCGATTTGGGTTGTTACGGGGCAAAACGAGATAGGCGGAGGACAGCTTTAAATGACCTTCTTTTTGCGGGTAATAGGTGAGAAAACGCGGCCCGCCCTGAAAGAAAACCTAACACGGGCGTAGCGTGCCgcccgtttctcccggcgctaaagg                                                       >  2:506157/1‑228 (MQ=255)
                                                                                                                                                                                           gCCGTTAAGCATACGCAGCGGCGGGTTCAATTCCCCTTTTTCACGCCCCCGCGCGCGGCCCAACATCACCCACACTTGCACCCTCTGCCCCTCCAGATTTGGGTTGCTACGGGGCAAATCGAGATAGACGTGGAAAAGCTTAAAAGACTCTTCTTTCGGCGAGGTAAAGGTGAACACACGCGGCCCCCCGTGAATGAAAACATAAGGCGGCAGAAGCGGGCCGTCGATTTTCCGCGGCCCTAACggggc                                                    >  1:373824/1‑247 (MQ=255)
                                                                                                                                                                                             cgtaaagcacccccAGCGGCGGTTTCTCTTCGCTTTTTTCCCGCCCCGGCGCGCGACCAAACACCCCCGACCCTGCCCCCATCTGCACATCCAGATTTGGGTTGCTCCGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCTTGTAATAGGTGACCACCCGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGc                                                <  1:559344/238‑1 (MQ=255)
                                                                                                                                                                                               ataggcctccgcggcgtggggTTCACTTCCCCTTTTTCCCGCCCCGGCCCCCACCAAAACACCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAACCGAGAAAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGACCACACGGGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGc                                                <  1:80206/231‑1 (MQ=255)
                                                                                                                                                                                                  aGCATACGCAGGGGCGGTTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACAAAACATCCCCGACACTGGAACCATCTGCACATCCAGATTTGGGTTGCTACGGTGAAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACCCGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTc                                           <  1:506157/251‑1 (MQ=255)
                                                                                                                                                                                                                cggGTTCACTTAGCCGTTTTCACGCCCCGGCGCGCGACCACACATCACCGACACTGGCACCCTCTGCCCATCACGATTTGGGTTGCTACGGTGCAAATCGAGATCGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTACTAGGTGAACACACGCAGCCCGCAGTGAATGAACACATAGCGCGGCAGTCGGGTGCAGTCGAGTACCTGCGACCATACCGGGTCAGTCAAACCATGCGCCAGACACt                             >  2:449182/1‑251 (MQ=255)
                                                                                                                                                                                                                  tgTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGGCAGGCAAGTCATGTGCCAGacacggg                           >  1:693530/2‑248 (MQ=255)
                                                                                                                                                                                                                   gTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCg                                                                 >  1:742204/1‑212 (MQ=255)
                                                                                                                                                                                                                       ccctcccctttttcccccccccGGCGCGCAACAAAACACCCCCGACACTGACACCTTCGGCACTCCAAGATTGGTGTTGCTACGGTCAAAATCGAGATAGTCGTGAAACAGGTAAATAGCCTCTTCTTTCGGCGTGTAAGAGGTAACCACCCGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGgcgc                       <  1:744830/233‑1 (MQ=255)
                                                                                                                                                                                                                              cTTTTTCACGCCACGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGACATGTGCAAGACACTGGGCGCGCAACgac               >  2:804514/1‑249 (MQ=255)
                                                                                                                                                                                                                                        gCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGGAGCGTGCCGTCGCTTTACAGCGGCTCAAACGGGTCAGGCAAGTCATTTGCCAGACACTTGGCGCGCAACGTAAGCAAATCt      >  2:682383/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ccccGGCGCGCGAACAAACATCCCCGCCACTGGCACTATCTGCACATCCAGATTGGGGTTGCTACGGTGAAAACCGAGATAGTCGTGAAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACCCGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTg     <  1:13233/251‑1 (MQ=255)
                                                                                                                                                                                                                                          cccGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTACCAGACAATGGGCGCGCAACGCAAGCAAAACTGc    >  2:632852/1‑251 (MQ=255)
                                                                                                                                                                                                                                            cGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGGTTTCCAGCGGCCCTAACGGGTCAGGAAAGTCATGTGCCAGACCCTGGGCGCGCAACGTCAGCACATCTGCtt  >  2:386266/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                
AAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCC‑‑AGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTT  >  NC_000913/4255838‑4256318

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