Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,975,737 G→T 8.3% A198S (GCG→TCG)  wecE → dTDP‑4‑dehydro‑6‑deoxy‑D‑glucose transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,975,7370GT 8.3% 307.4 / 10.6 92A198S (GCG→TCG) wecEdTDP‑4‑dehydro‑6‑deoxy‑D‑glucose transaminase
Reads supporting (aligned to +/- strand):  ref base G (42/42);  new base T (3/5);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTG  >  NC_000913/3975589‑3975886
                                                                                                                                                    |                                                                                                                                                     
tgtgGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGc                                                                                                                                                      >  1:730404/1‑150 (MQ=255)
 gtgGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCg                                                                                                                                                     <  2:331060/150‑1 (MQ=255)
   ggTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGtt                                                                                                                                                   <  2:730404/150‑1 (MQ=255)
    gTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTa                                                                                                                                                  <  1:1516844/150‑1 (MQ=255)
      agaagaTGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAAt                                                                                                                                                <  1:511006/150‑1 (MQ=255)
       gaagaTGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCACGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATc                                                                                                                                               >  2:982076/1‑150 (MQ=255)
         agttGCCGCTCAGGGCGTGATTTCCACTTACAAAGGGCCTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGCTTTATCGTTAATcag                                                                                                                                             <  1:710569/147‑3 (MQ=255)
         agaTGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGTTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGCTTTATCGTTAATcag                                                                                                                                             <  2:721615/150‑3 (MQ=255)
          gaTGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGCTTTATCGTTAATcagc                                                                                                                                            <  2:112567/150‑4 (MQ=255)
           aTGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGCTTTATCGTTAATcagcg                                                                                                                                           <  2:1224096/150‑5 (MQ=255)
            tGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACg                                                                                                                                          >  2:804591/1‑150 (MQ=255)
             gCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGa                                                                                                                                         >  2:1096977/1‑150 (MQ=255)
             gCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGa                                                                                                                                         >  2:971336/1‑150 (MQ=255)
               cGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGc                                                                                                                                       <  1:980754/150‑1 (MQ=255)
                  tCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCga                                                                                                                                    <  2:1547457/150‑1 (MQ=255)
                   cAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCgag                                                                                                                                   >  2:1280556/1‑150 (MQ=255)
                      ggCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGatc                                                                                                                                <  1:73735/150‑1 (MQ=255)
                         gTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGatcatc                                                                                                                             >  1:879886/1‑150 (MQ=255)
                         gTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGatcatc                                                                                                                             <  2:937457/150‑1 (MQ=255)
                           gATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCg                                                                                                                           >  2:1082048/1‑150 (MQ=255)
                             tGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTg                                                                                                                         >  1:910236/1‑150 (MQ=255)
                                ccACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGaaa                                                                                                                      >  1:1023848/1‑150 (MQ=255)
                                     tACAAAGGGCGTTCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACTCTGATTAACGATTTAGCTTTAATCGAACGAGCCGATATCATCCGTGAAAAGGGc                                                                                                                 <  2:440659/150‑1 (MQ=255)
                                     tACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGc                                                                                                                 <  2:422832/150‑1 (MQ=255)
                                     tACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGc                                                                                                                 <  1:616376/150‑1 (MQ=255)
                                      aCAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCa                                                                                                                <  2:1065984/150‑1 (MQ=255)
                                      aCAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCa                                                                                                                <  1:841550/150‑1 (MQ=255)
                                       cAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCAc                                                                                                               >  2:1195867/1‑150 (MQ=255)
                                       cAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCAc                                                                                                               <  1:669183/150‑1 (MQ=255)
                                             gCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCg                                                                                                         <  2:450600/150‑1 (MQ=255)
                                             gCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCg                                                                                                         <  1:971336/150‑1 (MQ=255)
                                             gCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCg                                                                                                         <  1:1096977/150‑1 (MQ=255)
                                              cGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGc                                                                                                        <  2:1248022/150‑1 (MQ=255)
                                                 gCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGc                                                                                                     >  2:787090/1‑150 (MQ=255)
                                                 gCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGc                                                                                                     >  2:905018/1‑150 (MQ=255)
                                                  cACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGcc                                                                                                    <  2:1278707/150‑1 (MQ=255)
                                                    cTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAg                                                                                                  >  1:1264786/1‑150 (MQ=255)
                                                    cTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAg                                                                                                  >  2:895835/1‑150 (MQ=255)
                                                         aaCCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGttctt                                                                                             >  1:1146049/1‑150 (MQ=255)
                                                                ggTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGt                                                                                      >  1:238659/1‑150 (MQ=255)
                                                                 gTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTc                                                                                     <  2:1264786/150‑1 (MQ=255)
                                                                    atatTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAgg                                                                                  >  1:109019/1‑150 (MQ=255)
                                                                      atTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTc                                                                                <  2:1146049/150‑1 (MQ=255)
                                                                      atTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTc                                                                                <  1:787090/150‑1 (MQ=255)
                                                                      atTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTc                                                                                <  2:637015/150‑1 (MQ=255)
                                                                          gctgctTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATa                                                                            <  2:1023848/150‑1 (MQ=255)
                                                                             gctTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAat                                                                         >  2:378122/1‑150 (MQ=255)
                                                                               tttAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAatat                                                                       <  1:1035477/150‑1 (MQ=255)
                                                                                ttAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAAtata                                                                      >  1:1105655/1‑150 (MQ=255)
                                                                                 tAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATAc                                                                     >  2:1631498/1‑150 (MQ=255)
                                                                                 tAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATAc                                                                     <  2:910236/150‑1 (MQ=255)
                                                                                       ccATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGgcg                                                               <  1:698181/150‑1 (MQ=255)
                                                                                        cATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGgcgc                                                              >  1:59792/1‑150 (MQ=255)
                                                                                          tGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGa                                                            >  1:460097/1‑150 (MQ=255)
                                                                                           gAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAGAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGat                                                           >  1:1006853/1‑150 (MQ=255)
                                                                                               ccAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTg                                                       >  2:1303524/1‑150 (MQ=255)
                                                                                               ccAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTg                                                       <  2:476581/150‑1 (MQ=255)
                                                                                                cAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTgg                                                      >  1:1582594/1‑150 (MQ=255)
                                                                                                   aaaCTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTc                                                   <  1:895835/150‑1 (MQ=255)
                                                                                                        acacGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCt                                              <  1:941605/150‑1 (MQ=255)
                                                                                                            ggcgggcgGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTAttt                                          >  1:225057/1‑150 (MQ=255)
                                                                                                            ggcgggcgGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTAttt                                          >  2:264400/1‑150 (MQ=255)
                                                                                                             gcgggcggTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAg                                                                                                      <  2:1147003/89‑1 (MQ=255)
                                                                                                             gcgggcggTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAg                                                                                                      >  1:1147003/1‑89 (MQ=255)
                                                                                                              cgggcggTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGa                                        <  2:1105655/150‑1 (MQ=255)
                                                                                                              cgggcggTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGa                                        <  1:1199739/150‑1 (MQ=255)
                                                                                                              cgggcggTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGa                                        >  1:1365790/1‑150 (MQ=255)
                                                                                                                ggcggTGAAGGCGGCGCGACGCTGGTTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATg                                      <  1:1060806/150‑1 (MQ=255)
                                                                                                                  cggTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTc                                    >  2:961091/1‑150 (MQ=255)
                                                                                                                    gTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCg                                  >  2:124372/1‑150 (MQ=255)
                                                                                                                    gTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCg                                  >  1:600050/1‑150 (MQ=255)
                                                                                                                     tGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAAACGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGa                                 <  1:1482958/150‑1 (MQ=255)
                                                                                                                         ggcggcGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTg                             >  1:1316545/1‑150 (MQ=255)
                                                                                                                         ggcggcGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTg                             >  1:1653638/1‑150 (MQ=255)
                                                                                                                          gcggcgCGACGCTGATTAACGATAAATCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTGGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGc                            <  1:1468099/150‑1 (MQ=255)
                                                                                                                           cggcgCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCa                           <  1:1497853/150‑1 (MQ=255)
                                                                                                                                cgACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTg                      <  1:905018/150‑1 (MQ=255)
                                                                                                                                 gACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGc                     <  1:1303524/150‑1 (MQ=255)
                                                                                                                                  aCGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCa                    <  1:1637771/150‑1 (MQ=255)
                                                                                                                                    ttcgATTAACGCTTTATCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATa                  >  2:1425802/4‑150 (MQ=255)
                                                                                                                                     tcgATTAACGCTTTATCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATAc                 >  2:885585/3‑150 (MQ=255)
                                                                                                                                     cTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATAc                 >  1:1113870/1‑150 (MQ=255)
                                                                                                                                      cgATTAACGCTTTATCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATAcc                >  1:1308747/2‑150 (MQ=255)
                                                                                                                                       gATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCt               >  1:1628877/1‑150 (MQ=255)
                                                                                                                                          tAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCtgtg            >  2:274892/1‑150 (MQ=255)
                                                                                                                                           aaCGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTa                                             <  2:1581795/116‑1 (MQ=255)
                                                                                                                                           aaCGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTa                                             >  1:1581795/1‑116 (MQ=255)
                                                                                                                                            aCGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTggg          <  1:378122/150‑1 (MQ=255)
                                                                                                                                            aCGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTggg          >  2:129476/1‑150 (MQ=255)
                                                                                                                                                tAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCa      >  1:681819/1‑150 (MQ=255)
                                                                                                                                                tAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCa      >  1:770843/1‑150 (MQ=255)
                                                                                                                                                tAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCa      <  1:124372/150‑1 (MQ=255)
                                                                                                                                                    gCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTg  >  1:635139/1‑150 (MQ=255)
                                                                                                                                                    gCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTg  <  1:958871/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTG  >  NC_000913/3975589‑3975886

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