Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,027,240 (A)7→6 intergenic (+12/+369) fre → / ← fadA flavin reductase/3‑ketoacyl‑CoA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,027,2340A.98.5% 273.1 / ‑3.1 67intergenic (+6/+375)fre/fadAflavin reductase/3‑ketoacyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (38/28);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGA  >  NC_000913/4027102‑4027368
                                                                                                                                    |                                                                                                                                      
gCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTGGGCGATGCGTTTGCATTTATCTGAGATATaaaaaccc                                                                                                                                 <  1:129864/140‑4 (MQ=255)
    gaagagCATGATATCTATATTGCCGGACGTTTTGAGATGGCAAAAATTGCCCGCGATTTGTTTTGAAGTGAGCGAAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGcc                                                                                                                            <  2:62315/138‑1 (MQ=255)
       gagCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCt                                                                                                                         >  2:315247/1‑140 (MQ=255)
         gCATGATATCTATATTGCCGGACGTTTTGAAATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGa                                                                                                                       <  2:744039/140‑1 (MQ=255)
          cATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGAc                                                                                                                      >  1:537945/1‑140 (MQ=255)
           aTGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACa                                                                                                                     <  2:471537/140‑1 (MQ=255)
              atatCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGc                                                                                                                  >  2:229636/1‑140 (MQ=255)
                    atatTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGgaag                                                                                                            >  1:73195/1‑140 (MQ=255)
                    atatTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGgaag                                                                                                            >  2:337031/1‑140 (MQ=255)
                    atatTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGgaag                                                                                                            <  1:423031/140‑1 (MQ=255)
                       ttGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAAc                                                                                                         >  2:1019260/1‑140 (MQ=255)
                         gCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACgg                                                                                                       >  2:333765/1‑140 (MQ=255)
                          ccGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTTTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGc                                                                                                      >  1:877465/1‑140 (MQ=255)
                          ccGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAACCCCGCCCCTGACAGGCGGGAAGAACGGc                                                                                                      <  2:1136833/140‑1 (MQ=255)
                          ccGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGc                                                                                                      >  1:234558/1‑140 (MQ=255)
                          ccGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGc                                                                                                      >  1:632449/1‑140 (MQ=255)
                          ccGGACGTTTTGAGATGGCGAAAATTGCCCGCGATATGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGc                                                                                                      <  1:627607/140‑1 (MQ=255)
                          ccGGACGTTTTGAGATGGCGAAAATTGCCCGCGATATGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGc                                                                                                      >  2:627607/1‑140 (MQ=255)
                           cGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCa                                                                                                     >  1:710765/1‑140 (MQ=255)
                             gACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAAc                                                                                                   >  1:727392/1‑140 (MQ=255)
                             gACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAAc                                                                                                   >  2:718031/1‑140 (MQ=255)
                             gACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAAc                                                                                                   >  1:910659/1‑140 (MQ=255)
                               cGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTa                                                                                                 >  2:95684/1‑140 (MQ=255)
                               cGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTa                                                                                                 <  2:710765/140‑1 (MQ=255)
                               cGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTa                                                                                                 >  1:226601/1‑140 (MQ=255)
                                 ttttGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTaaa                                                                                               <  1:480658/140‑1 (MQ=255)
                                 ttttGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTaaa                                                                                               <  1:333765/140‑1 (MQ=255)
                                 ttttGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTaaa                                                                                               <  1:1019260/140‑1 (MQ=255)
                                           gCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTg                                                                                            >  1:1029762/1‑133 (MQ=255)
                                           gCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTg                                                                                            <  2:1029762/133‑1 (MQ=255)
                                            cGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAAATAAACTGTTATTCAg                                                                                    >  1:605668/1‑140 (MQ=255)
                                             gAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGt                                                                                   <  1:95684/140‑1 (MQ=255)
                                             aaaaaTTCCCCGCGATTTTTTTTGCAGTGGGGAAAAGGCGGGGGAAGACCGCTTTTTTGGCGATGCTTTTCCATTTATCTGAAATAAAAAAACCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGt                                                                                   <  2:605668/139‑1 (MQ=255)
                                              aaaaTTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTTTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTg                                                                                  <  2:877465/140‑1 (MQ=255)
                                              aaaaTTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAACCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTg                                                                                  <  2:87788/140‑1 (MQ=255)
                                              aaaaTTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTg                                                                                  <  2:727392/140‑1 (MQ=255)
                                              aaaaTTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTg                                                                                  <  2:234558/140‑1 (MQ=255)
                                                aaTTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGc                                                                                <  1:1028795/140‑1 (MQ=255)
                                                   tGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTGGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCAtt                                                                             <  2:533776/140‑1 (MQ=255)
                                                     cccGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTa                                                                           >  1:861591/1‑140 (MQ=255)
                                                     cccGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTa                                                                           <  1:320552/140‑1 (MQ=255)
                                                        gcgATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGAt                                                                        <  2:632449/140‑1 (MQ=255)
                                                                      aGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTg                                                          >  1:1022370/1‑140 (MQ=255)
                                                                       gTGAGCGTAATGCGGGGGAAGATCCCCTGTTTGGCGATGCGTTTCCATTTATCTGAGATATAAAAACCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTgg                                                         <  2:706230/140‑1 (MQ=255)
                                                                          aGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAACCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCa                                                       <  1:811444/139‑1 (MQ=255)
                                                                          aGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGGAGAACGGCCACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGGATCTGGCaa                                                      >  2:990145/1‑140 (MQ=255)
                                                                          aGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGGAAAACGGCAACTAAACTGTTATTCAGTGGCCTTTAGATCTATGACGTATCTGGCaa                                                      >  1:258391/1‑140 (MQ=255)
                                                                          aGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCaa                                                      >  1:168407/1‑140 (MQ=255)
                                                                          aGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCaa                                                      >  2:172015/1‑140 (MQ=255)
                                                                          aGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCaa                                                      >  1:692644/1‑140 (MQ=255)
                                                                          aGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCa                                                       >  2:811444/1‑139 (MQ=255)
                                                                                    gcgGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCg                                            >  2:620230/1‑140 (MQ=255)
                                                                                    gcgGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGCAAGAACGGCCACTAAAATGTTGTTCAGGGGCATTTAGGTCCATTGCGGATCTGGCCAAACATTCCCg                                            >  2:1120996/1‑140 (MQ=255)
                                                                                                ccTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGt                                >  1:1147727/1‑140 (MQ=255)
                                                                                                ccTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGt                                >  2:82491/1‑140 (MQ=255)
                                                                                                 ttttttGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTc                               <  1:116876/137‑1 (MQ=255)
                                                                                                 cTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAACCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTc                               <  2:988996/140‑1 (MQ=255)
                                                                                                  tGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTcg                              >  1:54953/1‑140 (MQ=255)
                                                                                                     ttGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGt                           <  1:450917/140‑1 (MQ=255)
                                                                                                      tGGCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAAAGGCAAATAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTc                          >  1:27989/1‑140 (MQ=255)
                                                                                                       ggCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGt                           >  2:608686/1‑138 (MQ=255)
                                                                                                       ggCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGt                           <  1:608686/138‑1 (MQ=255)
                                                                                                        gCGATGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCtt                        >  2:1026052/1‑140 (MQ=255)
                                                                                                            tGCGTTTGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCg                    <  2:73195/140‑1 (MQ=255)
                                                                                                            tGCGTTTGCATTTATCTGAGATAAAAAAACCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCg                    <  2:905781/140‑1 (MQ=255)
                                                                                                                 ttGCATTTATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCt               <  2:27989/140‑1 (MQ=255)
                                                                                                                        tATCTGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATa        >  1:219522/1‑140 (MQ=255)
                                                                                                                            tGAGATATAAAAACCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTcc    <  1:455827/140‑1 (MQ=255)
                                                                                                                            tGAGATATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTcc    <  1:172015/140‑1 (MQ=255)
                                                                                                                              agaTATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGa  >  2:854371/1‑140 (MQ=255)
                                                                                                                              agaTATAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGa  >  1:642302/1‑140 (MQ=255)
                                                                                                                                    |                                                                                                                                      
GCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGA  >  NC_000913/4027102‑4027368

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